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In this article, we review the principles of macroscopic quantum electrodynamics and discuss 
a variety of applications of this theory to medium-assisted atom-field coupling and dispersion 
forces. The theory generalises the standard mode expansion of the electromagnetic fields in 
free space to allow for the presence of absorbing bodies. We show that macroscopic quan- 
tum electrodynamics provides the link between isolated atomic systems and magnetoelectric 
bodies, and serves as an important tool for the understanding of surface-assisted atomic relax- 
ation effects and the intimately connected position-dependent energy shifts which give rise to 
Casimir-Polder and van der Waals forces. 



PACS: 42.50.Nn, 42.50.Ct, 34.35.+a, 12.20.-m 

KEYWORDS: Macroscopic quantum electrodynamics, Van der Waals forces, Casimir forces, 
Spontaneous decay, Cavity QED, Atom-surface interactions, Purcell effect, Spin-flip rates, 
Molecular heating, Input-output relations 

Contents 

1 Introduction 3 

2 Elements of vacuum quantum electrodynamics 5 

2. 1 Quantisation of the electromagnetic field in free space 5 

2.1.1 Lagrangian formalism 5 

2.1.2 Maxwell's equations 7 

2.1.3 Lossless beam splitter n 

2.1.4 Casimir force between two perfectly conducting parallel plates 15 

2.2 Interaction of the quantised electromagnetic field with atoms 18 

2.2.1 Heisenberg equations of motion 22 

2.2.2 Spontaneous decay and Lamb shift 25 

2.2.3 Optical Bloch equations 27 

2.2.4 Jaynes-Cummings model [28 



i 

2 E-mail address: s.buhmann@imperial.ac.uk 



E-mail address: s.scheel@imperial.ac.uk 



1 



2 



S. Scheel and S.Y. Buhmann 



3 Macroscopic quantum electrodynamics 

3.1 Field quantisation in linear absorbing magnetoelectrics 

3.1.1 Huttner-Barnett model 

3.1.2 Langevin-noise approach 

3.1.3 Duality transformations 

3.2 Light propagation through absorbing dielectric devices 

3.3 Medium-assisted interaction of the quantised electromagnetic field with atoms 

3.4 Nonlinear quantum electrodynamics 



Atomic relaxation rates 

4. 1 Modified spontaneous decay and body-induced Lamb shift, local-field corrections 

4.2 Sum rules for the local density of states 

4.3 Heating of polar molecules 

4.4 Spin-flip rates [TO] 

Dispersion forces 74j 

5.1 Casimir forces 

5.2 Casimir-Polder forces 

5.2.1 Perturbation theory 

5.2.2 Atom in front of a plate 

5.2.3 Atom in front of a sphere 

5.3 Van der Waals forces 

5.3.1 Perturbation theory 

5.3.2 Two atoms inside a bulk medium 

5.3.3 Two atoms near a half space 

5.3.4 Two atoms near a sphere 

5.4 Relation between dispersion forces 

5.5 Thermal effects 



Cavity QED effects 

6. 1 Quantum-state extraction from a high-Q cavity 

6.2 Spherical microcavities 

6.2.1 Atom inside a spherical microcavity . . 

6.2.2 Atom outside a spherical microcavity . 



Ill 

111 
113 
113 
115 



7 Outlook 117 



A Dyadic Green functions 

A. 1 General properties 1 19l 

A. 2 Duality relations 

A. 3 Bulk material 

A. 4 Layered media: planar, cylindrical, spherical 

A.4. 1 Planarly layered media 

A.4.2 Cylindrically layered media 

A.4.3 Spherically layered media 



Macroscopic QED — concepts and applications 

A.5 Born series expansion 

A.6 Local-field corrected Green tensors 



References 



4 



S. Scheel and S.Y. Buhmann 



1 Introduction 

In this article we review the basic principles, latest developments and important applications of 
macroscopic quantum electrodynamics (QED). This theory extends the well-established quan- 
tum optics in free space (Sec. [2]i to include absorbing and dispersing magnetoelectric bodies in 
its Hamiltonian description. In that way, a connection is established between isolated atomic sys- 
tems (atoms, ions, molecules, Bose-Einstein condensates) and absorbing materials (dielectrics, 
metals, superconductors). This is achieved by means of a quantisation scheme for the medium- 
assisted electromagnetic fields (Sec. [3). 

We set the scene by reviewing the basic elements of quantum optics in free space (Sec. 
Beginning with the quantisation of the electromagnetic field in free space in a Lagrangian formal- 
ism (Sec. l2.1.TI ) and based on Maxwell's equations (Sec. 12.1.21) . we briefly review two important 
applications, the lossless beam splitter (Sec. 12.1.31 ) and the mode summation approach to Casimir 
forces (Sec. 12. 1.4b . We introduce minimal-coupling and multipolar-coupling schemes (Sec. I2.2I > 
and discuss important consequences of the quantised atom-field coupling, spontaneous decay 
and the Lamb shift (Sec. 12.2.21 ). For completeness, we briefly review optical Bloch equations 
(Sec. 12.2.31 ) and the Jaynes-Cummings model (Sec. 12.2.41 ). 

The main part of this review deals with macroscopic quantum electrodynamics (Sec. [3]) and 
its applications (Secs.|3H6|. Macroscopic quantum electrodynamics provides the foundations for 
investigations into quantum-mechanical effects related to the presence of magnetoelectric bodies 
or interfaces such as dispersion forces and medium-assisted atomic relaxation and heating rates. 
The quantisation scheme is based on an expansion of the electromagnetic field operators in terms 
of dyadic Green functions, the fundamental solutions to the Helmholtz equation (Sec. 13.1.21 ). We 
discuss the principles of coupling atoms to the medium-assisted electromagnetic field by means 
of the minimal-coupling and multipolar-coupling schemes (Sec. 13.31 ), the latter of which will be 
used extensively throughout the article. After deriving the basic relations, we first focus our at- 
tention on medium-assisted atomic relaxation rates (Sec. 2]). We present examples of modified 
spontaneous decay, near-field spin relaxation, heating and local-field corrections. As we fre- 
quently refer to a number of explicit formulas for multilayered media, we have collected some 
of the most important relations in the Appendix (App.lAl). 

The Kramers-Kronig relations provide a close connection between the relaxation rates (line 
widths) and the corresponding energy shifts (Lamb shifts). The Lamb shift already exists in 
free space where the bare atomic transition frequencies are modified due to the interaction with 
the quantum vacuum. Because the quantum vacuum, i.e. the electromagnetic field fluctuations, 
are altered due to the presence of magnetoelectric bodies, these energy shifts become position- 
dependent and hence lead to dispersion forces. We develop the theory of Casimir, Casimir-Polder 
(CP) and van der Waals (vdW) forces on the basis of those field fluctuations (Sec. [5]). Amongst 
other examples, we discuss under which circumstances the results based on mode summations 
and perfect boundary conditions (introduced in Sec. 12. 1.41 can be retrieved. 

Finally, we apply the theory of macroscopic quantum electrodynamics to strong atom-field 
coupling effects in microresonators (Sec. [6]). Here we discuss leaky optical cavities from a field- 
theoretic point of view (Sec. 16. U which provides insight into input-output coupling at semi- 
transparent cavity mirrors and generalises the Jaynes-Cummings model (Sec. 12.2.41 ). We further 
present an example for entanglement generation between two atoms that utilises surface-guided 
modes on a spherical microresonator (Sec. 16.2.21 ). 
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2 Elements of vacuum quantum electrodynamics 

Let us begin our investigations of quantum electrodynamics by revisiting some of the basics of 
this theory — QED in free space. There are essentially two ways of approaching the quantisation 
of the electromagnetic field. In the quantum field theory literature (see, e.g. [1]), the formal 
route is taken in which a Lagrangian is postulated that fulfils certain general requirements such 
as relativistic covariance (Sec. 12.1. 11 1. In order to stress the intimate connection with classical 
optics [2], we instead follow a second approach by sticking to classical Maxwell theory as long 
as possible before quantising (Sec. 12. 1.2V 

A simple application of the mode expansion that will be employed in this section is the de- 
scription of the lossless beam splitter (Sec. 12.1.3b which will be extended to lossy devices in 
Sec. 13.21 The mode expansion approach will also be used to discuss the Casimir force between 
two perfectly conducting plates (Sec. 12. 1.4V As we will later see in Sec. 15.11 this interpretation 
cannot be upheld if the rather severe approximation of perfectly conducting plates is being weak- 
ened. Instead, we will have to describe Casimir forces (and related forces) in terms of fluctuating 
dipole forces. 

Next, we consider the coupling of the quantised electromagnetic field to charged particles 
(Sec. l2.21 i. We will introduce the notion of Kramers-Kronig relations and discuss simple atom- 
field phenomena such as spontaneous decay, the Lamb shift, the optical Bloch equations and the 
Jaynes-Cummings model of cavity QED. These examples will play a major role in our subse- 
quent discussion of macroscopic QED (Sees. [4] and 0. 



2.1 Quantisation of the electromagnetic field in free space 

In this section, we briefly describe the theory of quantum electrodynamics in free space. We out- 
line both the usual Lagrangian formalism and the more ad hoc approach via Maxwell's equations 
that highlights the connections with classical optics. 



2.1.1 Lagrangian formalism 



Within the framework of U(l)-gauge theories, the coupling between the fermionic matter fields 
and the electromagnetic field is described by a gauge potential = ((j>/c, A) which is identified 
as the four- vector of scalar and vector potentials. In order to determine the dynamics of in 
the Lagrangian formalism, a Lorentz covariant combination in terms of derivatives of A^ has to 
be sought, the simplest of which is the combination 



L = J d A xC = -j- J d^xF^F^ = ^Jdtjd 3 



e E 2 (r,i)- — B 2 (r,t) 
Mo 



(1) 



with the (covariant) anti-symmetric tensor 

/ 

Ffj,v = d^A v — d v A^ = 






E x /c 


Ey/C 


E z /c 


E x /c 





-B, 


B y 


Ey/C 


B z 





—B x 


E z /c 


-By 


B x 






(2) 
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Recall that the contravariant components of a four-vector x^ can be derived from its covariant 
components x^ by contraction with the metric tensor = diag(l, — 1, — 1, —1), x^ = g^ v x v . 
The equations of motion that follow from the Lagrangian density (Q3, 

d^ v = , d^ vp °F pa = (3) 

\t lLV( " y \ completely anti-symmetric symbol], are equivalent to Maxwell's equations 

V-B(r.t) = 0, (4) 

VxE(r,i) = -B(r,t), (5) 

V-D(r,i) = 0, (6) 

VxH(r,i) = D(r,t). (7) 

They have to be supplemented with the free-space constitutive relations 

D(r,i) = e E(r,i), (8) 

H(r,i) - — B(r,i), (9) 
Mo 

that connect the dielectric displacement field D(r, i) with the electric field E(r, t) and the mag- 
netic field H(r, t) with the induction field B(r, t). 

Returning to the Lagrangian (Q]l, one constructs the canonical momentum to the four-potential 

as 

SL 

n M = -p . do) 

Its spatial components are proportional to the electric field, II = — eoE, whereas the component 
IIo vanishes due to the anti-symmetry of F M „. This means that there is no dynamical degree of 
freedom associated with the zero component of the momentum field. Hence, the dynamics of the 
electromagnetic field is constrained. Using the canonical momenta, one introduces a Hamiltonian 
by means of a Legendre transform as 



H=- I ill I </'•/■ 



£ E 2 (r, t) + —B 2 (r,t) 
Mo 



(ID 



An additional complication arises due to the gauge freedom of electrodynamics. From the 
definition of F^, Eq. (fj), it is clear that adding the four-divergence of an arbitrary scalar function 
A, 

A^A^ + d^A, (12) 

does not alter the equations of motion, i.e. Maxwell's equations. One is therefore free to choose 
a gauge function A that is best suited to simplify actual computations. Clearly, any physically 
observable quantities derived from the electromagnetic fields are independent of this choice of 
gauge. A particular choice that preserves the relativistic covariance of Maxwell's equations is 
the Lorentz gauge in which one imposes the constraint d^A^ = 0. In quantum optics, where 
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relativistic covariance is not needed because the external sources envisaged there rarely move 
with any appreciable speed, the Coulomb gauge is often chosen. Here, one sets 

= 0, V-A = 0, (13) 

which obviously breaks relativistic covariance. Hence, in free space there are only two indepen- 
dent components of the vector potential. The scalar potential is identically zero; this is actually 
a consequence of the requirement V • A = rather than a separate constraint. In the Coulomb 
gauge, the Poisson bracket between the dynamical variables and their respective canonical mo- 
menta reads 

{A(r,(),n(r',f)} = i 1 (r-r') (14) 

where 8 (r) denotes the transverse 6 function. With the relations E = —Il/eo and B = V x A, 
the Poisson bracket for these fields simply read 

{E(r, t), B(r', t)} = -- V x S ± (r - r') , (15) 

which serves as the fundamental relation between the electromagnetic fields. At this point, 
canonical field quantisation can be performed in the usual way by means of the correspondence 
principle. Upon quantisation, the Poisson bracket has to be replaced by (ih)^ 1 times the commu- 
tator and Hamilton's equations of motion have to replaced by Heisenberg's equations of motion. 

2.1.2 Maxwell's equations 

Instead of using the field-theoretic Lagrangian language, we adopt a slightly simpler approach 
to quantisation that keeps aspects of the classical theory for as long as possible. Maxwell's 
equations © and (0 can equivalently be expressed in terms of the vector potential [derivable 
from Eq. (0] 

E(r, t) = -A(r, t) , B(r, t) = V x A(r, t) . (16) 
The vector potential A(r, t) in the Coulomb gauge dT3l > obeys the wave equation 

AA(r,i)-iA(r,i)=0. (17) 
The solutions to Eq. (Tl\ can be found by separation of variables, i.e. we make the ansatz 

A(r,t) = 5^A A (r)« A (t) (18) 

A 

which amounts to a mode decomposition. The mode functions Aa (r) obey the Helmholtz equa- 
tion 

AA A (r) + — j-A A (r) = (19) 
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where we defined the separation constant as oj 2 for later convenience. One can read Eq. (fT9b as an 
eigenvalue equation for the Hermitian operator — A having eigenvalues oj x /c 2 and eigenvectors 
Ax(r). Because of the Hermiticity of the Laplace operator, the mode functions form a complete 
set of orthogonal functions, albeit strictly only in a distributional sense. Hence, 



J d 3 rAl(r)-Ay(v)=M x 6 xyi £ ^A A (r) ® A A (r') = 5 ± (r - r') 



(20) 



where A/a denotes a normalisation factor. 

The Helmholtz equation ( fT~9b is easily solved in cartesian coordinates. The solutions are 
plane waves A a (r) = (k)e* where the magnitude of the wavevector k obeys the dispersion 
relation k 2 = ui 2 /c 2 . For each wavevector k there are two orthogonal polarisations with unit 
vectors e CT (k) obeying e CT (k) • e CT ' (k) = 8 aa i and e CT (k) • k = 0. Hence, the sum over A has in 
fact the following meaning: 



5>E 



d 3 k 

(27r) 3 /2 



(21) 



In cylindrical or spherical coordinates the solutions to the scalar Helmholtz equation can be 
written in terms of cylindrical and spherical Bessel functions, respectively (see Appendix lA.4.21 
and lAA3"V 

The temporal part of the wave equation reduces to the differential equation of a harmonic 
oscillator, 



u x {t) + u 2 x u x (t) = 0, 



(22) 



with solutions u X (t) — e A u x . Combining spatial and temporal parts, we obtain the mode 
decomposition for the vector potential as (u> = kc) 



(7=1 



d 3 k 
(27r) 3 / 2 



e CT (k) 



-i(k-r— Ljt) 



(23) 



where we have explicitly taken care of the reality of the vector potential by imposing the condi- 
tion u ka (t) = u*_ ka (t). 

Writing the expressions dT6b for the electric field and the magnetic induction in terms of the 
vector potential ( |23l . 



E(r,t) 



B(r,t) 



(T=l 

2 



d 3 k 

(27t) 3 / 2 
d 3 k 



M kcr e 



e CT (k)cj 
[k x e CT (k)] 



:(k-i 



•jt) 



* -i(k-r-wt) 



*ka c 



the Hamiltonian (TTTb reads 



H = - 



<7,(T'—1 



(k-r— tut) 



(2tt) 3 / 2 

d 3 rd 3 kd 3 k' 

,,* „-i(k-r-wi 



: -(k-r— cut) 



-i(k-r— oj£) 



(24) 
(25) 



£o(e CT • e CT /)cjo;' H (k x e CT ) • (k' x e a >) 

Mo 



Mk'o-'e 



'.(k'-r— uj't) 



*k'<r' 



-i(k -r— u; t) 



(26) 
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Using the orthogonality of the polarisation vectors e a ■ e a i = S aa i as well as the relation (k x 
e a ) ■ (k x e a i ) — k 2 (e a ■ e a i ), and integrating over r and k' leaves us with 

2 

H = 2£o H / d 3 ku; 2 \u ka \ 2 . (27) 

cr=l ^ 

The complex-valued functions can then be split into their respective real and imaginary parts 

as 



<7kcr = V^o (uka + K a ) , pka = -iu^fe^ (u k( T - > (28) 
which finally yields the classical Hamiltonian in the form 

2 

H =\H f d*k(pi a +^ q L) ■ (29) 

In this way, we have converted the Hamiltonian ( TTTI ) of the classical electromagnetic field into an 
infinite sum of uncoupled harmonic oscillators with frequencies u) = kc. The functions q^ a and 
Pko- are thus analogous to the position and momentum of a classical particle of mass m attached 
to a spring with spring constant k = mu 2 . 

The conversion of a field Hamiltonian into a set of uncoupled harmonic oscillators is the 
essence of every free-field quantisation scheme. In cartesian coordinates it is equivalent to a de- 
composition into uncoupled Fourier modes (or alternatively into Bessel-Fourier modes if cylin- 
drical or spherical coordinates are used). Note that the introduction of mode functions with 
the completeness and orthogonality relations (|20T > circumvents the usual problem of having to 
perform field quantisation in a space of finite extent, followed by the limiting procedure to un- 
bounded space at the end of the calculation. The analogy with classical mechanics can be pushed 
even further by noting that the functions q^ a and p^„ obey the Poisson bracket relation 

{gw,Pk' CT '} = <Kk-k')<W ■ (30) 

Quantisation is then performed by regarding the classical c-number functions q^ a and ptu 
as operators in an abstract Hilbert space H, and by replacing the Poisson brackets ( |30l by the 
respective commutators times [3]: 

<?kcr l-> <7k CT , Pkcr^Pko-, [qka,Pk>cr>} = «M(k - k')<5 CTO -/ . (31) 
By returning to the complex amplitude functions, now with different normalisation factors, 



a a (k) 




which obey the commutation rules 

= 5(k - k')<W , (33) 



Mk),4,(k') 
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we can write the operator of the vector potential in the Schrodinger picture as 

k ^=tj (^/£ eCT[eik '^ (k)+e " ik ' rat(k)] • 

The plane-wave expansion ( f34b is a special case of the more general mode expansion 
A(r) = ]T[A A (r)a A + A A (r)a A " 



(34) 



(35) 



The amplitude operators a A and a x then obey the commutation rules 



a x ,a[ 



= S 



XX' 



Finally, by introducing the amplitude operators via Eq. (1321 ). the Hamiltonian 
into diagonal form 



- \ [ a\a x + - 



(36) 

is converted 
(37) 



where the second equality follows from application of the commutation relation (13 6t , The last 
term in ((37) is an infinite, but additive, constant, the quantum-mechanical ground-state energy. 

With the expansion (l35l l at hand, it is now straightforward to write down the mode expansion 
of the operators of the electric field and the magnetic induction as 



E( r ) = *S LUx [ A A( r )«A - A* (r)a A 

A 

B(r) - ^[VxA A (r)aj|VxAl(r)al 



(38) 
(39) 



Using these expressions, we arrive at the (equal-time) commutation relations for the electromag- 
netic field operators as 



E(r),B(r') 



:]T c,a {A A (r) • [V x A A (r')] + A A (r) ■ [V x A A (r')]} 



_!!vx5 1 (r-r') 



(40) 



where we have chosen a normalisation factor A/a = ^/(2£owa) as in Eq. d34l and used the 
orthogonality relation (1201 1. This commutator agrees with the canonical commutator implied 
by the Poisson bracket Sl5[ on imposing the correspondence principle. The commutation rule 
(l40l tells us that the quantised electromagnetic field is a bosonic vector field. Its elementary 
excitations, the photons, of polarisation a and wavevector k are annihilated and created by the 
amplitude operators and a*. 

Note that the operators of the electric field and the magnetic induction describe the electro- 
magnetic degrees of freedom alone. The operators of the dielectric displacement D(r) and the 
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Fig. 1. Simple one-dimensional model of a beam splitter of thickness d, consisting possibly of several layers 
of different materials. 

magnetic field H(r), which in free space are trivially connected to E(r) and B(r) via Eqs. ((HJ 
and (O, in general contain both electromagnetic as well as matter degrees of freedom. We will 
see later in the context of macroscopic quantum electrodynamics that the same commutation 
rules d4Qb can be upheld even in the presence of magnetoelectric matter. The proof of the validity 
of this commutation provides a cornerstone of macroscopic QED. 

2.1.3 Lossless beam splitter 

Many propagation problems involving classical as well as quantised light involves finding the 
eigenmodes of the geometric setup and expanding the electromagnetic fields in terms of those 
modes. The plane-wave expansion in empty space was the simplest case imaginable. Manipu- 
lating light using passive optical elements such as beam splitters, phase shifters or mirrors are 
classic examples of mode matching problems at interfaces between dielectric or metallic bodies 
and empty space that can be solved by mode expansion approaches. 

In most cases of interest, it is possible to restrict one's attention to a one-dimensional prop- 
agation problem by choosing a particular linear polarisation and considering one vector compo- 
nent of the electromagnetic field only. For example, let us consider light propagation along the 
x-direction in which case the electric-field operator turns into a scalar operator 

E{x)=i [ dkc\k\A(k,x)d{k) - i [ dk c\k\A* {k,x)a) (fc) . (41) 



The mode functions A(k, x) satisfy the one-dimensional Helmholtz equation 

d 2 

-r^A(k,x) + n 2 (x)k 2 A(k 1 x) = (42) 
dx A 



with a spatially dependent (real) refractive index n(x) = y/e(x). The simplest model of a loss- 
less beam splitter involves assuming a refractive index profile with piecewise constant refractive 
index (Fig.[T|i 
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Re 



-ik'X 



Akx 



rpf ^ikx 



R'e 



-d/2 d/2 



-d/2 d/2 



Fig. 2. A plane wave e* x with k > (left figure) or k < (right figure) impinges onto a beam splitter 
from the left or right, respectively, and splits into transmitted and reflected parts. 



The solutions to the Helmholtz equation ( |42T > with the refractive index profile ( |43T > are once 
again plane waves that can be constructed similar to the familiar quantum-mechanical problem 
of wave scattering at a potential barrier (Fig. If an incoming plane wave e tkx from the left 
(k > 0) impinges onto the barrier, it will split into a reflected wave R(ui)e~ lkx and a transmitted 
wave T(u))e l . Similarly, if a plane wave e %kx enters from the right (k < 0), it will split into a 
reflected wave R'(u})e~ lkx and a transmitted wave T'(oj)e lkx with as yet unspecified reflection 
and transmission coefficients R(u>), R'(uj), T(u>) and T'(u>), respectively. Hence, the mode 
functions A(k, x) can be written as 



A ^ x > = V dbc { I > i 1 fc < ' (45) 

where ^4 is a normalisation area. The transmission and reflection coefficients can be obtained 
by requiring continuity of the vector potential (the quantum-mechanical wave function) and its 
first derivative at the beam-splitter interfaces. This requirement is analogous to the well-known 
conditions of continuity in classical electromagnetism. The result can be found in textbooks (see, 
e.g., [2]) as 

T(lu) = 1 ~ r . . e < < n - 1 ^/ c , (46) 

R{lo) = -re- lduj/c + re mduj/c T(Lo), (47) 

where r = (n — l)/(n + l) = ^§^j- is the Fresnel reflection coefficient for p-polarised waves 
(see App. lA.4.i1 i. 

For a single dielectric plate, the transmission coefficients T(uS) and T'(uj) and the reflection 
coefficients R(w) and R'(ui) are identical. For multilayered dielectrics, this is not necessarily 
the case. There are, however, a number of physical principles that restrict the form of these 
coefficients. For example, Onsager reciprocity [4] requires the magnitudes of the transmission 
coefficients to be identical, |T(w)| = |T'(w)|. In Sec. |A. 11 we will give details how this can be 
seen from general properties of the dyadic Green tensor. Moreover, energy conservation (photon 
number conservation, probability conservation) dictates that the squared moduli of transmission 
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and reflection coefficients must add up to unity, 



\T(u)\ 2 + \R(uj)\ 2 = 1 



(48) 



The correctness of Eq. d48l ) can be immediately checked by applying Eqs. ( f46b and ( |47| >. 

We see from Fig.[T]that the electric field can be decomposed into its incoming and outgoing 
parts associated with the photonic amplitude operators a,(o;) and bi(ui), respectively. The total 
electric field is thus the sum of field components travelling to the right (k > 0) and to the left 
(k < 0), whose amplitudes transform as 



bi(uj) = T{uj)ai{uj) + R'{uj)a2{oj) , 
b 2 (uj) = R(u>)ai(u>) +T'(w)&2(w) . 



(49) 
(50) 



For the transformed amplitude operators to represent photons, they have to obey similar commu- 
tation rules as the untransformed operators, hence we must have 



bi{u))Mw') = 5 ii 5{u- t J) 



(51) 



Itfollowsthat|T(w)| 2 + |i?V)| 2 = \T'{uj)\ 2 + \R(lu)\ 2 = 1 andT(w)J?*(cj) +R'(lj)T'*(u) = 
0. These requirements can be fulfilled if we set T'(uj) = T*(uj) and R'(uj) = -R*{u). The 
transformation rules of the photonic amplitude operators can thus be combined to a matrix equa- 
tion of the form [a(cj) = [ai(cL>), a2(w)] T , h(uj) = [b\(uj), b2{oj)] T ] 



b(w) = T(uj) ■ a(w) 
where the transformation matrix 

T{lj) R{lj) 



-R*(lj) T*{uj) 



(52) 



(53) 



is defined up to a global phase. Because of its structure, T(uo) is a unitary matrix, and in par- 
ticular, T(lu) 6 SU(2) [5-10]. The unitarity of the transformation matrix reflects the energy 
conservation requirement. 

The input-output relations (l52l can be converted from a matrix relation to an operator equa- 
tion as 



where the unitary operator U is given by 



U = exp 



oo 

J ciu;[a t (w)] T ■ ■ a(o 



T(lo) = exp[-i*(w)] 



(54) 



(55) 



Using that operator, the quantum state of light impinging onto the beam splitter transforms as 

Qout = UQintft . (56) 
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This can be easily verified by noting that the expectation value of any operator that depends 
functionally on the amplitude operators a(u>) and aJ (ui) can be computed either by transforming 
the amplitude operators using the input-output relations d54l >. or by transforming the quantum 
state using Eq. (l56*l l. The input-output relations would then correspond to the Heisenberg picture, 
whereas the quantum-state transformation could be regarded as its Schrodinger picture equiva- 
lent. 

Note that, despite the fact that the matrix T(w) describes an SU(2) transformation, the uni- 
tary operator U in general does not. As the n-photon Fock space is the symmetric subspace 
of the n-fold tensor product of single-photon Hilbert spaces [11], the quantum-state transfor- 
mation (T56b can be regarded as a transformation according to a subgroup of SU(2n) where 
n is the total number of photons impinging onto the beam splitter. For example, in the basis 
{|0, 0), |1, 0), |0, 1), |2, 0) , 1 1 , 1) , |0, 2)} the unitary operator U has the matrix representation 



U = 



( 1 






















T 


-R* 
















R 
























T 2 


V2T*R* 


R* 2 













V2TR 


(\T\ 2 -\R\ 2 ) 


-y/2T*R* 




V o 








R 2 


-V2TR 


T* 2 


) 



(57) 



which is block-diagonal with respect to the Fock layers of total photon numbers (0,1,2). This 
again expresses photon-number conservation. The unitary matrix U thus has the structure of a 
direct product, 



U = U n . 



(58) 



The fact that the matrix transforming the quantum states acts on the symmetric subspace of a 
tensor product Hilbert space means that it can be constructed from permanents of the transmis- 
sion matrix T [12] that is responsible for the operator transformation ([52}. Let us define the set 
G n ,N of all non-decreasing integer sequences us as 



Gr, 



{u: : 1 < oji < . . . < uj n < N} . 



(59) 



Then the unitary transformation of an A^-mode Fock state \ni , . . . , tijv) with a total of n photons 
can be written as [ft = (l™ 1 , 2™ 2 , . . . , N nN )] [13,14] 



-1/2 



U\m,...,n N ) = I Y[rnl\ -^V&T[u\n]\ mi (u;),...,rn N (u)) (60) 



in which the rrii(u>) are the multiplicities of the occurrence of the value i in the non-decreasing 
integer sequence uj and = Yii TO i( w )- The notation T[a;|r2] thereby stands for the matrix 

whose row and column indices are drawn from the non-decreasing integer sequences u> and fl, 
respectively, and whose elements are taken from the transformation matrix T. For example, 



Macroscopic QED — concepts and applications 



15 



the matrix T[(l, 1)|(1,2)] contains the elements 
permanent of a matrix T which is defined as 



( 



Tn T12 
Tn T12 



) 



. The symbol per denotes the 



n 



p erT = E U T - 



(61) 



o-eSn »=i 



where is the (symmetric) group of permutations. 

For example, the probability amplitude of finding exactly one photon in each beam splitter 
output when feeding two photons into its input ports, is given by the permanent of the beam 
splitter transformation matrix itself, 



For a symmetric beam splitter with vanishing permanent, this probability is zero, and the Hong- 
Ou-Mandel quantum interference effect is observed [15]. The appearance of a matrix function 
such as the permanent in the context of single or networks of beam splitters is one particular 
example of the links that exist between quantum optics and matrix theory. 

For lossy beam splitters, neither the conservation law (08]) nor the direct product structure 
of the matrix representation U of the unitary transformation operator [Eq. ( 1571 )1 hold. They will 
have to be replaced by a generalised conservation law and a generalised unitary operator that 
include material absorption (see Sec. 13.21 ). 



One of the most intriguing consequences of quantising the electromagnetic field (or indeed any 
other field theory) is the existence of an infinite ground-state energy 



which is present even if no photon is. It might be argued that this vacuum energy would not be 
physically relevant as all photon energies can be referred to this base level. Because no absolute 
energy measurements can be done, only relative with respect to this ground-state energy, that en- 
ergy would not be measurable. This reasoning, however, is incorrect. To see why, it is instructive 
to look at the quantities the ground-state energy depends on. Inspection of Eq. (l63l reveals that 
is the mode structure itself that determines its magnitude. In other words, the base level from 
which photon energies are counted can be changed by altering the number and the structure of the 
allowed electromagnetic modes. One way to achieve that is to confine the electromagnetic field 
in a geometric structure with appropriate boundary conditions that limits the number of available 
modes, for example between two perfectly conducting parallel plates (see Fig. [3j. Because of the 
boundary conditions for the electromagnetic field at the plates surfaces, the modes perpendicular 
to the plates are discrete with wave numbers k z = rnr/d, n £ N. Hence, the lowest wave number 
that is supported in the interspace between the plates is oc 1 jd. After replacing the mode sum in 
Eq. d63l by an integral over wave vectors, we find that the ground-state energy scales as l/d 3 . 
Because the ground-state energy clearly decreases with decreasing plate separation d, there exists 



(1, 1\U\1, 1) = per T = TuT 22 + T 12 T 21 = \T\ 2 - \R\ 2 . 



(62) 



2.1.4 Casimir force between two perfectly conducting parallel plates 




(63) 



A 
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L 




Fig. 3. Electromagnetic modes confined between two perfectly conducting parallel plates of separation d 
lead to an attractive force. The figure on the right shows some of the allowed modes with the cavity walls. 



an attractive force between them. By dimensional arguments, this force per unit area L 2 must be 
proportional to 

F he 

^oc-^. (64) 

A more detailed calculation yields the correct numerical prefactor (we will follow the deriva- 
tion presented in Chap. 3.2.4 in Ref. [1]). First note that the ground-state energy in a box of 
volume L 2 d with L 3> d is 



£ (d) = l5> A = ^k A | = ^/ 



A A 



d 2 k\\ 



2 2\ 1/2' 



d 2 



(65) 



where we have used the fact that for k z > there are two possible polarisations a whereas there 
is only one independent polarisation if k z = 0. This expression is, of course, infinite. To render 
this expression finite, we subtract the contribution of free space, 





where the double counting of the polarisation state at n = does not influence the value of the 
integral. The ground-state energy per unit area L 2 is thus, using polar coordinates, 



E (d) - E {oo) _ he 
L 2 ~ 27 



3C: 




dkk\-+) \ k 2 + —^- dn\k 2 + —r-\.{61) 



This expression still seems to diverge for large wave numbers k and has to be regularised. For 
this purpose, we introduce a cut-off function f(k) such that 

J ^~ ' k <i & max (68) 

10, k 3> tax 
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Fig. 4. Smooth cut-off function f(k) that falls to zero for large enough wave numbers. 



where the cut-off wave number fc max could be chosen to be of the order of the inverse atomic 
size or, more appropriately, to correspond to a frequency larger than the plasma frequency of the 
material that the plates consist of (see Fig.[4]i. The physical background to the latter requirement 
is that the permittivity of a metal is for frequencies larger than the plasma frequency up well 
described by e(uj) «1- uip/uj 2 . Thus, for u> 3> wp even metals become transparent and fail to 
provide the required boundary conditions. 

After change of variables to u — d 2 k 2 /it, one obtains the convergent expression 



E (d) -Bb(oo) _ ftcTT 2 



L 2 



4d 3 



1 °° 7 

-F(0) + J2F(n)- J dnF(n) 



(69) 



with 



F(n) = I duy/ u + n 2 f Oifs/ u + n 2 /dj = / duy/uf (tt^/u/cC) 



(70) 



The expression in brackets in Eq. (1691 can be computed using the Euler-McLaurin resummation 
formula 



OO 

^(0) + £ F(n) - fdnF(n) = - £ ^-F^(0) 

n=l m=l ^ 



(71) 



where the Bk are the Bernoulli numbers. Since we have constructed the cut-off function such that 
/(0) = 1 and (0) = 0, the only non-zero contribution to Eq. {71} arises from F"'(0) = -4, 
together with B4 — —1/30. Hence, 



E (d)-E {oo) _ tt 2 He 
L 2 _ _ 720d3 



(72) 



which leads to a force per unit area as 
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geometry 


Casimir force 


Ref. 


planar 


n 2 hcL 2 
240d 4 


[16] 


cylinder 


0.02712ftcz 
i? 3 


[17] 


sphere 


0.0A5hc 
+ R> 


[18] 



Tab. 1. Casimir forces for parallel plates and cylindrical and spherical shells. 



which is precisely Eq. d64l i up to a numerical factor of 7r 2 /240 w 0.041. 

Similar calculations yield the Casimir forces for cylindrical and spherical shells as shown in 
Tab.Q] Note here that the Casimir forces in both planar and cylindrical geometries are attractive, 
whereas in case of a spherical shell it is repulsive. The latter result seems to contradict our 
intuition that a restriction of the number of modes always leads to an attractive force. In order 
to resolve this conundrum, one needs to look closer at the mode structure inside and outside a 
cylindrical or spherical shell. 

In the mode-summation approach that forms the basis of the calculations referred to above, 
one has to regularise the wave number integral at its upper limit by assuming a cut-off frequency 
above which the plates have to become transparent (for cylindrical and spherical shells one some- 
times assumes two half-cylinders or hemispheres whose separation provides the necessary reg- 
ularisation). This argument already suggests that the interpretation of the Casimir force as a 
mode restriction between perfect conductors cannot be upheld rigorously, and must be replaced 
by something that involves the dielectric properties of the plates. 

Let us interrupt the flow of the argument at this point and mention a classical analogue that 
can serve as an intuitive guidance to the problem of Casimir energies: the problem of determining 
altitudes on land. On literally all geophysical maps, the altitude of landmarks such as mountains, 
lakes, and human dwellings is given in terms of its altitude with respect to the average sea level. 
So one could say that the average sea level represents the level of the infinite ground-state energy. 
And in exactly the same way in which one is not interested in the altitude of a mountain with 
respect to the sea floor, we shall be content with measuring the photon energies from the infinite 
ground-state level. On the other hand, one might ask the question how the sea level can be 
properly defined given that there are tides, wind and waves that distort that level. As we will see 
later, it is exactly these fluctuations that are responsible for the Casimir effect in the quantum- 
mechanical setting. 

2.2 Interaction of the quantised electromagnetic field with atoms 

After we have determined how to quantise the electromagnetic field in free space, we will now 
couple external sources to the field and focus on the atomic degrees of freedom. For this purpose, 
let us begin again with classical Maxwell's equations which, in the presence of external sources, 
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read 

V-B(r,i) = 0, (74) 

VxE(r,i) = -B(r,t), (75) 

V-D(r,t) = p(r,t), (76) 

VxH(r,t) = j(r,t)+D(r,t). (77) 

The charge density p(r, t) and the current density j(r, t) fulfil the equation of continuity 

p(r,f)+V-j(r,i)=0 (78) 

which states that any change of the charge distribution within a region of space is accompanied 
by a flow of current across the boundary of that region. 

The charge density for an ensemble of point charges q a is given by 

p(T,t)=J2<l<*S[r-r a (t)] (79) 

a 

where r a (t) denotes their classical trajectory. From the continuity equation (l78l it then follows 
that the current density is 

j(r,t)=J2<io t r a (t)S[r-r a (t)]. (80) 

a 

In order to promote the charges to proper dynamical variables, we have to supplement Maxwell's 
equations with Newton's equations of motion for particles with mass m a , 

111 a Y a = Q a [E(r a ,t) + r a x B(r a ,t)] . (81) 

Introducing scalar and vector potentials as in free space, 

B(M)=VxA(r,i), E(r,t)=-A(r,i)-V0(r,t), (82) 

we obtain their respective wave equations in the Coulomb gauge V ■ A(r, t) = as 

A0(r,t) = --p(r,t), (83) 



AA(r,t)--L\(r,f) = -fx j(r,t) - £ O V0(r,i) 



c 2 



Equation ( |83] l is Poisson's equation with the solution 



4tt£ 7 ' |r-r'| 47re ^ |r - r a (t)\ 

where in the second equality we have used Eq. (|79l . The expression on the rhs of Eq. 
transverse current density which can be written in terms of the total current density as 



(84) 



(85) 



j>, ( )= jM )-v/, V ™ 



(86) 
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Fig. 5. Individual charged particles are combined to a coarse-grained system at position ya- 

using the continuity equation ( p78T >. 

The above equations of motion for the electromagnetic field and the charged particles can be 
derived from the classical Hamiltonian function 



// \ I il'-r 



e E 2 (r,t) + — B 2 (r,t) 
Mo 



Y] 7T— [Pa - q a A(r a )} 2 + V . 

^ 2m Q t-f 87re r Q - r Q * 

(87) 



in which the last term describes the Coulomb interaction between the charged particles. Note 
that the particle momentum p Q = m a r a + q a A(r a ) is different from the mechanical momentum 
m a r a due to the interaction with the electromagnetic field. 

The Hamiltonian d87b is referred to as the minimal-coupling Hamiltonian because the elec- 
tromagnetic field couples to the microscopic degrees of freedom of the individual charged parti- 
cles such as position and momenta. This microscopic description is often rather unwieldy, and 
an alternative approach in terms of global quantities is sought. A particularly important situa- 
tion arises if the individual charged particles constitute an ensemble of bound charges such as 
electrons and nuclei in an atom or a molecule. Let us introduce a coarse-grained charge dis- 
tribution p and current density j associated with that atomic system at centre-of-mass position 

Y A = J2a( m a/ m A)r a ( m A = E a m a ■ total mass ) (Fig.g), 



P(r, t) - Q^j S i r - r A(t)] , j - (^2 qc)j r A (t)6[r - r A (t)\ 



(88) 



Note that this charge density is zero for globally neutral systems. In order to make up for the 
difference between the actual charge distribution ( |79l ) and Eq. ( 1881 ), we define the microscopic 
polarisation field via the implicit relation 

V-P A (r,t)=-p(r,t)+p(r,t). (89) 

Adding this polarisation to the electric field, we define the modified displacement field as 



D(r,i) = £o E(r,t) + P(r,i) 



(90) 
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which obeys the modified Coulomb law 

V-D(r,f)=p(r,f). (91) 

Note that for globally neutral systems, the displacement field is a transverse vector field. From 
the implicit relation d89l . one can show that the polarisation can be written as [19] 



(92) 



(r a = r a — ta, relative particle coordinates). Analogously, the magnetisation field is introduced 
via the relation 

V x M A (r) - j(r) - j(r) - P A {v) (93) 
which can be written as [19] 



(94) 



As before, field quantisation is performed by replacing the relevant c-number quantities by 
Hilbert space operators and postulating their canonical commutation rules. In contrast to rel- 
ativistic quantum electrodynamics, the charged particles are not treated within second quanti- 
sation, i.e. in quantum optics electrons, atoms, molecules etc. cannot be created or annihilated. 
Instead, their quantum character is contained in their respective position and momenta, for which 
we postulate the canonical commutation rules 

[r„,pj] = iHSaa'I. (95) 

At the moment, it seems as if we have not achieved anything other than rewriting the charge 
density in terms of a new vector field that itself, by construction, depends on the original micro- 
scopic variables. To proceed, one either has to solve the microscopic dynamics explicitly which 
is only possible for sufficiently small systems, or one can invoke statistical arguments that relate 
the polarisation field causally to the electric field by means of a (in general nonlinear) response 
ansatz- The latter approach leads to the theory of macroscopic QED (Sec. [3]). 

A direct treatment of the microscopic dynamics can be considerably simplified by casting the 
atom-field interactions appearing in the minimal-coupling Hamiltonian (T87T > into its alternative 
multipolar form. To that end, we transform the dynamical variables by means of a Power-Zienau 
transformation O' = UOU^, where the unitary operator [20-22] 



U = cxp 



~ J d 3 rP A (r)-A(r) 



(96) 



depends on the polarisation d92| > and the vector potential (1351 ) of the electromagnetic field. Ex- 
pressing the Hamiltonian d87l i in terms of the transformed variables and applying a leading-order 
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expansion in terms of the relative particle coordinates, one obtains the electric -dipole Hamilto- 
nian for a neutral atomic system (cf. also Sec. 13.31) 



(97) 



Here, 

a. a 

is the electric dipole moment operator, and the transformed electric field 

E'(r)=E(r) + -Pi(r) (99) 

can be given in terms of the transformed bosonic operators a' x via an expansion analogous to 
Eq. ( 138) . The multipolar Hamiltonian is highly advantageous in comparison to the minimal 
coupling one (|87| | since the atom-field interaction H mt = — d • E(r^) consists of a single term. 
For this reason, we will exclusively employ it throughout the remainder of this section and drop 
the primes distinguishing the multipolar variables. 

2.2.1 Heisenberg equations of motion 

The electric -dipole interaction Hamiltonian H[ nt = — d-E(r^), can again be expanded in modes 
according to the description given above. In particular, the operator of the electric field strength 
is given by Eq. d38t . For the atomic system we choose to expand its free Hamiltonian and the 
dipole moment in terms of its energy eigenbasis |n). The atomic flip operators will be denoted 
by A mn = \m)(n\, and they obey the commutation rule 

A m n,Au = A m i5 n k — Akn§lm ■ (100) 

With these preparations, the electric-dipole Hamiltonian d97l > takes the form 

H = huj\ (a\a\ + 2 ) + X! ^nAm - i ^ uj mn A mn d mn ■ A(ta) 

A 71 m,n 

= ^ ^ ( a\ax + - J + ^ hw n A nn 



uj mn A mn d mn ■ Ax{r A )a\ + h.c. (101) 



m,n A 



This Hamiltonian governs the dynamics of the atom-field system via Heisenberg's equation 
of motion 



ih 



0,H 



(102) 
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Applying the commutation rules d36*b and j 100b . the equations of motions for the photonic am- 
plitude operators and the atomic flip operators read 

-Amn — ^^rnnArfayi -\- — ^ ^ ^ ^ (^nkd-nkA-mk km^-km A-kn) 

k X 

x (A x (r A )a x - A^(r A )4) , (103) 
ax = -iuJxax + A mn d mn ■ A* x (r A ). (104) 

m.n 

In most cases of interest it is sufficient to concentrate on two out of the potentially many atomic 
levels, a ground state \g) and an excited state |e) separated by a transition frequency u e — uo g = 
lo a . The three relevant atomic flip operators, the Pauli operators, will be denoted by a = |<7)(e|, 
Sr = \e) (g\, a z = |e) (e| — \g) (g\. Together with the identity operator I = |e) (e| + \g) (g\ in that 
two-dimensional Hilbert space, they generate the group SU(2). Finally, Heisenberg's equations 
of motion reduce in the rotating-wave approximation to 

& = -iu A a-~d-M+\r A )a z , (105) 

& z = ^•E< + '(r A )ff t +h.c, (106) 
n 

ax = -i^xax + ^d* • A* x {r A )& , (107) 
n 

where the positive-frequency part E^ + \r A ) of the electric field is given by the first term in 
Eq. ( |381 l. We can attempt to solve Eq. ( 1107t by formally integrating it, 



ax(t) = e-^ax + ' A aM / dt ' e^^-^otf) 



(108) 



The first term in this equation is the free evolution of the photonic amplitude operators whereas 
the second term is due to the interaction with the two-level atom. The time integral contains 
the solution to Eq. (1105b which itself is unknown. Thus, the integral can only be computed 
approximately. For this purpose, we split up the fast time evolution e _MAt from the atomic flip 
operator and define a slowly-varying quantity a(t) as 

d(t) = a{t)e lu)At . (109) 

Then, the time integral can be approximated as 



dt'e-^t-^ait') = J dt' e -*"*(*-*') e -*"At' 
o 

t t 
a(t) { dt' e -i*»(t-t') e -iu A t' =s .^ f dt i e <(uu-u»x)(t-t') (110) 
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Fig. 6. Functions s(x) and c(x) showing how they can be approximated by the functions ttS(x) and 
V(l/x), respectively. 

where the atomic flip operator has been taken out of the integral at the upper time. This is only 
possible if the amplitude operator d(t) is almost constant over the time scale \oja — ^a| _1 - This 
in fact also means that the interaction between the electromagnetic field and the two-level atom 
must not be too large. Because the atomic flip operator has been taken out of the time integral, all 
memory effects of the atom-field interaction have been neglected. This is known as the Markov 
approximation. 

The remaining integral can be easily evaluated as 



d£ 



(111) 

which we have split into its real and imaginary parts. The function s{u>a — u\) is sharply peaked 
at io\ = to a (Fig.©. If all quantities that contain either of these functions are averaged or cannot 
be resolved over sufficiently long times, then we can set 

. 1 



s(x) i — ► tt5(x) , c(x) 



(112) 



[V: principal value] thereby introducing little error. 

If we re-insert the formal solution to Eq. ( 1107b into Eqs. ( 11051 ) and ( 11061 l. we obtain the 
equations of motion of the atomic flip operators in the Markov approximation as 



= -r(i 



fatEWM-d-faEi-iM-d*, 



(113) 
(114) 
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where E {lcc (rA,t) denotes the freely evolving electric field operator which, for example, de- 
scribes the action of an external (classical) driving field. The symbols V and Su) are abbreviations 
for the following objects: 

r = ^^^|A A (r A )-d|^(^- WA ), (115) 

A 

Soj = ^Y. V {—^- ) |AAM-d| 2 , (116) 
whose relevance will become clear in the next section. 



2.2.2 Spontaneous decay and Lamb shift 



The equation of motion for the population inversion operator & z , Eq. (II 14b . can be solved easily 
if no external electric field is present. In this case the equation of motion reduces to 

& z = -T(l + a z ). (117) 

Rewriting the inversion operator in terms of the projectors onto the excited and ground states, <j z 
= |e)(e| — \g)(g\ = cr e e — &gg, we find for the excited-state projector the simple relation 

a ee = -ra ee ^ & ee (t) = e- r( - t - t ' ) a ee {t r ) . (lis) 

Hence, the quantity V determines the rate with which a two-level atom decays spontaneously 
from its excited state to its ground state. 

Using the expansion d38l ), the decay rate Eq. dl 15t can be written as 

r = % E d ' ® &-\v A )\0) ■ d*5(u A - cox) , (119) 

A 

which has to be understood in such a way that the 5 function is placed under the mode sum. 
Hence, the rate with which the atom loses its excitation depends on the strength of the vacuum 
fluctuations of the electric field at the frequency of the atomic transition. In a certain sense, spon- 
taneous decay can be viewed as stimulated emission driven by the fluctuating electromagnetic 
field. Using the plane-wave expansion d34b . we obtain the well-known result for the spontaneous 
decay rate in vacuum 

oo 

2ir ^ f u} 2 dui f 7 ^fkv . . w^ldl 2 

r « - ¥ E J J M^-K- ^ A - U ) = ^LL . (120) 

CT 

From the above, it should have become clear that the rate of spontaneous decay can be mod- 
ified by altering the mode structure of the electromagnetic field. We have seen previously in 
Sec. 12.1.41 that the presence of boundary conditions for the electromagnetic field modifies the 
mode structure and hence the strength of the vacuum field fluctuations at the atomic transition 
frequency loa- As a simple example, let us consider a radiating dipole located close to a per- 
fectly conducting mirror. Depending on its orientation with respect to the mirror surface, its rate 
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(b) 



d' 



d 



Fig. 7. An atomic dipole d in front of a perfectly conducting mirror and its image d' for (a) parallel and 
(b) perpendicular orientation of the dipole with respect to the mirror. 



of spontaneous decay is either completely suppressed or doubled with respect to its free space 
rate (Fig. |7J. Suppression occurs when the dipole is parallel to the mirror and hence the dipole 
and its image are antiparallel and cancel each other [Fig. 13 a)]; doubling follows for perpendicu- 
lar orientation where the dipole and its image are parallel [Fig. 0b)]. 

The quantity 8uo that arises in the context of weakly interacting systems in the Markov ap- 
proximation, Eq. dl 16b , induces a shift of the atomic transition frequency, the Lamb shift. Rewrit- 
ing the Lamb shift using the plane-wave expansion ( f34t . we find 



Su> = 



1 



h 2 c 3 



J2 V 



dui 



(27r) 3 uja — to 



dil e 

2eo' 



dto- 



(121) 



which is clearly infinite. This is yet another artefact of quantum theory in free space which can 
be remedied by mass renormalisation (for details, see e.g. Ref. [23]). For our purposes it is 
sufficient to argue that the 'bare' atomic transition frequency uja is unobservable because the 
interaction with the electromagnetic vacuum field can never be switched off and hence the only 
observable quantity is the renormalised frequency uja + 8u>, However, in the following section 
we will again show how the presence of boundary conditions and, in particular, dielectric matter 
can modify the Lamb shift by an additional finite (and thus measurable) amount. 

At this point it is perhaps interesting to observe that spontaneous decay and the Lamb shift 
are intimately connected by a causality relation analogous to the Kramers-Kronig relations that 
we will encounter in the next section. In fact, it is easy to see that T and Suj, taken as functions of 
the atomic transition frequency uja, form a Hilbert transform pair. Rewriting the integral (111 lb 
as 



lim / dt* e *(«A-u»x)(t-*') 

i— ooj 




dre i(u 



dr e l(wA ^ A)T 6(r) , 



(122) 



one observes that this is nothing else than the Fourier transform of the Heaviside step function 
<d(t). This in turn can be interpreted at the causal transform of the function f(ui) = 1, and 
the functions s(x) and c(x) are its real and imaginary parts. Due to the definition of a causal 
transform, and by Titchmarsh's theorem [24], s(x) and c(x) and therefore T and 8uj are Hilbert 
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transform pairs and mutually connected via Rramers-Kronig relations. Hence, knowledge of the 
spontaneous decay rate at all frequencies gives access to the Lamb shift and vice versa. 



2.2.3 Optical Bloch equations 



In this section we will return to Heisenberg's equations of motion for the atomic quantities in 
Markov approximation, Eqs. (II 13b and (11 14b . and solve them under the assumption of an external 
driving field prepared in a single-mode coherent state | a) with frequency uj\ . Hence, we set 



a\T&£jAr A ,t)\a) = iuj x A x (r A )ae 



-iU!\t 



a \^Uei r A,t)\a) = -iu x A* x (r A )a*e 



(123) 



In a frame that co-rotates with the angular frequency u>\, Heisenberg's equations of motion re- 
duce to 



—cr H <7 Z , 

2 2 



iSa 

-r(i + & z ) - - 



(124) 
(125) 



where <5 = uj\ — ui A — Sui is the detuning and f2 = 2 ^ a A\(r A ) ■ d denotes the Rabi frequency. 

The density operator of any two-level (or spin- 1/2) system can be written in terms of the 
Pauli operators as 



I + u- & 



(126) 



where u = (it, v, w) T is a real vector with norm |u| < 1 and & is the vector of Pauli matrices. 
Converting Eqs. (1124b and ( 1125b into equations of motion for the Bloch vector u one finds the 
matrix equation 



(127) 



For negligible spontaneous decay T, the Bloch equations take on the form of the gyroscopic 
equations 
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u = (3 x u 

with j3 = (— Oj, £Ir, 5) T . Its solution is then given by 



u(t) 
v(t) 
w(t) 



n 



sn 
p + n 2 
n 2 



im{Vs 2 + m) , 

1 - cos( \/S 2 + QH) 



■1 



5 2 + tt 2 



1 - cos( y/8 2 + Q 2 t) 



(128) 

(129) 
(130) 
(131) 
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Hence, the time evolution of the Bloch vector on time scales t <C 1/r is described by precession 
of the Bloch vector with frequency \J ft 2 + il 2 along the surface of a cone with opening angle 

±arccos[«V(<5 2 + ^ 2 )]. 

In the long-time limit, t — > oo, the solution of the optical Bloch equations reaches its steady- 
state value. In the stationary state, when u = 0, the solution to the full Bloch equations is 

SQ _ jffi _ 

^- + S 2 + (^) 2 \ + 5 2 + {±) 2 \ + 5 2 + {±) 2 

In the weak driving limit, 51 <C (r, S), the induced atomic dipole moment u st + iv s t takes the 
form 

Q 

u st + W st ~ —z rprrr (133) 

o — il / 2 

whose real and imaginary parts fulfil the Kramers-Kronig relations when integrated over the 
mode frequency uj\. Recall that both the spontaneous decay rate T and the Lamb shift Su> (via 
the detuning S) are contained in (11331) . In this weak-coupling regime, the induced dipole is thus 
described by a linear susceptibility. Its imaginary part, being proportional to the spontaneous 
decay rate T, describes a loss channel for the incident field. The concept of linear-response 
functions and their role in the quantisation of the electromagnetic field in the presence of mag- 
netoelectric matter will be detailed in Sec. [3] 



2.2.4 Jaynes-Cummings model 

In a sense the opposite limit to the case described in Sec. 12.2.31 is obtained if one considers a 
situation in which the mode structure of the electromagnetic field has been altered in such a way 
that discrete field modes interact with the atomic system. We have previously seen in connection 
with the Casimir effect (Sec. I2.1.41 i that this can be achieved in resonators of Fabry-Perot type 
where the allowed modes in a cavity of length d have discrete wave numbers k z = mr/d. The 
half distance between two neighboring modes, c/ (2d), is called the free spectral range. 

If a two-level atom with transition frequency loa is almost resonant with one of the cav- 
ity modes, one can treat the coupled atom-field system approximately by a single-mode model 
described by the Jaynes-Cummings Hamiltonian [25-27] 

H = huia^a + -Itljao'z — fig {aa^ + o^cr) . (134) 

The coupling constant g, which we have given the dimension of a frequency, can be read off 
from Heisenberg's equations of motion as g = -^pd • A\(ta)- Because we have assumed near 
resonance between atomic transition and the relevant cavity field mode, we have employed the 
rotating-wave approximation and subsequently dropped counter-rotating terms in the Hamilto- 
nian ( 11341 i. 

The assumption that effectively only a single field mode interacts with the two-level atom 
implies a sharply peaked, comb-like, density of field modes. This requires a discretisation of the 
modes inside the cavity that can only be achieved with (almost) perfectly reflecting cavity walls. 
In reality, the material making up the cavity mirrors shows some transmission, part of which is 



Macroscopic QED — concepts and applications 



29 




Fig. 8. Equivalent potential model of a cavity with highly reflecting walls (left panel). The transmission 
coefficient is sharply peaked at the cavity resonances (right panel). The line widths of these resonances 
decrease with increasing reflectivity of the cavity walls. 



of course wanted in order to be able to probe the cavity field from the outside. In effect, each of 
the cavity mirrors can be treated as a beam splitter (Sec. 12.1.3b . For one-dimensional light prop- 
agation the equivalent potential encountered by the vector potential is a double barrier (Fig. [8]). 
Analytical solutions for the transmission coefficient have been first obtained in Ref. [28]. How- 
ever, the problem becomes simpler by assuming that the barriers are S functions with strength 
g [29]. Then the transmission coefficient near one of the cavity resonances loc can be written as 
r 

7» ~ ^— (135) 

where the line width V is inversely proportional to the barrier height g. Associating the bar- 
rier height with the squared index of refraction of the mirror material, one finds that the better 
the reflective properties of the mirrors [recall that r = (^fe — + 1)] the narrower the 

resonances. 

The Jaynes-Cummings model is one of the few exactly solvable models of interacting quan- 
tum systems. The Hamiltonian ( 1134t is in fact block diagonal in the basis {|n, e), \n + 1, g)}, 

fr ( \n,e) nu + \u A -gVn + 1 \ f \n,e) \ 

U \\n + l,g) ) n \ -gV^+T (n + l)w-^ A J \ \n + l,g) J ' 
Its eigenfrequencies are 

w n ,± = ( n+ I ) lo± ^A n (137) 



2 2 



where we have defined the Rabi splitting A n = ^/ S 2 + fi^ which depends on the detuning 6 = 
u>a — <jJ and the n -photon Rabi frequency 0„ = 2gy / n + 1. The eigenstates \n, ±) of the Jaynes- 
Cummings Hamiltonian are superpositions of the unperturbed eigenstates {\n, e), \n + 1, <?)}, 



\n, +) \_f cos6„ -sin9„ \ / \n + l,g) 
jn, — ) / \ sin9 n cos0„ J \ \n,e) 

where the rotation angles are 



(138) 



sin9„ = — , cos0„ = — ; (139) 

y/(A n - 5Y + ni y/(A n 5Y + ni 
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Fig. 9. Energy level diagram of the resonant Jaynes-Cummings model. Apart from the collective ground 
state |0, g), the states {|n, e), \n+ 1, g)} are doubly degenerate. This degeneracy is lifted by the interaction. 

For resonant interaction, 6 = 0, the unperturbed eigenstates are pairwise degenerate (Fig. |5J. 
This degeneracy is lifted by the atom-field interaction. The level splitting A„ depends on the 
number of photons. Note that even if there is initially no photon present, the exact eigenstates 
will be split by an amount fin = 2g, the vacuum Rabi splitting. This splitting, or 'dressing', of 
the bare energy levels is equivalent to the Lamb shift which we found in Sec. 12.2.21 However, 
it should be noted that the vacuum Lamb shift was due to the interaction with electromagnetic 
modes of all frequencies, whereas in the Jaynes-Cummings model the level shift arises from the 
interaction with a single discrete mode that has been selected by the resonator. 

Because the Jaynes-Cummings Hamiltonian can be explicitly diagonalised, the unitary evo- 
lution operator is also known explicitly and reads 



U(t) = e 



_ -iHth _ mi a 



= e 



^ 2 |0, ff )<0, 9 | + ^£, 



\n,a)(n,a\ 



a=± n=0 



oo 

= e iUAt / 2 \0,g){Q,g\ + £ e -H»+i/*)«t J e -<A„*/2| nj + \^ + | + e ^/2^ _)( nj _| 



e^ At/2 \0, g)(0,g\ + £ e -i{n+l/2)uit f 



cos^ +i^-sin%i 



\n + l,g)(n + l,g\ 



cos ~ i'KZ s ^ n ~i 



f^] \n, e)(n, e| + ifc sin %^ (\n + l,g) (n, e| + |n, e) (n + 1, ,g|)} ■ 

(140) 



The last equality in Eq. ( |140| > expresses the unitary time evolution in terms of the unperturbed 
eigenstates. 

An important special case is the dispersive limit in which the detuning 6 is large compared 
to the relevant ?i-photon Rabi frequency, 6 3> f2„. In this limit, the level splitting can be approx- 
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imated by A„ ~ S + 0^/(2$) = 6 + 2g 2 (n + 1)6, and the exact eigenstates are approximately 
the unperturbed eigenstates, \n, +) ~ \n, e) and \n, — ) ~ \n + 1, g). Under this approximation, 
the unitary evolution operator can be written as 



U ~ e^ 2 \0,g)(0,g\ + J2, 



-i(n+l/2)ujt 



-^/ 2 +s 2 (»+ 1 )/^|n, e) (n, e| + e^/ 2 + 92 (" +1 '/^|n + 1, 5 ) (n + 1, ff | 



or, with the free Hamiltonian Hq = huia'a + \TujJa^z, 



(141) 



1 „ 

U ~ exp I ——H t 

h 



exp ( - l -^-(n + 1)^ |e)(e| + exp f^-nt } \n) 



(142) 



It is instructive to note that the evolution operator the dispersive limit, Eq. ( 1142b . is quadratic in 
the interaction strength g. This means that it results in an effective nonlinear atom-field interac- 
tion. In order to investigate this claim in more detail, let us rewrite Eq. (11421 ) in the following 
form. The term in square brackets contains, apart from a linear Stark shift of the excited state 
|e), a factor e~' 9 na * t / s which is the result of an effective nonlinear Hamiltonian 



H. 



off 



n—a' aa z 
o 



(143) 



Comparing Eq. d 143b with the Jaynes-Cummings Hamiltonian ( 1134b we see that the trilinear 
Hamiltonian ( 1143b does not appear in the original Jaynes-Cummings model. The appearance 
of such an effective Hamiltonian is solely due the far off -resonant interaction. This is just one 
example of a generic nonlinear interaction such as those studied in Sec. 13.41 
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3 Macroscopic quantum electrodynamics 



Having established the framework of quantum electrodynamics in free space, we are now in a 
position to generalise the theory to magnetoelectric background materials. Before we go ahead 
with our program, let us first discuss the intrinsic difficulties associated with magnetoelectric 
media. 

Let us assume we wanted to naively extend a plane-wave expansion of the electromagnetic 
field to include dielectrics. We would then try to replace the plane wave solutions e ik r to the 
Helmholtz equation by e in where n = n(ui) is the index of refraction of the dielectric. In 
order to conform with standard requirements from statistical physics, the refractive index must 
be a complex function of frequency n{ui) = r)(w) + tK(ui) that satisfies the Kramers-Rronig 
relations 



Due to the inevitable imaginary part of the refractive index, the plane waves are generically 
damped. That in turn means that they do not form a complete set of orthonormal functions 
needed to perform a Fourier decomposition of the electromagnetic field. The consequences of 
this failure are quite severe; either one insists on bosonic commutation rules for the photonic 
amplitude operators ci\ and d\ which subsequently lead to wrong commutation relations between 
the operators of the electric field E(r) and the magnetic induction B(r), or one postulates the 
correctness of the latter and ends up with amplitude operators in their Fourier decomposition that 
do not have the interpretation of annihilation and creation operators of photonic modes. 

The reason for the failure of this naive quantisation scheme is easily found. The introduction 
of the index of refraction n(uj) means that there exists an underlying (microscopic) theory that 
couples the free electromagnetic field to some dielectric matter, the effect of which is taken into 
account only by means of the response function n(ui). In doing so, the matter-field coupling is 
hidden from view but is nevertheless present. The damped plane waves e mk r have therefore to 
be regarded as eigensolutions of the combined field-matter system, and not of the electromagnetic 
field alone. 



The above arguments necessarily lead one to consider the electromagnetic field interacting with 
an atomic system coupled to a reservoir that is responsible for absorption. An explicit matter-field 
coupling theory that achieves field quantisation in dielectric matter on the basis of a microscopic 
model has been developed by Huttner and Barnett [30, 3 1 ] (Sec. 13.1.11 ). This Hamiltonian model 
can be generalised to an effective Langevin noise model (Sec. 13.1.21 ) which forms the basis of the 
remainder of this article. 



oc 



oc 




(144) 



3.1 Field quantisation in linear absorbing magnetoelectrics 



3.1.1 Huttner-Barnett model 



Historically, the first successful attempt at quantising the electromagnetic field in an absorbing 
dielectric material is due to Huttner and Barnett [30,31]. They considered a Hopfield model [32] 
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of a homogeneous and isotropic bulk dielectric in which a harmonic oscillator field representing 
the medium polarisation is linearly coupled to a continuum of harmonic oscillators standing for 
the reservoir (first line in Fig.lTot. Such a model leads to an essentially unidirectional energy flow 
— from the medium polarisation to the reservoir — which means that the energy is absorbed. 
Strictly speaking, because a single harmonic oscillator is coupled to a continuum, the revival time 
is infinite, hence an excitation stored in the continuum of harmonic oscillators will not return to 
the medium polarisation in any finite time. The overall system of radiation, matter polarisation, 
reservoir and their mutual couplings are regarded as a Hamiltonian system whose Lagrangian 
reads 



L= d 3 r£= d 3 r (£ em + £ mat + £ int ) 



where 



£o 
2 



c 2 (V x A) 



C 



mat 



'Y^,-c 2 Y 2 



)• 



DC 

-a (A • X + cj)V • x) - J duo v(ui)X ■ Y u 



(145) 

(146) 
(147) 

(148) 



Here £ em and £ ma t are the free Lagrangian densities of the radiation field and the matter, respec- 
tively, where <fi and A are the scalar and vector potentials in the Coulomb gauge (V • A = 0), 
and X and Y w the medium and reservoir oscillator fields with density \i, respectively. In the in- 
teraction part, A n t, a is the electric polarisability and the medium-reservoir coupling constants 
v(w) are assumed to be square integrable. 
Upon introducing the canonical momenta 



„ dC 
II = = e A , 

dA 

dC ■ 
f = — - = ua — aA , 

ax 

CL = -=7— = - v(uj)x , 



(149) 
(150) 
(151) 



one can perform the Legendre transformation and construct a Hamiltonian H = H em + H mat + 
Hi nt . In Fourier space, 

2 

A(r) ^ A(k) = ^^ A (k)e A (k), (152) 

A=l 

2 

X(r) ^ X(k) =X||(k)e k + ^X A (k)e A (k), (153) 



A=l 
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are the longitudinal and transverse components of the vector potential and the matter polarisation, 
respectively. Similar decompositions are made for all other fields. 
As in free space, one introduces mode amplitudes according to 



a A (k) 

&*(k) 
b X (k,u) 



2hkc 

JJ_ 

2fuj 

JL 

2hu 



kcA x {k) + — n A (k) 
uX x (k) + -Q x (k) 



-icvX\(k,Lu) + -Q x (k,uj) 



(154) 
(155) 
(156) 



where 



u 2 



du> 



v 2 (lo) 



(157) 



The expressions ( 1157b reflect the level shifts due to the interaction between fields. Indeed, we 
have encountered such shifts already in vacuum QED (Lamb shift, dressed energy levels in the 
Jaynes-Cummings model etc.). Similar decompositions can be made for the longitudinal fields 
which we will not consider here [33]. 

The amplitude operators are then promoted to Hilbert space operators with the usual bosonic 
commutation relations 



a A (k),4,(k') = s xx ,6(k-k'), 

St 



Mk),& A ,(k'; 



b\(k,Lu),b\>(k' ,lo' 



8 xx ,S(k-k'), 

5 X y6(k-k')5(LU~L0'). 



(158) 
(159) 
(160) 



The transverse Hamiltonian can be expressed in terms of the annihilation and creation operators 

as 



with 



H = H c , 



2 

E 

A=l 
2 

E 



d 6 khkc a\(k)a x (k) 



d 3 k\tujb[(k)b x {k)+ / dujhojb[(k,uj)b x (k,u}) 



(161) 



(162) 



+- I du V(u) b\(k)+b x (-k) b\(-k,uj) + b x (k,Lu) 



ih x - 



d 3 kA(k) a A (-k)+a A (k) b{(k)+b x (-k) 



(163) 



(164) 
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Fig. 10. Two-step diagonalisation of the Huttner-Barnett model. In the first step, the polarisation field and 
the harmonic-oscillator heat bath form dressed-matter operators. These are then combined in the second 
step with the free electromagnetic field to form polariton operators. 

where V(u>) = [v(us)/ ix\{lo /lj) 1 ^ 2 and A(fc) = [d>a 2 / (/iceofc)] 1 ^ 2 ■ The Hamiltonian is clearly 
bilinear in all annihilation and creation operators, and can therefore be diagonalised by a Bogoli- 
ubov (squeezing) transformation, that is, by a linear transformation involving both annihilation 
and creation operators. In the present context, the procedure is known as a Fano-type diagonalisa- 
tion [34]. The diagonalisation is performed in two steps. In the first step, the matter Hamiltonian 
H m &t is diagonalised first (second line in Fig.fTOli, leading to 

2 00 

£mat = J2 J d3k J dujhLoB{(k,Lu)B x (k,u;). (165) 







In the second step, the dressed-matter operators B\(k, uj) and E>\(k 7 uj) are combined with the 
photonic operators to the diagonal Hamiltonian 

2 00 

H = J2 fd 3 k f diuhu;fl(k,u)h(k,Lo). (166) 
x=i J { 

Diagonalisation of the longitudinal field components can be achieved analogously. Adding 
the resulting expression to Eq. ( 1166t and Fourier transforming gives 



H = J frfcf(r,u) -f(r,w) (167) 

which is depicted in the last line in Fig. [10] Due to the bosonic commutation relation of the 
amplitude operators, the commutation rule for the new dynamical variables are 

f{r,cj),P(r',u')] = S(r - r')5(uj - uf) . (168) 
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Inverting the Bogoliubov transformation that has led to the polariton-like operators C\ (k, u>) 
and C\ (k, ui) and subsequent Fourier transformation leaves one with an expression for the vector 
potential A(r) and the polarisation field P(r) in terms of the dynamical variables f(r, w) and 
F(r, u). The expansion coefficients turn out to be the dyadic Green tensor for a homogeneous 
and isotropic bulk material with a dielectric permittivity that is constructed from the microscopic 
coupling parameters a, v(uo) and m [30, 35]. Later, this theory has been extended to inhomoge- 
neous dielectrics where Laplace transformation techniques have been used to solve the resulting 
coupled differential equations [36]. However, neither of these expressions for the resulting fields 
contains any hints towards their underlying microscopic model, so it seems quite natural to start 
from the source-quantity representation of the electromagnetic field instead. 



3.1.2 Langevin-noise approach 

From now on, we leave the microscopic models behind and concentrate on the phenomenological 
Maxwell's equations, assuming that the relevant response functions such as dielectric permittivity 
and magnetic permeability are known from measurements. Maxwell's equations of classical 
electromagnetism, in the presence of magnetoelectric background media read, in the absence of 
other external sources or currents, 

V-B(r) = 0, (169) 

V x E(r) = -B(r) , (170) 
V-D(r) = 0, (171) 

V x H(r) = D(r). (172) 

They have to be supplemented by constitutive relations that connect the electric and magnetic 
field components. Assuming for a moment that the medium under consideration is not bian- 
isotropic, we can write 

D(r) = £o E(r) + P(r) , H(r) = — B(r) - M(r) (173) 

Mo 

where P(r) and M(r) denote the polarisation and magnetisation fields, respectively. 

Polarisation and magnetisation are themselves complicated functions of the electric field 
E(r) and the magnetic induction B(r). Assuming that the medium responds linearly and lo- 
cally to externally applied fields, the most general relations between the fields that are consistent 
with causality and the linear fluctuation-dissipation theorem can be cast into the form 

oo 

P(r, t) = e J dr Xe (r, r)E(r, t - r) + P N (r, t) , (174) 
o 

oo 

M(r, t) = — fdr x m (r, r)B(r, t - r) - M N (r, t) (175) 
Mo J 
o 

where Pu(r, t) and Mo(r, t) are the noise polarisation and magnetisation, respectively, that are 
associated with absorption in the medium with electric and magnetic susceptibilities x c (r, r) and 
Xm(r,r). 
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The Fourier transformed expressions ( 1174b and ( 1175b convert the constitutive relations ( 1173b 
into 

D(r,w) =e e(r,w)E(r,cj)+P N (r,cj), H(r, w) = Ko/s(r, «)B(r, w)-M N (r, w) , (176) 
[k(t,uj) = /i _1 (r,o;)] where 

oo oo 

e(r, u) = 1 + J dr Xo (r, r)e iuT , «(r, w) = 1 - / dY Xm (r, r)e MT (177) 
o o 

are the relative dielectric permittivity and (inverse) magnetic permability, respectively. An im- 
mediate consequence of the causal relation (1177b is the validity of Kramers-Kronig (Hilbert 
transform) relations between the real and imaginary parts of the susceptibilities in Fourier space, 

oo oo 

Re x(r, w) = -T 5 / aw , Imx(r,u) = V \ aw ; . (178) 

IT J LO' — UJ 7T J L)' — UJ 

— oo — oo 

Using the expressions ( 1177b for the dielectric permittivity and the magnetic permeability, 
Maxwell's equations for the Fourier components can thus be written as 

V-B(r,w) = 0, (179) 

VxE(r,w) = B(r,w), (180) 

e V- [e(t,u)E(t,u)] = p N (r,w), (181) 

Vx[ K (r,w)B(r,w)]+4e(r,w)E(r,w) = /i j N (r,w). (182) 

Here we have introduced the noise charge density 

ftf(r,w) = -V-P N (r,w) (183) 
and noise current density 

j N (r, u) = -iwP N (r, u) + V x M N (r, to) , (184) 

respectively, which by construction obey the continuity equation. 

Equations ( |181b and (1182b now contain source terms. Hence, the electromagnetic field in 
absorbing media is driven by Langevin noise forces that are due to the presence of absorption 
itself. Moreover, the particular combination in which the noise polarisation and magnetisation 
enter Maxwell's equations, Eq. (1184b . suggests that dielectric and magnetic properties cannot 
be uniquely distinguished. For example, one might include the magnetisation in the transverse 
polarisation in which case the constitutive relations ( 1 173b would have to be altered. This simple 
observation implies that the constitutive relations in the present form cannot be the fundamental 
relations. Instead, the noise current density appears as the fundamental source of the electromag- 
netic field. This becomes even more apparent if one allows for spatial dispersion that makes the 
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dielectric response functions nonlocal in configuration space. It is therefore expedient to rewrite 
Eqs. (USB and <fl82t as 

uj 2 

V x V x E(r, to) - -5-E(r, to) = i/iowj(r, uj) , (185) 

and to consider the most general linear response relation between the current density and the 
electric field in the form of a generalised Ohm's law as 



j(r,u;) = J dVQ(r,r»-E(r»+j N (r,w) 



(186) 



where Q(r, r', to) is the complex conductivity tensor in the frequency domain [37]. 

The Onsager-Lorentz reciprocity theorem demands the conductivity tensor to be reciprocal, 
Q(r, r', to) — Q T (r' , r, uj). The two spatial arguments must be kept separate in general, except 
for translationally invariant bulk media in which the conductivity only depends on the difference 
r— r',i.e. in this caseQ(r, r',oj) = Q(r—r',oj). We assume that, for chosen to, the conductivity 
tensor is sufficiently well-behaved. By that we mean that it tends to zero sufficiently rapidly as 
|r — r'| — > oo and has no non-integrable singularities. However, S functions and their derivatives 
must be permitted to allow for the spatially nondispersive limit. The real part of Q(r, r', uj), 

ct^v^uj) =ReQ(r,r',uj) = ^[Q(T,r',uj) + Q + (r , ,T,to)] , (187) 

is connected with the dissipation of electromagnetic energy and for absorbing media, as an inte- 
gral kernel, associated with a positive definite operator [37]. Under suitable assumptions on the 
causality conditions satisfied by its temporal Fourier transform Q(r, r', t), the conductivity ten- 
sor is analytic in the upper complex uj half-plane, satisfies Kramers-Rronig (Hilbert transform) 
relations, and obeys the Schwarz reflection principle 

Q(t,v',-uj*) = Q*(r,r». (188) 

We now identify the current density ( 1186b as the one entering macroscopic Maxwell's equa- 
tions in the frequency domain. The medium-assisted electric field thus satisfies an integro- 
differential equation of the form 

uj 2 f 

V x V x E(r, uj) - ^E(r, uj) - i[i uj / dV Q(r, r',to)- E(r', to) = ifi uj^(r, uj) . (189) 

The unique solution to the Helmholtz equation (1189b is 

E(r, uj) = i^uj J dV G(r, r', to) ■ j N (r', uj) (190) 

where G(r, r', to) is the classical Green tensor that satisfies Eq. J 1 89b with a tensorial 6 function 
source, 

uj 2 f 
VxVxG(r,s,uj)-—G(r,s,uj)-ifi Q uj / dV Q(r, r', uj) ■ G(r', s, uj) = J(r-s) (191) 
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together with the boundary conditions at infinity. It inherits all properties such as analyticity 
in the upper complex to half-plane, the validity of the Schwarz reflection principle, as well as 
Onsager-Lorentz reciprocity from the conductivity tensor, viz. 

G(r',r,w) = G T (r,r», (192) 
G(r,r',-w*) = G*(r,r',w). (193) 

In addition, the Green tensor satisfies an important integral relation that can be derived as follows. 
The integro-differential equation ( 11911 ) can be rewritten as 



J d 3 s H(r, s, w) ■ G(s, r', to) = S( 



(194) 



where the integral kernel H(r, r', to) = V x V x S(r — r') — to 2 /c 2 8(r — r')—i/j,QOjQ(r, r', lo) 
is reciprocal, from which Eq. ( 1 1921 ) follows. With that, the complex conjugate of Eq. ( 11941 ) reads 



J d 3 s G + (r, s, u) ■ H + (s, r',u) = 5{r - r') 



(195) 



If we now multiply Eq. (1194b from the left with G + (s', r, u>) and integrate over r, then multiply 
Eq. (11951 ) from the right with G(r' , s', oj) and integrate over r', and finally subtract the resulting 
two equations, we find that for real lo the integral equation 

fi u J d 3 s J d 3 s' G(r,s,uj) -(t(s,s',w) ■ G + (s',r',w) = ImG(r,r» (196) 

holds [38] (see also App.lAl. 

Up until this point, all our investigations regarded classical electrodynamics. In order to 
quantise the theory, we have to regard the Langevin noise sources j^(r, uj) as operators with the 
commutation relation 

jVr^J^rV)] =— %-wV(r,r>), (197) 

J 7T 



which follows from the fluctuation-dissipation theorem associated with the linear response ( 1186b . 
In this way, the operator of the electric field strength is given by the operator-valued version of 
Eq. (fT90l as 



OO 

E(r) = J dwE(r,u) + h.c. , E(r,w)=i/i w J d 3 r' G(r, r', u) ■ j N (r', u>) 

o 



(198) 



The consistency of this quantisation procedure can be proven by checking the fundamental 
equal-time commutation relation between the operators of the electric field and the magnetic 
induction. Using Faraday's law, Eq. ( 11801 ), we find the frequency components of the magnetic 
induction field as 

OO 

B(r) = J dt 1 )B(r,w)+h.c, B(r,w) = ^ Vx J d 3 r' G(r, r', lo) ■ j N (r', w) . (199) 

o 
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Hence, the equal-time commutator reads, on using the commutation relation (1197b and the inte- 
gral formula ( 11961 >. as 



E(r),B(r') 



ttsqc z 



:V,/ X 



du id G(r,r' ' ,u>) . 



(200) 



where the second equality follows from the Schwarz reflection principle. Using the analyticity 
properties of the Green tensor in the upper complex oj half-plane, we then convert the integral 
along the real oj axis into a large semi-circle in the upper half-plane. From the integro-differential 
equation (1191b and the properties of the conductivity tensor Q(r, r', oj) we find the asymptotic 
form of the Green tensor for large frequencies as 



G(r,r» H ~ °-£^(r-r') ) 

0J Z 

so that the equal-time field commutator takes its final form of 
E(r),B(r') 



-Vx5(r-r') = --Vx<5 x (r 



(201) 



(202) 



The striking feature is that the field commutator ( 12021 i is exactly the same as in free-space 
quantum electrodynamics [Eq. (l40ll. despite the presence of an absorbing dielectric background 
material. This fact reinforces the view that the fields E and B represent the degrees of freedom 
of the electromagnetic field alone and have little to do with any material degrees of freedom. The 
apparent discrepancy with the notion of the expansion ( 1198b as a medium-assisted electric field 
is resolved by interpreting the Green tensor as the integral kernel of a projection operator onto 
the electromagnetic degrees of freedom. Finally, the Langevin noise currents j^(r,u>) can be 
renormalised to a bosonic vector field by taking the 'square-root' of the tensor er(r, r', oj) (which 
exists because of its positivity in case of absorbing media). Writing 



er(r,r',a;) = J d 3 s K(r, s, oj) ■ K+ (r 1 , s, t 
and defining 



j N (r,w) = ( — 

1 7T 



1/2 



dV K(r, r',Lu) ■ f(r',w) 



where 



f(r,o/),fVV)l = S(u - lj')S(t - r') 



(203) 



(204) 



(205) 



the expansion J 198b of the frequency components of the operator of the electric field strength 
finally becomes 



E(r, to) = i/iou 



— I d 3 r' 



d 3 s G(r, r', oj) ■ K(r', s, u) ■ f (s, u) . 



(206) 
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Hamiltonian: In order to complete the quantisation scheme, we need to introduce a Hamil- 
tonian as a function of the Langevin noise sources jn(i", u>) and jjj(r, uS) or, equivalently, in 
terms of the bosonic dynamical variables f (r, ui) and F(r, ui). Imposing the constraint that the 
Hamiltonian should generate a time evolution according to 

j N (r,u),H =MnO,w), (207) 

the Hamiltonian must be of the form [38] 

oo 

H = irJdto J d\ J dVjk(r,w) • p{v,v' \oj) •j N (r' » (208) 
o 

where p(r, r', ui) is the inverse of the integral operator associated with er(r, r', a;). In terms of 
the bosonic dynamical variables, the Hamiltonian is diagonal, 

oo 

H= I duo I d 3 rhujiHr,uj) -f(r,w) (209) 



which is its most commonly used form [33,39,40]. Perhaps surprisingly, it closely resembles its 
free-space counterpart, Eq. (l37T i. in that it is bilinear in its dynamical variables. The reason behind 
this behaviour is that any linear reponse theory can be derived from an underlying microscopic 
model that is bilinear in its constituent amplitude operators which, after a suitable Bogoliubov- 
type transformation, leads to a Hamiltonian of the form (I2091 >. An example is provided by the 
Huttner-Barnett model of a homogeneous, isotropic dielectric (Sec. l3.1.TT l. 

Spatially local, isotropic, inhomogeneous dielectric media: We now apply the general theory 
to some special cases that are of practical importance. Let us begin with the simplest, and histori- 
cally first, example of a spatially nondispersive, isotropic and inhomogeneous dielectric material 
that shows no magnetic response. The neglect of spatial dispersion makes the conductivity tensor 
Q(r, r',u>) strictly local, so that tr(r, r',ui) = er(r, u))S(r — r'). Furthermore, isotropy means 
that er(r,cj) = cr(r,oj)I. If we then make the identification cr(r, ui) — eoo;Imx(r, ui) where 
x(r, ui) is the dielectric susceptibility, Eq. (1206b becomes 



E(r, u) = hl-^i f d 3 r' y/lm x(V,lj)G(t, r', u) ■ f(r', w) (210) 

V 7T£ C l J 

which yields the well-known quantisation scheme for a locally responding dielectric material 
[33,39,41-46]. 

Spatially dispersive homogeneous bulk media: As a second example, we consider an in- 
finitely extended homogeneous material for which Q(r, r', us) is translationally invariant [38,47]. 
That is, it is only a function of the difference r — r'. In this case, we represent <r(r, r', oj) as the 
spatial Fourier transform 

<x(r,r» = <r(k,u,y k -( r - r '). (211) 
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A similar decomposition can be made for the integral kernel K(r, r', lo). For an isotropic 
medium without optical activity, the Fourier components <x(k, lo) can be written as [37] 

<x(k, lo) = «7|| (Jfc, + ffJ .(A, (i ^) (212) 

and similarly for K(k, lo), where the expansion coefficients have to be replaced by their positive 
square-roots <r^ 2 (k, lo) and a^ 2 (k, lo) 

Clearly, the decomposition ( 1212b is not unique as <r(k, lo) can be equivalently decomposed 
into 

cr(k, w) = (T|| (jfe, w)J-kx 7 (fc, w)I x k (213) 

where 

7 (fc, w) = [*j_(k, lo) - mj (fc, w)] /fc 2 . (214) 

Since both mi (fc, w) and crj_(fc, w) have to be real and positive to yield a positive definite integral 
kernel <r(k, oj), the new variable -f(k, u>) is real, too. However, it does not have a definite sign. 
If, on the other hand, one forces -f(k, u>) to be positive, then the tensor 

K'Qi, lo) = o\ /2 (k, uj)I ± 7 1/2 (fc, uj)I x k (215) 

is the positive square-root of the integral kernel <r(k, ui). The kernels K'{k, lo) and K(k, lo), 
despite being different, are related by a unitary transformation [38]. 



Spatially local magnetoelectric media: The local limit of the above theory has a rather inter- 
esting structure. If one assumes that the functions mi (fc, lo) and ^(k, lo) vary sufficiently slowly 
with k and possess well-defined long-wavelength limits lirufe_>o ctm (k, lo) = <J\\(lo) > and 
rimfc_>o 7(fc, lo) — j(lo) > 0, one finds the approximation 

er(r, r', lo) = a\\ (lo)S(t - r') - j(lo)V x <5(r - r') x V' . (216) 

The full conductivity tensor associated with that real part (real and imaginary parts are related by 
a Hilbert transform) is then 

Q(r, r', lo) = Q (1) (w)5(r - r') - Q^ 2) (lo)V xS(r-r')x V' (217) 
with the identifications 

qW(lo) = -te w[eM-l] , Q {2) {lo) = -*k [1 - k{lo)] /lo , (218) 

where e(lo) is the dielectric permittivity and /i(lo) — 1/k(lo) the (para-)magnetic permeability 
of the medium. Note that the requirement j(lo) > implies that Im k(lo) < for lo > 0, from 
which it follows that fj,(uj — > 0) > 1 [4]. This in turn means that this theory can only describe 
paramagnetic materials. Diamagnetic materials are intrinsically nonlinear as their response func- 
tions themselves depend on the magnetic field and thus are excluded from a linear-response 
formalism. 
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The noise current density that is derived from the kernel 

K'(r, r',uj) = <j\ /2 {uj) S(r - r') =p j 1/2 (uj) V x <5(r - r') (219) 

[which follows from Eq. (I216H can be decomposed into longitudinal and transverse parts ac- 
cording to j N (r , uj) = j N M (r , uj) +j N j_ (r, uj) with 



JN||(r,w) 



Jn_l(,i",w) 



wV Im ^(w)fj|(r,w) , 
J ^^-ujy/lme(uj){j_(r, uj) =F i<J V 



'lm^(w) } 



(220) 



(221) 



The distinction between longitudinal and transverse components of the noise current density 
(and subsequently the bosonic dynamical variables) is essentially a projection formalism, and 
the f|| m ( r : are termed projective variables [38]. 

Another, more frequently used decomposition is obtained by redistributing the electric part of 
the transverse noise current density. In this way, two new sets of independent bosonic variables, 
f e (r , oj) and f m (r , uj), are introduced that lead to an equivalent decomposition of the noise current 
according to [33,40] 



j Ne (r,w) = 



jNm(r,w) 



y^^w v/lm e(r, uj) f e (r, uj) , 



' fiK _ 
/\/— V x 



' Im /i(r, uj) 
Kr,w)| 2 



fm(r, uj) 



(222) 
(223) 



in which the possible spatial dependencies of the dielectric permittivity and the paramagnetic 
permeability have been reinstated (see Ref. [38] for details). The Hamiltonian ( |209t takes the 
form 



d r / doj fku fl(r, uj) ■ f\(r, uj) . 



A=e,m 



Finally, the electric field ( 1206b and the magnetic induction can be written as 

E(r,w) = 

B(r,w) = 

with the definitions 
G e (r,r» 



^ / dVG A (r,r>).f A (r>), 
- V / dV [VxG A (r,r»].f A (r» 



., -Ime(r,w)G(r,r',w), 
tr V tt£o 



r Cr r' f /l - ^ \ % Im ^ r ^) 
Lr m (T,r ,UJ) — j— -, rry 

c V 7re lM(r,w)| 2 



G(r,r',w) x V 



(224) 

(225) 
(226) 

(227) 
(228) 
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where G(r, r', u>) is the usual classical Green tensor satisfying Eq. (11911 ). The latter Helmholtz 
equation condenses to (see also App.lAl 



V x /t(r, lu) V x G(r, r', u) je(r, w)G(r, r', u) = S(r - r') . 



(229) 



For completeness, we mention that the integral relation ( 11961 > can be cast into the form 

J d 3 sG A (r,s,w)-G}(r',s,w) = - ^ImG(r,r». (230) 



A— e,m * 



Alternatively, the electric and magnetic noise current densities ( 1222b and ( |2231 l can be recast 
into the form of noise polarisation and magnetisation fields 



JncO, cj) = -icjP N (r, u) , j Nm (r, w) = V x M N (r, w) . 
Then, the electromagnetic fields read in terms of these noise fields as 



E(r,w) 



dVG(r, r',w)-P N (r',a;) 



-ifiQU d 3 r' G(r,r»x V -M N (r',w). 



B(r,w) = ifiouj / dV[VxG(r,r»]-P N (r» 



— Mo / ^ 



j3„/ 



V x G(r,r',u;) x V ■M N (r», 



D(r,w) = -i — e(r,w) / dV G(r,r',w) x V -M N (r',w) 



/ dV [e(r, w)G(r, r', w) + J(r - r')] • P N (r', w) , 



H(r,w) 



/i(r,w) 



dV 



dV [V x G(r,r',cj)] -P N (r',w) 



V x G(r,r J ,u>) x V 
M(r,w) 



+ 6(r - r') 



(231) 



(232) 



(233) 



(234) 



M N (r». (235) 



Statistical properties: Thermal expectation values of the electromagnetic field can be obtained 
from those of the dynamical variables. Assuming the electromagnetic field in thermal equilib- 
rium with temperature T, it may be described by a (canonical) density operator 



Qt = 



e -H F /(k B T) 
tle -H F /(k B T) 



(236) 
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[ks' Boltzmann constant]. Thermal averages (. . )t = tr[. . . qt] of the dynamical variables are 
thus given by 



(f A (r,w)®f A /(r / J w / ))r 
(f](r,o;)® f A ,(r', u/))t 



(fA(r,w)) T = = (fj[(r >W )) T , 

n th (uj)Sx\>8(r - y')5(uj - a/) , 



= <ft(r, W )®fJ,(rV)>T, 



(f A (r,u;)®ft,(r>')>T = [tithH + l]*AV«(r-r / )*(w-a/), 



where 



nth(^) 



E, 



-mhuj/(k B T) 



V P -mfiw/(fc B T) p hui/{k B T) _ 1 



(237) 
(238) 
(239) 
(240) 



(241) 



is the average thermal photon number. This translates into the statistical properties of the elec- 
tromagnetic fields as follows: 





<E(r 




= = <Et(r,w)) 


T , 






(242) 


(E(r,w) 


g)E(r', 


W')>T 


= 0=(E t (r,w)( 


3Et(r',a/)) T , 






(243) 


(Et(r,o;) 


8>E(r', 


o/))t 


h uj 2 

= n t h(w)-rImG(r,r»(S(w 

ire c A 




), 


(244) 








(E(r,w)€ 




&/)>t 


r / -, 

= — Kh(w) + 

7T£o 


1] -^ImG(r,r', 

cr 


uj)5(u) — lo' 


(245) 




(B(r 


, w))t 


= = (Bt(r,w)) 


T , 






(246) 


(B(r,w)< 


E> B(r', 


a/)>T 


= = (B t (r,w) 


»Bt(r>')>T, 






(247) 


Bt(r jW )( 


3 B(r', 


w'))t 


= n t h (w Jim 


V x G(r,r',w) 


x V 


8(u- 


to'), (248) 


B(r,w) <8 


B^r', 


o/))t 


= — Kh(^) + 


1] Im [v x G(r 


r',w 


x V 


(5(w — u> ) . 



(249) 

These expressions will be needed for the calculation of relaxation rates (Sec. and dispersion 
forces (Sec. [5}. 



3.1.3 Duality transformations 

An important symmetry of the Maxwell's equations in free space is duality where interchanging 
electric and magnetic fields yields the same differential equations. Here we will show that this 
type of symmetry can be established even within the framework of macroscopic quantum elec- 
trodynamics. At first, we consider macroscopic QED without external charges and currents. We 
group the fields into dual pairs and rewrite Maxwell's equations as 



9_f 1\ (y/J^D 
dt \-l QJ 



(250) 
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where the constitutive relations are combined to 

Va^dA = i (V^v\ ( \ 

^B) c + \^ MJ ■ 

A general rotation T>{0) in the space of dual pairs can be written as 

;)*-»<)• zS)- 

It is easily checked that Maxwell's equations ( 1250b in free space with constitutive relations ( 125 lb 
is invariant under rotations of the form ( 12521 ). 

In the presence of spatially local magnetoelectric materials, the constitutive relations in fre- 
quency space can be further specified to 



V7^D\ = l(e 0\ (y/e^E\ (1 Q\( Vm^Pn 
^Bj c \0 n) \^UJ \0 fx) V^moM n 



(253) 



Invariance of the constitutive relations (12531 l under the duality transformation (1252b requires that 

e* 0\_ , m /e 0\ i, . _ /e cos 2 6» + sin 2 6» (/x - e) sin6>cos6»\ 
a*V~ W V« A*J W ~ W-e) sine cos e sin 2 + /i cos 2 9) {Z ™ } 

which can be fulfilled in two ways. The first is obtained if the dielectric permittivity of the 
material equals its magnetic permeability, e = /i. This is achieved in free space as well as by 
certain metamaterials, for example by a perfect lens with e = fi = — 1 [48]. In this case duality 
is a continuous symmetry which holds for all angles 9. 

Generally, duality holds only for discrete values of the rotation angle, 9 = mr/2 with n E Z. 
In this case, the transformation results in 

e\*_ /cos 2 9 sin 2 9\ /e\ / yffiQp^ \ _ ( cos 9 /zsin#\ / ^//TqPn 



\i) ^sin 2 ^ cos 2 9J\^)' ^/i M N / \- £ 1 sin9 cos 9 J ^^(jqMn 

(255) 

It should be remarked that, not only are Maxwell's equations invariant under the discrete duality 
transformation, but also the Hamiltonian that generates them. To see this, one can derive the 
transformation properties of the dynamical variables from the relations (1255b which read for 

9 = nn/2 as 



f e \ _ / cos 6? —i(fi/\fj,\)wo.Q\ ( f e 

W ~ l-*(|e|/e)sinfl cos9 ) {{„ 



(256) 



These transformations obviously leave the Hamiltonian ( 1224b invariant. Combining all relevant 
transformation relations, we can collect the duality relations for all electromagnetic fields, the 
linear response functions, and dipole moments, in Tab. |2] These relations allow one to establish 
novel results for magnetic (electric) materials and atoms in terms of already known results from 
their corresponding dual electric (magnetic) counterparts. Finally, duality does not only hold 
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Partners 



Transformation 



E, H: 
D, B: 
P, M: 

Pa,M a 
d, rh: 

Pn,Mn 



E* = c^ H, 
D* = c£ B, 
P* = M/c, 
P£ = M A /c 

d* = m/c, 
P* = mMn/c, 



H* = -E/(c Mo ) 
B* = -D/(ce„) 



M* = -cP 



= -cP.4 



rh* = cd 




= -«(M/lMl)fm, 



a* = (3/c 2 , 



M* = -cP N /e 

[1* = E 

(3* = c 2 a 



Tab. 2. Effect of the duality transformation. 



on the operator level in macroscopic QED without external charges and currents, but can also 
be established for derived atomic quantities. It is shown in Ref. [49] that dispersion forces as 
well as decay rates are all duality invariant, provided that the bodies are stationary and located 
in free space and that local-field corrections are applied when considering atoms embedded in a 
medium. The duality invariance of dispersion forces is further discussed and exploited in Sec. [5] 



In a previous section (Sec. 12.1.3b we developed the theory of quantum-state transformation at 
lossless beam splitters which accounts for a unitary transformation between the photonic am- 
plitude operators associated with incoming and outgoing light. Unitarity is directly related to 
conservation of photon number during the beam splitter transformation. It is already intuitively 
clear that in the presence of losses, i.e. absorption, photon-number conservation and thus unitar- 
ity cannot be upheld, at least not on the level of the photonic amplitude operators. Having said 
that, because we have constructed a bilinear Hamiltonian of the electromagnetic field even in the 
presence of absorbing dielectrics, Eq. (12241) . there will be a unitary evolution associated with the 
medium-assisted electromagnetic field, but not with the (free) electromagnetic field alone. 

In order to see how the restricted evolution emerges, we consider again a one-dimensional 
model of a beam splitter that consists of a (planarly) multilayered dielectric structure surrounded 
by free space (see Fig. [TJ. As opposed to a mode decomposition in the lossless case, we seek 
the Green function associated with the light scattering at the multilayered stack. In this one- 
dimensional model, the Green function reduces to a scalar function which can be constructed by 
fitting bulk Green functions at the interfaces between regions of piecewise constant permittivity 
[50]. Knowledge of the Green function amounts to knowledge of the transmission, reflection and 
absorption coefficients associated with impinging light of frequency u>. Similar decompositions 
can be made for three-dimensional structures with translational invariance [51]. 

It turns out that the input-output relations d52l have to be amended by a term associated with 
absorption in the beam splitter, 



where gi{ui) denote (bosonic) variables associated with excitations in the dielectric material (de- 



3.2 Light propagation through absorbing dielectric devices 



fa(w) = T(u>) ■ &(w) + A(u) ■ g(w) , 



(257) 



48 



S. Scheel and S.Y. Buhmann 



vice operators) with complex refractive index n(ui) = rj(oj) + i/t(w), and A(ui) the absorption 
matrix. This expression is a direct consequence of the expansion of the electromagnetic field 
operators in terms of the dynamical variables. It is shown in Ref. [50] that the device operators 
are integrated dynamical variables over the beam splitter and read 



d/2 



91,2 M 



2cX±(d, uS) 



^in(uj)ujd/ (2c) 



dx 



-in{uj)u:x I c 



-d/2 



where 



K (u,)ud/c \ sinh[K{u)ud/c] ± sm[rj{u})u;d/c 



k(u>) vi 1 ^) 
The transmission and absorption matrices obey the relation 

T(w) ■ T+(cu) + A(cu) ■ A+{uj) = I 



f{x,u) (258) 



(259) 



(260) 



which serves as the generalisation of the above-mentioned energy conservation relation d48b . 
Equation ( 12601 ) says that the probabilities of a photon being transmitted, reflected or absorbed add 
up to one. Hence, photon numbers and thus energy is conserved only if one includes absorption. 
For a single plate of thickness d surrounded by vacuum, the matrix elements read [33, 50] 



Tii(w) = T 22 (w) 

T 12 (uj) = T 21 (oj) 
Ah(lu) = A 21 (uj) 



— iujd/c 



r(w) 1 - iiHe 2m(,Jj)wd/c L>Hi 2 (w) 
e -^' i/c ii(w)e m(w) " ci/c L»(w)t 2 (w), 



Ai 2 (w) = -A m {u) = v / ^R^e- lwd /( 2 ^ 1 ( w )D( w )yTM 



1 - e 



in(uj)ujd/c 



(to) 



(261) 
(262) 

(263) 
(264) 



The interface reflection and transmission coefficients are functions of the index of refraction and 
are defined as 



n(u>) — 1 
n(uj) + 1 ' 



*i(w) = 



1 + 



*2(W) = 



2n(uo) 
1 + n(cj) 



(265) 



and the factor D(uo) = [1— r 2 (( J S)e 2m ^°' >u " i ' c ]~ 1 accounts for multiple reflections inside the plate. 
These coefficients are special cases of the generalised Fresnel coefficients for p-polarisation and 
normal incidence (see App. |A.4t . 

The input-output relations ( 1257b translate into a generalised quantum-state transformation 
formula. For this purpose, we need to look into an enlarged Hilbert space for the electromagnetic 
field and the dielectric object. With the four-dimensional vectors 6i(uS) — [a(w), g(w)] T and 
f3(u>) = [b([j), h(w)] T , the input-output relations can be extended to a unitary matrix transfor- 
mation of the form 



P(u) = A(tu) ■ 6l(uj) 



(266) 
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where the unitary 4 x 4-matrix A(uj) is an element of the group SU(4) and can be expressed in 
terms of the transmission and absorption matrices as [52] 



AM 



-S(u) 



T(u) 



T(w) C{lo) 



A(u) 
S-^uj) ■ A{uj) 



(267) 



where C(u>) = y/T(oS) ■ T+(cu) and S(u) — y A(lu) ■ A + (uj). Given a density operator g as 
a functional of the input operators a, g- ln = Qi n {ct(uj), or(u)], the transformed density operator 
of the photonic degrees of freedom alone is then [52] 



^ {§ m [A+(w) • a(u), A T H • dt(w)] } 



(268) 



where tr( D ) denotes the trace over the device variables. As a first illustrative example, we con- 
sider the transformation of coherent states at an absorbing beam splitter. If we assume that both 
the incoming electromagnetic field as well as the beam splitter are prepared in two-mode coherent 
states | a) and |/3) with respective amplitudes a and (3, application of the input-output relations 
(I257l i [or equivalently, the quantum-state transformation ( I268H reveals that the outgoing fields 
are prepared in a two-mode coherent state 



lot 



= \T-a + A- (3) 



(269) 



Hence, the transformed amplitudes are determined not only by the transmission matrix T but 
also by the absorption matrix A [33]. 

This theory has wide-ranging applications that include nonclassicality studies of light prop- 
agation through optical elements and entanglement degradation in optical fibres [53,54]. We 
mention here two important results relating to propagation of two-mode quantum states of light 
through optical fibres. Consider a two-mode squeezed vacuum state with squeezing parameter £ 
being sent through identical optical fibres of length I that are held at a temperature T. We regard 
the optical fibres as essentially one-dimensional objects whose effect on the quantum states of 
light propagating through them can be described by the input-output theory presented above. The 
optical fibres are characterised by their absorption length ^ a b s which implies that we approximate 
their transmission coefficient associated with them by |T| 2 = e~'/ iabs . Furthermore, a nonzero 
temperature T gives rise to a mean thermal photon number n t h = [e nu ^ kBT " > — 

As a two-mode squeezed vacuum is a Gaussian state, it is fully characterised by its first and 
second moments, the mean and covariances of its Wigner function or characteristic function, 
respectively. For the quantum state under consideration, the mean is zero and its covariance 
matrix is 



/ 


c 





s 





\ 





c 





— s 






s 





c 







V 





— s 





c 





c = cosh 2£ , 
s = sinh 2£ . 



(270) 



The covariance matrix T transforms under an arbitrary completely positive map as 



r^r' = ArA T + G 



(271) 
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where G is a positive symmetric matrix and A an arbitrary matrix, provided that the resulting 
covariance matrix T' is a valid covariance matrix, i.e. obeys Heisenberg's uncertainty relation 
r' + > [£: symplectic matrix]. The non-orthogonality of the matrix A is a direct con- 
sequence of dissipation, and G is the additional noise as required by the fluctuation-dissipation 
theorem. The general structure (127 U clearly also follows from the application of the quantum- 
state transformation formula (12681 l. 

Using the input-output relations ( 1257b . we find that the entries in the covariance matrix T 
change according to [55] 



c^ c |T| 2 + |i?| 2 + (2n th + l)(l- |T| 2 - |i?| 2 ), s^ s |T| 2 , (272) 
which translates into a transformation of the covariance matrix as 



|T| 2 T + [|i?| 2 + (2n th + 1)(1 - \T\ 2 - |i?| 2 )] I 



(273) 



from which the matrices A and G can be read off as A — \T\I and G = [\R\ 2 + (2n t h + — 
\T\ 2 - \R\ 2 )}I, respectively. 

The entanglement content of a two-mode squeezed vacuum state, expressed in terms of its 
negativity [56], is En ~ 2£. At zero temperature, and neglecting coupling losses into the fi- 
bres, the maximal amount of entanglement (in the limit of infinite initial squeezing) that can be 
transmitted through the fibres is i?Ar imax = — In (l — e~ l / lsbB ). For finite initial squeezing, and 
at finite temperature T, the state becomes separable, i.e. it loses all its entanglement, after the 
separability length 



T ^abs n 
IS = — hl 



1 



2n t h 



(274) 



Note that at zero temperature this separability length is infinite, i.e. in this case it is always 
possible to transmit some entanglement over arbitrary distances. 



3.3 Medium-assisted interaction of the quantised electromagnetic field with atoms 

Up until now, our emphasis has been on devising a quantisation scheme for the electromagnetic 
field in the presence of magnetoelectric background materials, but without external sources. The 
starting point had been the definition of derived electromagnetic quantities such as the displace- 
ment field D and the magnetic field H in terms of the matter-related quantities P and M, the po- 
larisation and magnetisation fields [Eq. (I173H . However, the definition of the latter is not unique 
as parts of them may be associated with magnetoelectric background media, whereas other parts 
could be attributed to additional atomic or molecular sources not included in the background mat- 
ter. It is these additional sources that give rise to a successful theory of microscopic-dielectric 
interfaces. 

In close analogy to the free-space case, the medium-assisted electromagnetic field can be 
coupled to an atomic system consisting of non-relativistic spinless charged particles via the 
minimal-coupling scheme. The dynamics of the combined atom-field system is governed by 
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the Hamiltonian 



00 * 
H = E / d\ J d^%(r,b,).k(r,a,) + X; ^"^"^ 

A— e,m q a 



where 



E 5 T"' - 1 + / rf 3 ^A(r)^(r) (275) 



Pa{t) = y ^q a 5{v-r a ) (276) 



is the charge density of the particles. The vector and scalar potentials A(r) and <j>(r) of the 
medium-assisted electromagnetic field are expressed in terms of the dynamical variables f (r, oj) 
and f*(r, oj) as 

00 00 
A(r) = f doj^E x (r,oj)+h.c. , — V0(r) = J E" (r, + h.c. . (277) 



The total electromagnetic fields are the sums of the medium-assisted fields and the fields 
associated with the atomic system, 

£{r) = E(r)-V^(r), (278) 
B(r) = B(r), (279) 
V(r) = D(r)- £o V^(r), (280) 
H(r) = H(r), (281) 

where 

h(r)= I dV Pf ] (282) 
J 47T£ |r-r'| 

is the scalar potential of the charged particles. In those cases in which the atomic system consists 
of sufficiently localised particles such as in an atom or molecule, it is expedient to introduce 
shifted particle coordinates f a = r a — ya relative to the centre of mass f a = J2 a ( m a/ rn A)i'a 
[with the total mass wia = J2 a m a]- Expanding the vector and scalar potentials A(r) and </>(r) 
around the centre of mass, the Hamiltonian ( 12751 ) simplifies for globally neutral atomic systems 
[<L4 = Ea la = 0] to me electric dipole Hamiltonian 



H = Hp + Ha + Haf 



(283) 
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where 



cm 



H F = d 3 r dLuhw{l(r,u})-{ x (r,Lu), (284) 



o 



' 2m Q ^ 87re r Q - r Q / 
H AF = d • V0(r) I - V -*^p Q • A(f A ) + V A A 2 (f^) . (286) 

I r — r_4 ZiTTLq ZiTTLq^ 



Recall that the electric dipole moment is given by Eq. 

Equations d284t — (1286b could form the starting point for investigations into interaction of the 
medium-assisted electromagnetic field with atomic systems in the long-wavelength approxima- 
tion. However, as for free-space QED, the treatment can be considerably simplified by trans- 
forming using the alternative multipolar coupling scheme. To this end, we again apply a Power- 
Zienau transformation [20-22] 



U = exp 



l - J d 3 rP A (r)-A(r) 



(287) 



where the polarisation is still defined by Eq. d92l ), but the vector potential A is given by 
Eq. ( 1277b in the presence of magnetoelectrics, as opposed to the expansion (f35b valid in free 
space. The unitary operator transformation O' = U OU^ with the operator ( 12871 ) leads to trans- 
formed dynamical variables f^(r, oj) as [57] 

fl(r,o;) = f A (r,a;) + ^ />r'Pi(r') • G+(r,r» (288) 



and thus transformed electric fields ( T99b . The magnetic induction field remains unchanged by 
this transformation, B'(r) = B(r), because it clearly commutes with the operator of the vector 
potential. Other quantities that remain unchanged are those that depend solely on the atomic 
position operators r' a = r a such as the atomic scalar potential, <f> A { v ) — <Pa(y), and the polar- 
isation and magnetisation fields. The atomic momentum operators, however, transform as [57] 

Pa = Pa - 9«A(r a ) - J d 3 r S a (r) x B(r) (289) 
where we have defined the auxiliary vectors 

E Q (r) = ga Q (r)-^V^(r) + ^(r), (290) 

i 

Q (r) = r a J ds s8(r — ta — sr a ) ■ (291) 
" 
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With these preparations, we can now write down the Hamiltonian ( 12751 ) in terms of the trans- 
formed variables as 

oo 



- jd\ P^r) • E'(r) + Y, ^ k + / ^ 2. (r) x B(r) 



(292) 



It should be mentioned that, due to the unitary nature of the Power-Zienau transformation, the 
commutation relations between the transformed atomic position and momentum operators r' a and 
p^, as well as between the transformed dynamical variables f' x (r,u) and fJJ (r,oj) are unchanged. 
The fields E'(r) and B'(r) have to be thought of as being expanded in terms of the transformed 
dynamical variables in the same way as the untransformed fields are expanded in terms of the 
untransformed dynamical variables. 

In the long- wavelength approximation, one replaces the function S(r — f a — sr a ) by its value 
at s = 0, 6(r — ta), so that the polarisation and the auxiliary fields reduce to 

P A (r) = dS(r-i A ) , Q (r) = -f a 6{r-r A ) , E Q (r) = q a & a ( r ) + -^P A (r) . (293) 

2 2m A 



Upon using these simplifications, we obtain the multipolar Hamiltonian ( 12921 ) in long-wavelength 
approximation as 

H = H' F + H A + H' AF 



where 



oo 



X—e.m 



H'a = 



h 'af 



W 2 1 

^ 2m- ^ 



2e 



d 3 rP A (r), 



(294) 

(295) 
(296) 



-d- &(i A ) -&!•&(? A ) +Y 

oc 

3 



8mA 



where 



- / 1 \ "* Qa a 

in - > r Q 

2 ^ m a 



d x B'(r A ) 



x Po 



8m 

a 

^- • [d x B'(r A ) 
mA L 



r Q x B'(r A ) 



(297) 



(298) 



is the magnetic dipole operator, = ^ a p' a , and t>a = Pa K/fii)^. 

At the moment, the magnetic dipole moment operator comprises only the contribution from 
the angular momentum. Spin can be included via a Pauli interaction term in Eq. ( 12751 ) leading to 



Qa 



2m, 



-r a x p a + j a s c 



(299) 
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(s a : particle spin, j a : gyromagnetic ratio) [58]. 

The scheme can be easily extended to include more than one atomic subensemble. In the 
minimal-coupling scheme, this leads to interatomic Coulomb interactions in Eq. ( 12751 ) which is 
why the multipolar-coupling scheme is strongly preferable. In this case, Eq. ( 12941 ) generalises 
to [130] 



H> - H' F + ^ [ h 'a + h 'af 



(300) 



where the dynamics of each atom is given by a Hamiltonian of the form ( 12961 ) and each atom 
couples individually to the electromagnetic field via coupling Hamiltonians of the form ( 1297b . 

As in the free-space case, we will henceforth treat all atom-field couplings within the frame- 
work of the multipolar coupling scheme and drop all primes denoting multipolar variables. 



3.4 Nonlinear quantum electrodynamics 

In all previous (as well as all subsequent) sections we have concentrated on magnetoelectric 
background materials whose dielectric (and magnetic) response to an external perturbation can be 
described within the framework of linear-response theory. That is, polarisation and magnetisation 
fields are linearly and causally related to the primary electromagnetic fields. For purely dielectric 
media, this means that [cf. Eq. dl74H 

oo 

P(r,f) = e J dr X (r,r)E(r,i-r)+P N (r,<) (301) 
o 

where x(r, r) is the dielectric susceptibility. However, many materials show nonlinear behaviour, 
i.e. their dielectric response has to be described by a polarisation with a component that depends 
quadratically (cubically, quartically etc.) on the external electric field. In free space where 
dispersion and absorption are disregarded and a mode expansions of the electromagnetic field 
can be used, the effect of these nonlinear polarisations is to add interaction Hamiltonians that 
are cubic (quartic, quintic etc.) in the photonic amplitude operators (see, e.g. [59]). These 
effective Hamiltonians arise from off -resonant interactions with atomic systems in rotating-wave 
approximation [see, e.g. Eq. ( 11431 )1. Our aim is to extend this theory to nonlinear processes in 
absorbing matter where mode expansions generally do not hold. 

For this purpose, we consider an interaction Hamiltonian in its most general normal-ordered 
form corresponding to a medium. Using the abbreviation k as a short-hand for the collection 
of spatial and frequency variables k = (rfc,Wfe), this interaction Hamiltonian can be written 
as [60] 

i? NL = f dld2d3a ijk (l,2,3)fj(l)f j (2)f k (3) +h.c. (302) 
with an as yet unknown coupling tensor a^(l, 2, 3) that will eventually have to be (linearly) 

(2) 

related to the second-order nonlinear susceptibility x\j k - The integration initially ranges over all 
frequencies, even including those frequencies that would not guarantee energy conservation. 
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The Hamiltonian ( 1302b . together with the Hamiltonian H^ of the linear theory (we include 
the index L here to distinguish it from the nonlinear interaction Hamiltonian), is now used to 

construct the time-dependent Maxwell equations such as Faraday's law V x E(r) = — B(r) as 



V x E(r) = - — 

in 



B(r),ff L +ff, 



NL 



(303) 



The fields E(r) and B(r) have to be thought of as being expanded in terms of the dynamical 
variables f (r, u) and ft (r, oS). Since Faraday's law is valid irrespective of the presence of matter, 

and by definition [B(r), Hi]/ (ift) = B(r), we must have 



B(r),ff 



NL 



(304) 



which is a condition that is needed in the next step. We write Ampere's law V x H(r) = D(r), 
using Faraday's law, as 



V x V x E(r) = -MoD(r) = -/x D L (r) - // P N L(r) 
where we split up the total displacement field D(r) into its linear part, 



(305) 



D L (r) =e e(r,w)E(r) 



(306) 



and some nonlinear polarisation Pnl (r). Heisenberg's equations of motion then imply that 



VxVxE(r) = g 



Mo 

~ K 2 



D L (r),H L 
D L (r),i? NL 



Mo 

n 2 



t> L (r),H L 



NL 



Mo 

h 2 



P NL (r),i?i 



(307) 



where we kept only those terms that are at most linear in the coupling tensor cty fe. The terms 
that have been left out have to be included into higher order nonlinear processes. The first term 
on the rhs of Eq. ( 13071 ) is by definition equal to the lhs of the same equation. The second term on 
its rhs vanishes because of the constraint ( 13041 ). The remaining two terms have to satisfy 



E>L (!■),£ 



NL 



pN L (r),ff L 



(308) 



which yields a solution for the nonlinear polarisation Pnl(i"). Note that the double commutator 
has been reduced to a single commutator as it turns out that a general solution would have to 
include functionals that commute with whose contributions can be shown to diverge [61] and 
thus have to be discarded. 

The way to solve Eq. ( 13081 ) is to view its rhs as being the Liouvillian generated by £Zl, i.e. 
£l» = Hi], whose inverse can be formally written as 



PN L (r) = --£L 1 



1 

ih 



D L (r),#] 



NL 



(309) 
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The action of the inverse Liouvillian on an operator O is given by 

oo 

fc£6 = Urn j dte- st e^ t / h de i "^ t/h (310) 
o 

which can be checked by direct calculation, with the result that 

oo 

PnlW = -4 lim f dte- st e- i6 ^ H \t) L (r),H Nh \ e l "^ t/h . (311) 
ih s^o J L J 



Before we continue solving this equation, we remark that due to the decomposition of the lin- 
ear displacement field ( 13061 > into a reactive part and a Langevin noise contribution Pi,(r), the 
nonlinear polarisation also contains a contribution, pj^ (r)= — (ifi,) -1 /^ 1 [p£ ( r )i ^l]> triat 
disappears identically with vanishing absorption and can thus be regarded as the nonlinear noise 
polarisation. The result of the Liouvillian inversion can be cast into the form [60,61] 

P N L,m(r) = / d0 d2d3 ^^-a njk (0, 2, 3)^U(r, w)G m „(r, 0)/ i (2)/ fc (3) 

+ i r d0dld3 V Im£ (° ) a: ( 1;0;3) ^ £(r ^ )Gmn ( r;0) /f( 3) / a(1) 

in V 7r J u>i — u>3 c l 



-h.c + <' m (r) (312) 



with 



in \ IT J L02 + ^3 

+^V^ / d0dld3 Vlm£(0) «L fc (l ; 0,3) ( 5(r-s)/t(3),A(l) 
in \ n J lu\ — LO3 

+h.c. (313) 

and the notation = (s, w). 

From Eqs. ( 1312b and (13131 ) it is hard to see how the coupling tensor ctijk has to be related to 
the nonlinear susceptibility. Instead, we invoke comparison with the definition of the nonlinear 
polarisation from classical nonlinear response theory [59], 

t 

P NL ,™(M)=£o / dT'dr"x^ s (r J t-r , ,i-r")^r(r,T')-B s (r ) r")+P£ ) m (r,i) (314) 

—00 

and introduce slowly varying electric fields whose (non-overlapping) amplitudes are centred at 
the mid-frequencies flj with £1 = fij + £1 2 such that 

P N L,m(r, n,t) = eox£Ur, ^o, "i, ^)S r (r, tl^E^r, fi 2 ) + P$} m (r, tl ,t) . (315) 
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Expressing the electric fields in Eq. ( 1315b in terms of slowly varying dynamical variables and 
comparing with the slowly varying version of Eq. ( 1312b yields the sought after relation between 
a ijk and x fjl as [60,61] 



/ o x ft /Ime(s 2 ,fi 2 )lme(s3,fi 3 ) 

aijfc(r, J2 , s 2 , S2 2 , s 3 , S2 3 J = — 



xif M (r,fi ) 



Ime(r 2 ,^o) 

xil(r,n2,n 3 ) Gm . (ri S2) ^ 2)Gnfe(rj S3; 



(316) 



where Hu(r, Qo) are the cartesian components of the Helmholtz operator, or equivalently, the 
integral operator associated with the inverse dyadic Green function. 

Reinserted into the nonlinear interaction Hamiltonian ( 13021 ). and combined with the dynami- 
cal variables, Eq. ( 1316b will yield a formulation in terms of electromagnetic field operators. What 
becomes immediately clear, though, is that i?NL will not be of the form i?NL oc x\- k ElEjEk 
as in standard nonlinear optics. Mathematically, the reason is that there is no third Green tensor 
in Eq. ( 13161 ). Instead, an inverse Green function (or Helmholtz operator) has to be dealt with 
which, by formally expanding it into a power series, will lead to additional contributions to the 
standard nonlinear interaction. The physical reason for this behaviour has to be sought in the 
fact that, from a microscopic point of view, an effective nonlinear interaction does not take place 
in free space but rather inside the absorbing medium where the local electric field is altered by 
the presence of the dielectric material. To account for that, local-field corrections such as those 
discussed in Sec. lA.6l have to be included which automatically yields additional contributions to 
the nonlinear interaction Hamiltonian that are not of the standard form. 
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4 Atomic relaxation rates 



The first set of applications we consider in this review regards the theory of atomic transition 
rates. In this section we discuss the influence of dielectric bodies towards atomic relaxation and 
heating rates, and some of their experimental ramifications. 

To begin, it is necessary to study the dynamics of internal atomic degrees of freedom in 
the presence of absorbing magnetoelectric matter. In Sec. 13.31 we derived the Hamiltonian of 
the system comprising the medium-assisted electromagnetic field, the atomic system and their 
mutual interaction (here taken in the electric dipole approximation) as 



H = Hp + Ha + Haf 



with 



oo 

H F = J <$ T J dwhjji\(T,w) -f A (r, 



(317) 



(318) 



A— e,m 



the Hamiltonian of the medium-assisted electromagnetic field [Eq. ( I295H in the presence of 
magnetoelectric bodies. The free atomic Hamiltonian Ha [Eq. d296H is expanded into atomic 
energy eigenstates \n), 



(319) 



with the corresponding eigenenergies huj n , and the atomic flip operators A T 
the commutation rules 



A M ,A n 



— 3lmA kn 3knA m l . 



\m)(n\ obeying 



(320) 



The electric -dipole interaction Hamiltonian Haf which, in multipolar coupling, is given by 
Haf = — d • E(r^) [first term in Eq. ( [297) 1. is expanded in terms of the energy eigenstates as 



Haf = - ^ d mn • ^(ta)A\ 



(321) 



where d mn — (m|d|n) are the matrix elements of the dipole operator. Recall that we have 
dropped all primes on the electromagnetic field operators that indicate the multipolar coupling. 

In similar fashion to the free-space theory outlined in Sec. 12.2.11 the internal atomic dynamics 
is governed by the solution to the coupled set of Heisenberg's equations of motion (u> mn = 



A,. 



h L 



mn- r *-mn 



| ^ J duj (d nk A mk - d km A k7 }j ■ E(r A , of) 

k q 

+E f (rA, to) ■ (dnkAmk - d km A kn ^j 



(322) 
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and 



fx(r,o;) = - 



(323) 



fx(r, w),H = -iwf\(T,u) + -^d mn • G* x (r A ,r,uj)A mn . 
We formally integrate Eq. d323t as 

t 

fx(r,w,t) = e- iwt f A (r,o;) + / dr e'^^d™ • GJ(r A , r, w)i ron (r) (324) 



and reinsert this formal solution into Eq. ( 1322b and obtain 



oo 



+e 4 "*E t (r A ,u;) • d„ fc i mfe (£) - d fcm i fen (i) I + Z mn (i) 



E(r A ,w) 
(325) 



where 



2- / dijjuj 1 \ dr 

M.J o 



e-^-^A^A^T) - e^-^A^A^it) 
i-^-^Ak^Auir) - e^-^A^A^t) 



d nk ■ lmG(r A ,r A ,uj) ■ dy 



d kn ■ lmG(r A ,r A ,uj) • dy 

(326) 



is the zero-point contribution to the internal atomic dynamics due to the second term in Eq. (1324b 
A self-consistent solution is obtained if Eq. ( 1325b is formally integrated as 



oo E 

f i m „(0) + ^^ Jdujdrt 



x < e 



11 o o 
d„fcA mfc (r) - d m fcA fen (T) •Efr^o;) 



+e' i " r E t (r A ,w) • d„ fc A mfe (r) - d mfe A fcn (r) 



(327) 



reinserted into Eq. (1325b . and (thermal) expectation values being taken. The u mn are the shifted 
atomic transition frequencies. Using the thermal expectation values for the electromagnetic field 
operators given in Sec. 13.1.21 Eqs. ( I242b -( I245I ). what remains is a set of coupled differential 
equations for the atomic quantities, 



(A mn)T — 2LJ mn (j4 mn )7 1 4~ (Z mn )x "h (^mn} r J 



(328) 
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e -iu(t-r) , e iu(t-r) 



where the thermal contributions (f mn ) T read (n th = [e^l^^ - l]" 1 ) 

OO t 

(T mn ) T = f duu 2f hh(u) J dr 

k,l g 

(A m i) T d nk ■ JmG(r A ,r A ,Lu) ■ d k i - {Ai k ) T d nk ■ Im.G(r A ,r A ,cu) ■ d t 
(A k i) T d km ■ lmG( ■ ImG^r ^,1^4,0;) • d /fc 



x < e 



i£fe„(t— r) 



(329) 



recall Eq. (13261 >. These sets of equations form the basis for all following investigations in internal 
atomic dynamics. 



4.1 Modified spontaneous decay and body-induced Lamb shift, focaf -field corrections 

For weak atom-field coupling, these expressions can be further evaluated using the Markov ap- 
proximation (see Sec. 12.2.11 in which the expectation values at times t can be related to those at 
the upper limit t of the time integrals as 



(A mn (r))T * e-^»(*-^(i mn (t)) T • 
The time integrals themselves are approximated as in the free-space theory by 



J dr e^-^XH ~ tt8(uj - u mn ) + iV 



1 



(330) 



(331) 



Introducing the atomic density matrix elements <r mn — (m\&\n) = (A nm ), the set of differ- 
ential equations reduces to 



-r n <7 nn (i) + X r fe „CT fcfc (i) , 



iuj rnn (T m -|- T n ) 



(t) , m^n. 



(332) 



(333) 



These equations resemble closely those obtained in free-space theory. One can identify indi- 
vidual decay rates T nk from state |n) to \k), and total loss rates T n — J^k^nk of a level |n). 
The individual rates contain zero-point [superscript (Z)] as well as thermal [superscript (T)] 
contributions, T nk = T$ + T^}, with 



T< nk = ^^fcB(w„ fc )d nfe ■ lmG(r A ,r A ,Q nk ) • d kn , 
^nk = -j^-^lk^nk ■ JmG(r A ,r A , \& n k\) ■ dfen [S(w„fe)nth(w„fe) + 6(w fe „)n t h(wfc r , 



(334) 



(335) 
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Similarly, the shifted atomic transition frequencies ui mn = u) mn + 5u> m — 5u n depend on the 
frequency shifts 5ui m = J\ ^nk where the zero-point and thermal contributions to the level 
shift induced by a single level \k) read 



8u> 



(Z) 
nk 



/ <W 

TTh J 





d nk ■ ImG^(r A ,r A ,u}) ■ d fai 
ui n k — LO 



Soj( nk = I diouj 2 d nk ■ lmG{r A ,v A , lo) ■ d kn 

TTtl . 



nth(w) nth{uj) 



l^nk — U> ^nk + 



(336) 



(337) 



In Scol^}, we have explicitly used the scattering part G^ s \r A , r A , uj) of the dyadic Green func- 
tion because we assume that the vacuum-induced Lamb shift has already been included into the 
bare atomic transition frequencies io mn . 



Spontaneous decay of a two-level atom: As an instructive example, let us consider a two- 
level atom with energy levels \g) and |e) separated by an energy hui A — hio eg which we assume 
to contain the free-space Lamb shift. At zero temperature, the spontaneous decay rate T = 
of the excited state |e) is then given by Eq. ( I334t as [62,63] 



2uj 2 

r= ^d-ImG(r A ,r il ^).d* (338) 

where we have neglected the shift in the atomic transition frequency lo a . This is the same result 
one would obtain using perturbation theory, i.e. Fermi's Golden Rule. We have noted previously 
that the rate of spontaneous decay in free space is proportional to the strength of the vacuum 
fluctuations of the electric field [Eq. dl 191 )1 which, by Eq. ( 12451 . is now seen to be proportional 
to the imaginary part of the dyadic Green function, 

<0|E(r A ,w)®Et (1^,^)10) = —^$-ImG(r A ,r A ,u A )5(uj-u A ), (339) 

7T£o cr 

which underpins our interpretation of w 2 Im G(r A , r A , lo a ) as the local density of states. The rate 
T of spontaneous decay in vacuum, Eq. fll201 >, is recovered by inserting the free-space Green 
tensor G^ (r, r', u) into Eq. (1338b . 

To give a nontrivial example of how Eq. d3381 > can be used for investigating atom-surface 
interactions, we consider a two-level atom placed near a planar dielectric half-space with permit- 
tivity e(ui). The Green function for this structure is known (see Appendix lAb. and the result in 
the limit z a oj/c<^. 1 reads [64,65] 

r = T 4 f 1 + w) f -^) 3 ,! m 1? 1 ) |2 + OizA 1 ) (340) 
8 V |d| / \u! A z A J \e(uj A ) + l\ 2 



Hence, a dipole oriented perpendicular to a planar surface decays twice as fast as a dipole parallel 
to it. Note the difference to the spontaneous decay rate near a perfect mirror which approaches 
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finite values (2ro for perpendicular dipole orientation, for parallel dipole orientation) in the on- 
surface limit. In reality, however, nonradiative decay processes cause the near-field spontaneous 
decay rate to diverge. 

Another important example, which has been used to investigate modified spontaneous decay 
inside a dielectric host medium [33,63], is that of an atom in a spherical microcavity of radius 
Rcav (Fig. n~TT >. We use the dyadic Green function presented in Appendix lAl and note that, if 
the atom is located at the centre of the microcavity, only the TM-wave spherical vector wave 
functions N pm i(fc) do not vanish. Inserting the result into Eq. ( 13381 > yields [33,63] 

r = r [l + Rerf {1 = 1)] (341) 

where the reflection coefficient r 22 for I = 1 takes the form 

r 22 (; = x) = [j + Pin + 1) - W 2 n - p 3 n 2 /{n + l)]e l P 

p sin p — p(cos p + in sin p) + ip 2 n cos p — p 3 (cos p — in sin p)n 2 j(n 2 — 1) 

(342) 

[n = n(uj) = \J e(u>) and p = R cw u>a/c]. In the near-field limit, i.e. when the size of the 
microcavity is much smaller than the atomic transition wavelength, we can expand in powers of 
P = RchvOJa/c <C 1 and obtain [see also Eq. ( I595H 

^ T f3(e-l)/ c \ 3 9[4e 2 -3e-l]/ c \ . 9e 2 n ^. .1 

r ^ roIm {w(^j + 5( 2£ + i) 2 \ ^) +1 wtw + 0{r - ] ) 

(343) 

[e = e(uj A ), n = ^e{uj A )]- 

The leading terms in both examples, Eqs. ( 13401 ) and ( 13431 l. are proportional to the inverse 
cube of the atom-surface distance. This is attributed to resonant energy transfer to the absorb- 
ing dielectric surroundings which is a nonradiative decay process [65-69]. The next-to-leading 
terms, the induction terms, are proportional to the inverse atom-surface distance and correspond 
to absorption of real photons. The nonradiative decay rates make it impossible to hold atoms or 
molecules near dielectric or metallic surfaces and use their internal states for coherent manip- 
ulation. The correct quantum-statistical description of the dielectric bodies also imply that the 
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suppression of spontaneous decay of a dipole parallel to a mirror surface (see Sec. 12.2.2b cannot 
be observed. 

As the near-field contributions are proportional to the imaginary part of the permittivity, 
lme(uJA), these terms will vanish if absorption can be disregarded at the atomic transition fre- 
quency u A . In this case, the spontaneous decay rate is modified to [63,70] 



r 



3e(uj A ) 
2e(u A ) + 1 



n(u> A ) 



(344) 



The results for a small spherical cavity can also be used to study the decay of atoms embedded 
inside a medium. In this real-cavity model, the cavity implements the local-field correction 
arising due to the difference between the macroscopic field and the local field experienced by the 
atom. Using the local-field corrected Green tensor from Appendix I A. 61 one can show that, for 
weakly absorbing media, Eq. (1344b remains valid for arbitrary geometries when written in the 
form 



Tloc 



/ 3e(uj A ) 
\2e(uj A ) + l 



(345) 



where T is the uncorrected decay rate. For absorbing media, the more general result 



r = rv 



c 2 heo 



d • Im 



3e{u A ) 
2e(lu a ) + 1 



G {s \r Al r A ,u A ) 



■A* 



(346) 



holds, where Tc is the near-field rate (13431) [71] 



4.2 Sum rules for the local density of states 

One of the principal cornerstones of macroscopic QED is the validity of Kramers-Kronig rela- 
tions (1178b for the response functions such as the dielectric permittivity, the magnetic permeabil- 
ity or the generalised conductivity. Equally important are sum rules that follow from asymptotic 
limits of these Hilbert transforms. Examples for sum rules in optics are the optical theorem 
that relates the imaginary part of the forward scattering amplitude to the total scattering cross 
section [2], or the Thomas-Reiche-Kuhn oscillator strength rule [72]. A particularly relevant 
relation for our purposes is the integral relation [73] 

oo 

J dcu [r)(u) - 1] = (347) 
o 

satisfied by the real part t](lu) of the complex index of refraction n(uj) — \J e(io)fi(uj) = rj(ui) + 
in(uj). This follows from the superconvergence theorem for Hilbert transform pairs [74]. It 
means that the propagation properties of light in a dielectric medium are redistributed in such 
a way that, averaged over the whole frequency axis, the (real part of the) refractive index is the 
same as in vacuum. 

Because macroscopic QED is inherently based on linear response theories, it is conceivable 
that quantities such as spontaneous decay rates would obey certain sum rules, too. Recall from 
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Sec. 12.2.21 that the rate of spontaneous decay and the Lamb shift form a Hilbert transform pair. 
We will now try to derive a sum rule for the rate of spontaneous decay, and investigate integrals 
of the form 



oo 

r — Tq f , 6irc 



duj—^= diu—^d-lmG^HrA^A^)-^, (348) 
To J war 



where To is the free-space decay rate ( 11201 ) and (ta , ?a , the scattering part of the dyadic 
Green function. Hence, we are seeking to compute the integral 

oo oo 

J duj^lmG is \rA,r A ,iu) = Im J dw^G^(r A ,r A ,w) (349) 

-oo 

where we have used the Schwarz reflection principle (11931 . 

For the sake of definiteness, we imagine the radiating atom being placed inside a dielectric 
structure such as the spherical microcavity sketched in Fig. QT| In fact, any generic situation 
will involve an atom in free space surrounded by some dielectric material, although not neces- 
sarily in a spherically symmetric way. The scattering part of the dyadic Green function for any 
arrangement of dielectric bodies can always be expanded into a Born series (see Appendix lA.5b 
as 

G^(v a ,t a ,uj) = [ d 3 S ' X (s',Lu)G( Q \rA,s',Lu)-G( \s',r A ,Lo) 

c J Vx 

d 3 s'd 3 s" X (s', u) X (a", U )Gf°)( ri , s', u) ■ G<«V, ■»", ■ G (0) (s", r A ,u) 
+ ... (350) 

where G(°'(r, r',o;) is the free-space Green tensor and x(r,o;) the dielectric susceptibility of 
the material surrounding the radiating atom. The integration extends over the total volume of the 
body, but excludes the location of the atom. The free-space Green tensor can be read off from 
Eqs. ( 1534b and (1535b . setting q(uj) = ui/c. The first thing to note is that all terms containing the 
5 function involving the location ya of the atom do not contribute to the Born series. 

As a function of u), the free-space Green tensor has single and double poles at co = whose 
contributions to the integral (1349b give rise to contributions that can be computed as follows. 
Concentrating on the first term in the Born series expansion, we need to look at contributions of 
the form f(u>) = J2 n c «x( s ': u})e 2lup / c /uj n with n = 1, 2, 3 whose residues at uj = are [75] 



n-1 



Res/HU = £ 7 -!* J2 ^x^—V.O) . (351) 
^— ' (n — 1)! ' ml 

n—l v 1 m—0 

These terms arise from a short-distance or, equivalently, low-frequency expansion of the Green 
tensor. It is seen from Eq. ( 1534b that they can be traced back to the longitudinal part of the 
free-space Green tensor. Hence, all these pole contributions, in any order of the Born series, 
arise from either purely longitudinal terms oc tJ^G^- '" (r A , s', ui) ■ ■ ■ G' -"' (s" +1 ) , r A , ui) or 
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from those in which one (and only one) of the longitudinal Green tensors has been replaced by 
a transverse part G^°' >± (s^ , s(' i+1 ) , uj) [75]. Subtracting these terms from the total scattering 
Green tensor, 

G {s) '(r A ,v A ,u) = G^ s \v A ,v A ,w) 

■ f d s s'x(s',co)G^\(r A ,s',co)-G^\( S ',v A ,uj) 
Jv x 

■ [ d 3 S ' X (s'^)G( )ll(r A ,s'^).G(°) ± (s',r A , W ) 

~ f d 8 8 ' x ( g »G0 1 (r yl)S ' |W ).G™l(s' I r A|W ) 
c Jv v 

d 3 s'd 3 s" ■■■ , (352) 

Vv 



^ 2 



c 2 



or 



c 2 



2\ 2 
0T 



and noting that all other contributions containing non-negative powers of oj vanish after contour 
integration, it becomes obvious that 

oo 

J du^lmG {s) ' {r A ,r a ,lo) = 0. (353) 
o 

This in turn means that we can define a modified spontaneous decay rate V which is constructed 
from G^ (r A , r Al uj), that obeys the sum rule 

oo 

p/ p 

dui— — -= 0. (354) 
1 o 



From its construction, it is clear that V excludes dipole-dipole interactions. Because nonradiative 
contributions to the decay rate have been subtracted from T, we can also interpret the sum rule 
(13541 ) as a conservation of the integrated local density of states associated with photonic final 
states. Recalling our discussion following Eq. ( 1347b this means that a magnetoelectric medium 
merely redistributes the photonic density of states across the frequency axis, but remain at the 
free-space level on average. 

Finally, we make the connection to local-field corrections. The local-field corrected single- 
point Green tensor, obtained after discarding the terms containing i? C av in Eq. ( 1595b , is one 
example of a modified Green tensor ( 13521 ) that leads to a valid sum rule for spontaneous decay 
rates. In this way, a connection is established to the theory in Ref. [76] which is strictly valid 
only in nonabsorbing materials. 



4.3 Heating of polar molecules 

In a previous section we argued that at very small atom-body distances, nonradiative decay will 
dominate. Given that this would require the atom or molecule to be held at distances much 
smaller than the (optical) wavelength, this effect will hardly be seen in a controlled experiment. 
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Spec. 


£ c (GHz) 


w c (THz) 


/z e (10- 3U Cm) 


M'e(lO-^C) 


T r ( S ) 


T V (S) 


LiH 


222 


42.1 


19.6 


60.5 


2.1 


25 


CaF 


10.5 


18.4 


10.2 


172 


3,400 


4.7 


BaF 


6.30 


14.1 


11.7 


285 


7,200 


1.8 


YbF 


7.20 


15.2 


13.1 


195 


4,400 


4.1 


LiRb 


6.60 


5.55 


13.5 


21.4 


4,900 


128 


NaRb 


2.03 


3.21 


11.7 


12.6 


70,000 


1,400 


KRb 


1.15 


2.26 


0.667 


1.89 


6.7xl0 7 


120,000 


LiCs 


5.80 


4.92 


21.0 


28.4 


2,600 


80 


NaCs 


17.7 


2.94 


19.5 


21.4 


330 


580 


KCs 


92.8 


1.98 


8.61 


6.93 


62 


12,000 


RbCs 


0.498 


1.48 


7.97 


4.41 


2.5 xlO 6 


63,000 



Tab. 3. Properties of various diatomic radicals (electronic ground state, rotation and vibration constants, 
dipole moment and its derivative at equilibrium bond length, reduced mass) and life times r r and r v of their 
rovibrational ground states against rotational and vibrational heating at room temperature (T — 293K) in 
free space [78]. 

The situation is changed dramatically when the transition wavelength in question is very large, 
say longer than a centimetre. This is the case for polar molecules since the spacings between 
neighbouring rotational and vibrational levels are relatively large. Molecules whose projection A 
of the total orbital angular momentum L vanishes (A = 0) are best described by Hund's coupling 
case (b) [77]. In this scheme, the molecular eigenstates IS", N, J, M)\v) are characterised by the 
quantum numbers J and M of the total angular momentum and its projection on the space-fixed 
x-axis, the total spin quantum number S (S), the rotational quantum number N (N = J — S) 
and the vibrational quantum number v. For deeply bound states, the rotational and vibrational 
eigenenergies are [77] 

E N = hB e N(N+l), N = 0,1,... E v =hw e (v + ±),v = 0,l,... (355) 

where B c and uj c are the rotational and vibrational constants, typical values of which are listed 
in Tab. [3] 

Even at room temperature, the first few rotationally and vibrationally excited states may 
be considerably populated due to heating out of the ground state, placing severe limits on the 
coherent manipulation of polar molecules in their ground states. According to Eqs. (1335K the 
respective heating rate is given by 

T = ^ Wfc nth(wfeo)d fe • Im G(r A , r A , uj k o) • d k0 , (356) 

k 

where frequency shifts [Eqs. ( 1336b and ( I337H can typically be neglected. For a molecule in free 
space, use of the Green tensor (1536b leads to [78, 79] 

k k u 
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In order to evaluate these rates, the relevant dipole matrix elements need to be determined. For 
rotational transitions, this can best be done using Hund's case (a) basis \S, A, S, il, J, M), where 
A, S and denote the projections of L, S and J onto the internuclear axis (f2 = A + £). In this 
basis, one has [77] 



d mn = (njM\d\n'j'M') = Li e (njM\u\n'j'M') 



/'c 



dr 1 - v +1 )— + (u^ 1 + u +1 + u° e 

V mn 771 \ win 1 mn) win z 



(358) 



where fi c is the molecular dipole moment at the equilibrium internuclear separation, u = f/|f | 
and 



Cn = (-1) A/ -V(2J+1)(2J' + 1) ( 



J 1 J'\( J I J' 

—M q M'J \-n fl' 



(359) 



Relating the two bases via 

s 



S,N,J,M) = J2 (-1) J " S V27VTT ( ^ S n ^)\n,J,M) (360) 



one finds 



dofc ® d fe o = lull (361) 

fe 

for the molecules under consideration. For (ro-)vibrational heating, the dipole matrix elements 
follow from [77] 

(vQJM\d\v'n'JM') = ^' c {nj M\u\VL' J' M'){v\q\v') , (362) 



<-» M >^t (363) 



(fi' c : derivative of the dipole moment at equilibrium bond length, m: reduced mass) to be 

^d ot «d M = ^ /. (364) 



127rmu; c 

fe 

Combining the above results, one can calculate life times tq = T^ 1 of polar molecules in 
free space against rotational and vibrational heating out of the ground state at room temperature 
(T = 293K) which are listed in Tab. [3] Since the rotational transition frequencies lie typically 
well below the maximum of the thermal spectrum (17THz at T = 293K), rotational heating 
mostly affects light molecules like LiH whose transition frequency is largest. For this molecule, 
rotational heating severely limits the life time of the ground state to about 2 seconds. Vibrational 
heating, on the contrary, mostly affects the fluorides whose vibrational transition frequencies are 
very close to the peak of the thermal spectrum; associated lifetimes lie in the range of a few sec- 
onds only. The rotational and vibrational heating rates exhibit strong temperature-dependences 
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via the thermal photon number, so the impact of both heating channels can be considerably re- 
duced by lowering the environment temperature. 

Just like spontaneous decay, heating can be considerably enhanced when molecules are 
placed close to surfaces. Using the decomposition G = G^ + G^ s > of the Green tensor into 
its free-space part and the scattering part accounting for the presence of the surface and recalling 
Eqs. (1361b and ( 1364b . the ground state heating rate ( 13571 ) can be written as [78] 



r(^) = r 



i 



2nc 

WfeO 



Im (r A ,r A ,uj k o) 



(365) 



(za- atom-surface separation), where G^ is given in App. IA.4I Note that due to our neglect 
of the frequency shifts, Eqs. d336t and d337| >, the heating rate separates into the temperature- 
dependent factor T as given by Eq. ( 13571 ) and a purely position-dependent part. 

In the nonretarded limit, za \ \/ £{uj n k)\u>nk /c <C 1, the position-dependence is approximately 
given by [78] 



r(*u) 







1 



Ime(wfco) 



For metals with Drude permittivity 

, ,2 



1 



uj(uj + £7) 

and for sufficiently small transition frequencies ujko *C 7 < <^p, one has 



7 



2w P W fcO 



(366) 



(367) 



(368) 



and the condition za | y/s(u) n k)\uj n k /c <C 1 is not even valid for very small atom-surface separa- 
tions. Instead, the heating rate is well approximated by the empirical formula [78] 



n:,> = r„( i + jf + f| 

Z A Z A 



3c 

T 



2tJ P W fcO 



(369) 



which is valid for distances za < z c (z c being the critical distance for which surface-induced 
heating becomes comparable to free-space heating). 

In the opposite retarded limit, ZAUJ n k/c ^> 1, the ground-state heating rate is approximately 
given by [78] 



= To 



1 



1 



ZzA^kO \ 1 



sin 



2zA^nk 



(370) 



where the second equality holds for good conductors. The distance-dependence is thus governed 
by attenuated oscillations away from the surface where the oscillation period equals twice the 
molecular transition wavelength. 




Fig. 12. Heating rates for CaF and NaCs as a function of distance from a gold surface (uip = 1.37 x 
10 16 rad/s, 7 = 4.12x 10 13 rad/s). Solid lines: total heating rate. Dotted lines: vibrational excitation rate. 
Dashed lines: rotational excitation rate [78]. 



The spatial dependence of the molecular heating rates over the entire distance rates is shown 
in Fig. [12] for CaF and NaCs. It is seen that vibrational heating dominates for CaF and results 
in rapid oscillations of the heating rate as a function of distance. For short distances, rotational 
heating strongly increases and begins to contribute to the total heating rate. For NaCs, rotational 
heating slightly dominates for moderate distances, although the oscillations associated with vi- 
brational heating are still manifest in the total heating rate. At distance smaller than about 10/im, 
rotational heating becomes strongly dominant. The results show that molecular heating can 
strongly increase in close proximity to surfaces, thus placing severe limits on the miniaturisation 
of molecular traps. 
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L_l?j substrate 
I77 transverse 




Fig. 13. Schematic set-up of a wire trap creating a confining trapping potential for low-field seeking atoms 
(left figure) (picture taken from Ref. [85]). Typical experimental set-up using a Z-shaped wire and a reflec- 
tive gold surface (right figure) [picture courtesy of E.A. Hinds]. 



4.4 Spin-flip rates 

In cold-atom physics, where microengineered magnetoelectric or metallic structures are de- 
signed to magnetically trap ultracold atoms (e.g. 87 Rb) in a well-defined Zeeman sublevel of 
their respective hyperfine ground states (cf. Fig. [13}, typical transition frequencies range from 
100 kHz ... 10 MHz (for reviews, see [80-82]). The equivalent free-space wavelengths are thus 
on the order of several metres. As typical atom-surface distances are in the sub-millimetre range, 
the atoms are in the deep near field of the relevant (magnetic) transitions. 

In complete analogy to electric-dipole transitions, we can investigate magnetic -dipole transi- 
tions using the approximate interaction Hamiltonian ( |297t 



Hai 



-m-B(r A ). 



(371) 



From the statistical properties of the electromagnetic field and our previously used arguments 
regarding Fermi's Golden Rule we already know that the spin transition rate is proportional 
to the strength of the (thermal) magnetic field fluctuations, T oc (B(r^, lu)® B^r^, oj a ))t- 
If we assume an atom to be in its electronic ground state, the magnetic moment vector rh is 
proportional to the electronic spin operator S (the nuclear spin operator I is smaller by a factor 
of m e /m p , i.e. the ratio between electron and proton mass; we also assume that the ground state 
has L = 0). It can then be shown that the transition rate between two magnetic sublevels |i) and 
I/} is [84] 



2mo 



m ■ Im 



V x G{r A ,r A ,uj A ) x V 



[nth(ku) + 1] 



(372) 



(/|S|i)-Im VxG(t a ,t a ,ua)xV' ■ <i|S|/> [nthM + 1] 



[m = (/|m|i), fiB- Bohr magneton, g$ ~ 2], Note that this rate corresponds to the rates associ- 
ated with electric transitions as given by Eqs. ( 13341 ), ([335} by means of a duality transformation 
d <-> m/c, e <-> /i, cf. the transformation properties ([524} and ( 1525b of the Green tensor given 
in App. IA.2l An experiment using 87 Rb atoms in their \F = 2, 771f = 2) hyperfine ground state 
has revealed spin flip lifetimes on the order of seconds for distances between 20 . . . 100/im [85]. 
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Fig. 14. Trapping lifetime r = 1/r as a function of the atom-surface distance z (right figure). Experimental 
data taken from Ref. [85], theoretical curve taken from Ref. [84]. The geometry is depicted on the left. The 
wire consisted of a Cu core (inner radius ai = 185 /im) and a 55 /im thick Al cladding. 



Figure [14] shows the experimental data together with the theoretical predictions according to 
Eq. ( |372t where the Green function for a three-layered cylindrical medium was used [84] (see 
also App. |A.4.2| >. 

It is instructive to investigate certain asymptotic regimes in which analytical approximations 
to the expected spin flip lifetime can be given. For this purpose, we consider an atom at a 
distance d away from a planar metallic surface with skin depth 5 and with thickness h. Its Green 
tensor can be found in App. IA.4.11 The skin depth 8 is related to the dielectric permittivity 
by e{u)) = 2ic 2 / (cu 2 S 2 ) . Together with the transition wavelength A, there are four different 
length scales involved, of which A is by far the longest and can be taken to be infinite. From the 
remaining three length scales one can find three experimentally relevant extreme cases that can 
be summarised in the following expression [86, 87]: 



r 

To 



1 



ana 



( 

35 ' 
S^d 

~2~' 

S 2 d 2 



2h 



S^d,h, 
S,h>d, 
, 5 >rf > h. 



(373) 



Here, To = 1 /To is the trapping lifetime in free space which, for a transition frequency u>a = 
2ir 400kHz, amounts to 3 • 10 25 s [88]. For fixed values of (d, h) there are two distinct functional 
dependencies on the skin depth S. These two scaling laws can be supported by rather intuitive 
explanations. If the skin depth is larger than the remaining length scales (i.e. the material is 
more dielectric), the magnetic field fluctuations weaken. As their strength is proportional to 
the imaginary part of the permittivity and hence the squared inverse skin depth, the lifetime 
increases quadratically with the skin depth 5 [second and third lines of Eq. d373H . If, on the 
other hand, the skin depth is the smallest length scale, the effective volume diminishes from 
which fluctuations emanate which leads to a lifetime increase with decreasing skin depth [first 
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Fig. 15. Left panel: Trapping lifetime as a function of skin depth for a thin substrate layer (h — lfim, dashed 
line) and a half-space (solid line). The atom-surface distance was chosen as d — 50/im, the transition 
frequency is to a = 2tt 560kHz [87]. Right panel: Spin flip lifetime r near a superconducting Nb slab 
(T c = 9.2K) as a function of atom-surface distance. Dashed line: two-fluid model; symbols: Eliashberg 
theory with various elastic scattering rates [89]. 



line in Eq. ( 13731 )1. This in turn means that there is a pronounced lifetime minimum when the skin 
depth is on the order of the atom-surface distance (Fig. [15]). To either side of this minimum, the 
expected lifetime increases drastically. However, choosing a surface material with larger skin 
depth limits its capability to generate strong enough magnetic traps. On the other hand, natural 
materials with skin depths below 50 /im at room temperature and 1 MHz are impossible to find. 

A possible solution is to make use of superconductors on the assumption that due to the 
supercurrent magnetic field fluctuations are shielded away from the vacuum-superconductor in- 
terface, and hence spin transitions are suppressed. For a bulk superconductor with conductivity 
<j(ijj) = ct'(lj) + ia"(u>) and cr"(uj) <C cr'(ui), we can rewrite the first line of Eq. ( 13731 ) as [89] 



to V 3 / n t h + 1 ^ c / 



(374) 



In the London two-fluid model [90] one assumes that the two types of charge carriers, normal 
and superconducting, react to an external field according to Ohm's law j„ = er rl E and the London 
relation = E, respectively. As a function of temperature, the fraction of normal conducting 
electrons follows the Gorter-Casimir expression n n (T)/no = (T/T c ) 4 [91]. With the plasma 
frequency ujp and the elastic scattering rate 7 of the electrons, the conductivity can be written in 
the form 



a(u>) = e up \ - \ 



1 \T< 



1 



(375) 



Although the two-fluid model does not fully capture the rich dynamics of superconductors (as 
it neglects coherence effects and dissipation), it gives a relatively accurate and intuitive picture 
of the strength of the magnetic field fluctuations. More elaborate models such as the Eliashberg 
theory prove that with regards to spin transitions the two-fluid model is perfectly adequate [89]. 
In Fig. [15] we show the distance-dependence of the spin flip lifetime for the two-fluid model of 
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superconducting Nb (dashed line) and the corresponding results for an Eliashberg calculation 
with varying elastic scattering rates (symbols). 
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5 Dispersion forces 

Dispersion forces such as Casimir forces between bodies [16,92], Casimir-Polder (CP) forces 
between atoms and bodies [93, 94] and van der Waals (vdW) forces between atoms [94, 95] 
are effective electromagnetic forces that arise as immediate consequences of correlated quantum 
ground-state fluctuations. The total Lorentz force on an arbitrary macroscopic or atomic charge 
distribution characterised by a charge density p and a current density j occupying a volume V is 
given by 



/ d\ fp(r)E(r)+j(r) x B(r) 

Jv 1 



(376) 



To see how this quantum force acquires a nonzero average even in the absence of external 
electromagnetic fields due to correlated zero-point fluctuations, let us consider the example of 
a neutral, stationary ground-state atom A. Expressing the atomic charge and current densities 
in terms of polarisation and magnetisation [recall Eqs. (l92l and d93ll. the Lorentz force can 
equivalently be represented as [57, 96, 97] 



3„ 



P A (r)-E(r)+M^(r)-B(r) 



(377) 



(Va: derivative with respect to the atomic centre-of-mass position). The quantum averages of 
both the electric and the magnetic field vanish in the absence of external fields, (E) = (B) = 0, 
and for an unpolarised and unmagnetised atom so do the atomic polarisation and magnetisation, 
(Pa) = (Mi) = 0. In the absence of correlations, this would imply that the average net force 
on the atom is vanishing, (F) = 0. However, both the electromagnetic and atomic fields are 
subject to nonvanishing zero-point fluctuations, (E 2 ), (B 2 ), (P\), (M- A ) ^ 0. These quantities 
are mutually correlated and thus lead to a nonvanishing dispersion force: For a nonmagnetic 
atom, one can show that this force is given by ( [98], cf. also Sec. 15.2. fl below) 



F = ^V A / duuj 2 lm\a(uj)trG( s \r A ,r A ,u;)} (378) 
4tti J 

■Va/ dw u? [lma(o;)tr Re G (s) (r A ,r A ,oj) + Re a(w)tr Im G (s) (r A , r A , u)) . 
■Jo 



47TI 



In accordance with the fluctuation-dissipation theorem, the real and imaginary parts of the atomic 
and field response functions (i.e. atomic polarisability and the Green tensor of the electromag- 
netic field) represent the reactive and fluctuating behaviours of these systems. Hence, the first 
term in the above equation correspond to the atomic zero-point fluctuations (Im a) giving rise to 
an induced electromagnetic field (Re G), while the second term is due to the fluctuations of the 
electromagnetic field (Im G) giving rise to an induced atomic polarisation (Re a). Both fluctua- 
tion sources thus contribute equally to the CP force. This interpretation is in contrast to the mode 
summation picture in which only the electromagnetic field fluctuations are taken into account, 
and the only role of the matter is to provide the perfect boundary conditions. 

It is worth noting that for small separations, dispersion forces are primarily due to the atomic 
zero-point fluctuations which interact via the instantaneous Coulomb interaction. Dispersion 
forces where first postulated by J. D. van der Waals [99] and theoretically analysed by F. London 
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[95] and J. E. Lennard-Jones [93] in this nonretarded limit. The zero-point fluctuations of the 
transverse electromagnetic field become important in the retarded limit of large separations, as 
was shown by H. B. G. Casimir and D. Polder [94]. To honour these two major steps, one often 
uses the notion vdW forces for all nonretarded dispersion forces on atoms and the term CP force 
for fully retarded ones — in contrast to the naming convention adopted throughout this work. 



5.1 Casimir forces 

The Casimir force (for general literature and reviews, see Refs. [23,97, 100-104]) on a body of 
permittivity e(r, u>) and permeability /i(r, uj) occupying a volume V can be found by calculating 
the ground-state expectation value of the Lorentz force ( 1376b 



F = J^r {({0}| [p(r)E(r') + j(r) x B(r')] |{0})} r( 



(379) 



(|{0}): ground-state of the body-assisted electromagnetic field) where the coincidence limit 
r' — > r must be performed in such a way that unphysical divergent self-force contributions 
are discarded after the vacuum expectation value has been evaluated. Using the relations 



p{r,u) = -£oV-{[£(r,a;)-l]E(r,a;)} + ^V-] N (r,w), 
j(r,w) = -iuJSo[e(T,uj) - l]E(r,w) 

+V x |/t [l - «(r,w)]B(r,w)} + j N (r,w), 



(380) 



(381) 



the field expansions (|225b , (|2261 l and the expectation values ( I242b -( |249l i, one finds for a homo- 
geneous body at zero temperature [105] 



h r f°° / ' uj 2 

F = - d 3 r du> — V-ImG (s) (r,r, 
n Jv Jo V c 



-- d 



tr-i I 



V x V x — - 



ImG (s) (r,r,w) X V'J 



/■oo / e2 

•J <%( 5_V-G( s )(r,r,iO 



- tr<^ I x 



V x V x +\ 

c 1 



G( s )(r,r,^)x V'jj, (382) 



where the coincidence limit has been performed by replacing the Green tensor with its scattering 
part. Here and in the following, the gradients V and V are understood to act on the first and 
second arguments of the Green tensor, respectively. 

Alternatively, the Casimir force can be equivalently expressed in terms of a surface integral 
rather than a volume integral. To that end, one makes use of the relation 



i 5(r)E(r)+j(r)xB(r) = V-f(r)-£ ^ E(r) x B(r] 



(383) 
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with the Maxwell stress tensor 

f (r) = e E(r) ® E(r) + Mo ^(r) ® B(r) - \ [e E 2 (r) + ^ X B 2 (r)l I . (384) 
This can be used to rewrite Eq. 13791 in the form 

F= / da- ({0}|{£ O E(r)®E(r')+^o lB ( r )® B ( r ') 

- | [ £Q E(r) ■ E(r') + ^(r) • B(r')] J}|{0}) (385) 

[note that the last term in Eq. ( 13831 l does not contribute in the stationary case]. Evaluating the 
field expectation values according to Eqs. (1242b — ( |249T >, one finds that (T = 0) [106] 

,2 



dui 


f 




av 


1 








2" 


c 2 



da ■ 



— ImG(r, r,w)-Vx ImG(r, r, uS) x V' 





1 

2 tr 



da ■ 



av 



£ 2 <- ' 

\ G(r, r, «C) + V x G(r, r, if) x V 



£ 2 <- ' 

\ G(r, r, if) + V X G(r, r, if) x V 



(386) 



The general expressions ( 1382b and ( 1386b can be used to calculate Casimir forces between 
bodies of arbitrary shape. In particular, Eq. (13861 ) in connection with the Green tensor for planar 
multilayered dielectrics (App. lAATT l immediately leads to the famous Lifshitz formula [92] for 
the Casimir force between two dielectric half-spaces [107]. In the limit of perfectly conducting 
plates, one recovers the mode summation formula ( 1731 ). 



5.2 Casimir-Polder forces 

Similarly to the Casimir force, the Casimir-Polder force on a single atom (for a general overview, 
see Refs. [23,26, 100, 108, 109]) can be calculated as an effective Lorentz force starting from 
expression ( 1376b , where the electromagnetic field acts on the atomic charge and current distribu- 
tions. In general, this force is a time-dependent quantity which can only be found by solving the 
coupled atom-field dynamics [57,96,97]. In this section, we restrict our attention to the force 
on a ground-state atom with the body-assisted field being in its vacuum state, and we ignore the 
effect of motion on the CP force. For this stationary problem, the CP force can alternatively be 
derived from the atom-field coupling energy, following the approach originally used by Casimir 
and Polder [94]. 

Such an approach can be justified by means of a Born-Oppenheimer approximation by as- 
suming that the fast internal (electronic) motion effectively decouples from the slow centre- 
of-mass motion. To see this, we express the total multipolar Hamiltonian (1292b in terms of the 
centre-of-mass momentum pa = X) a gAP" and the internal momenta p a = p Q — (m a /mA)pA 
to obtain 

H = + Hf + H F + H AF , (387) 
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where 



rrint 



p 2 

fa 

2m„ 



^- J d?r Pi(r) = ^ricj fc i fcfc 



(388) 



[Afefe = is the internal atomic Hamiltonian associated with the electronic motion and the 

field Hamiltonian H-p is given by Eq. ( 12951 ). The influence of atomic motion on the atom-field 
interaction can be discarded by formally letting tua — > oo, leading to 



H 



AF 



d 3 r P A {r) ■ E(r) 



d 3 r M A (r) ■ B(r) 



E — 



d 3 r@ a (r) x B(r) 



(389) 



where the three terms describe the electric, paramagnetic and diamagnetic interactions with the 
electromagnetic field, respectively. Having thus separated the internal and centre-of-mass motion 
as far as possible, we can apply the Born-Oppenheimer approximation by integrating out the 
internal motion for given values of f a and p^, leading to an effective Hamiltonian for the centre- 
of-mass motion, 



2rriA 



E + AE. 



(390) 



Here, E is the energy of the uncoupled atom-field system and AE is the energy shift due to the 
atom-field coupling Haf- Since we have neglected the influence of motion on the atom-field 
interaction, the energy shift does not depend on so that we may write 



AE = AE + AE(ta), 



(391) 



where AEq is the well-known (free-space) Lamb shift (1121b [23]. By means of the commutation 
relations \ta, Pa] = «W, the effective Hamiltonian thus generates the following equations of 
motion for the centre-of-mass coordinate: 



x 1 
m A YA = tt 

in 



m A rA,H eS 
1 



Pa, 



F = m A VA = tt m A rA,H eS 
in L 

where the CP potential 

U(r A ) - AE(r A ) 



-V a U(ta) 



(392) 
(393) 

(394) 



is the position-dependent part of the energy shift, as suggested by Casimir and Polder. In many 
cases, the centre-of-mass motion is effectively classical and position and velocity in Eqs. d3921 > 
and d393t can be treated as c-number parameters. In the following, it is not necessary to distin- 
guish the classical from the quantum case and we drop the operator hats for convenience. 
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(a) 



(b) 



/ 

/ 

/ 1 



Fig. 16. Schematic representation of single-photon interactions (a) and the second-order energy shift (b). 
Solid lines represent atomic states and dashed lines stand for photons. We do not distinguish electric and 
magnetic interactions. 



5.2.1 Perturbation theory 

For weak atom-field coupling, the ground-state energy shift ( 13911 ) can be obtained from a per- 
turbative calculation. While the coupling Hamiltonian d389t presents very general basis for this 
purpose, the calculation can often be simplified by applying the long-wavelength approxima- 
tion (which is valid provided that the atom-body separations are large with respect to the atomic 
radius) and neglecting the weak diamagnetic interaction, so that [cf. Eq. (I297H 

H AF = -d • E(r A ) - rh • B(r A ) (395) 

This Hamiltonian being linear in both atomic and field variables, AE is to leading order given 
by the second-order energy shift 

AE= y 2 (0\H AF \I){I\H AF \0) : (396) 

where \G) = |0) |{0}) denotes the (uncoupled) ground state of H 1 ^ + Hp and the relevant inter- 
mediate states \I) = \k)\l\(r, u>)) are those where the atom is excited, with a single-quantum 
excitation of the field being present, |lx(r, u)) = f|(r, w)|{0}). The formal sum in Eq. ( I396l l 
thus represents discrete summations over A and the vector index as well as integrations over r 
and lo. Recalling the field expansions ( 12251 and ( 1226b . the matrix elements of the electric and 
magnetic dipole interactions are found to be 

(0|({0}|d-E(rA)|lx(r,a;))|A;) - d ofc ■ G x (r A , r, u), (397) 

/ni/rmi - ■£./ n H / u, n m ofc ■ x G\(r A , r, w) 

(0 ({0} m-B(r A ) l A (r,u;)} \k) = ; ^ (398) 

The matrix elements and the second-order energy shift are schematically depicted in Fig. [T6] 
As the energy shift is quadratic in H A p, it contains purely electric, purely magnetic and mixed 
electric-magnetic contributions. The latter are proportional to dofc<g)mfco, so for nonchiral atoms, 
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they can be excluded by means of a parity argument (d is odd and rh is even under parity). After 
substitution of Eqs. ( 1397b and ( 13981 . we thus only retain the purely electric and magnetic terms, 
which upon using the integral relation ( 12301 > leads to 



AE ■■ 



Mo \ ■* 

7T ^ 



du 



^kO + ^ 



uj 2 d 0k ■ lmG(r A ,r A ,oj) • d 



/,(] 



m fe ■ V x Im G(r A ,r A ,v) x V' ■ m fe0 



(399) 



The Casimir-Polder potential d394t can be extracted from this energy shift by discarding the 
position-independent contribution associated with the bulk part of the Green tensor and 
only retaining its scattering part G^ s \ The result can be further simplified by writing Im G = 
(G — G*)/(2i), making use of the Schwarz reflection principle, Eq. ( 11931 ), and transforming 
the integrals along the real axis into ones along the purely imaginary axis (cf. Ref. [57]). The 
resulting ground-state CP potential is given by [57, 66, 96, 98, 100, 1 10-1 13] 



U(v A ) = U e (r A ) + U m (v A ), 



(400) 



with 



U e (r A ) 



U m {r A ) = 



2ir 



/"CO 

/ dU 2 tr[ a (^)-G^(r A ,r A ,iO] 
Jo 



f3(i0-VxGW(r A ,r A ,i0xV' 



d£tr 



2?r 



(401) 



(402) 



denoting the electric and magnetic parts of the potential and 

e— *u 

2 



2 m—^, ujkodok ® dfco 
lim — > 

e^0 h ' 



lim — 

3H (J 



^fco|doA;| 2 



2 

feO 



/3(w) 



lim - 



— VjJ€ 



I = a(u))I 



(403) 



J k0 



lim-V 



Wfco mo/fc 



^feO 



I = P(u)I 



(404) 



being the atomic polarisability and magnetisability. The second lines of equalities in Eqs. d40U - 
(14041 ) above hold for isotropic atoms. 

When considering CP forces on atoms that are embedded in a body or a medium, one has to 
account for the fact that the local electromagnetic field interacting with the atom differs from the 
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macroscopic one employed in our derivation. This difference gives rise to local-field corrections 
which can be implemented via the real-cavity model by assuming the atom to be surrounded by a 
small free-space cavity (recall the discussion in Sec. 14. II ). In order to apply this procedure to our 
results ( 140 11 1 and ( 14021 i. one has to replace the Green tensor by its local-field corrected counterpart 
as given in App. IA.6l which, after discarding position-independent terms, results in [58, 1 14] 



Ue (r A ) = ^ / 



U m (v A ) = ^±J™d€/3(tf) 



3e(»0 



+ 1 

3 



ttGW(T A ,T A) %Z), (405) 
tr [V x G( s \r A ,r A ,i€) x V'j , (406) 



where e(lu) — e(r A ,bj) and fi((j) = fi(r A ,u>) denote the permittivity and permeability of the 
host body at the position of the atom. For an atom situated in free space, the local-field correction 
factors are equal to unity and one recovers Eqs. ( 140 It and ( 14021 >. 

It is instructive to study the behaviour of the total CP potential under a duality transformation, 
which in this case amounts to a simultaneous global exchange a «-> (3/c 2 and e «-» ^ (see 
Tab.[2]i. Using the associated transformation laws (15241 i and (15251 l of the Green tensor as given in 
App. IA.2l one sees that the duality transformation results in an exchange U e <-> U m of the free- 
space potentials ( 140 U and (14021 i. so that the total potential is invariant with respect to a duality 
transformation [49]. The transformation laws (160 \\ and ( 1602b for the local-field corrected Green 
tensor imply that the same is true for an embedded atom, provided that local-field corrections are 
taken into account. The duality invariance of the CP potential is very useful when considering 
specific geometries: Once the electric CP potential of an atom in a particular magnetoelectric 
environment has been calculated, that of a magnetic one can be obtained by simply replacing 
a — > (3 /c 2 and exchanging e <-» fi. 



5.2.2 Atom in front of a plate 

Let us apply the general results to an isotropic atom which is placed above (z A > 0) a magne- 
toelectric plate of thickness d, permittivity e(uo) and permeability see Fig. [IT] Using the 
Green tensor of this very simple geometry as given App. IA.4.11 the electric CP potential (140 11 1 
reads [96,115] 



87r Jo 



-2k z za 







r s - 









(407) 



where r s and r p are the reflection coefficients of the half space for s- and p-polarised waves 
and k z = Im k z , with k being the wave vector of these waves in free space. By virtue of the 
duality invariance, we can obtain the magnetic potential from the above expression by replacing 
a — * (3 /c 2 and exchanging e /x, which is equivalent to an exchange of r s and r p [cf. Eq. ( I55U 
in App.EAB [58,116]: 



U m (z A ) 



8^2 



2 m) 



dn. 



-2k z za 



(408) 
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Fig. 17. An atom interacting with a magnetoelectric plate. 



Perfect mirror: A perfect mirror is realized for a perfectly conducting plate (e — > oo) or 
an infinitely permeable one (p, — > oo), in which cases the reflection coefficients are given by 
=pl, respectively. In this particularly simple case, the K 2 -integrals in Eqs. ( 1407b and 



1408b can be performed, so that the total CP potential ( 1400b reads 



U{z A ) = T 



l6ir 2 EQZ 



1 



(409) 



where the upper (lower) is valid for the perfectly conducting (infinitely permeable) mirror. This 
result includes the famous Casimir-Polder potential of a polarisable atom in front of a perfectly 
conducting wall [94]. In the retarded limit z A ^> c/bj m ^ (w m j n : minimum of all relevant atomic 
transition frequencies), the ^-integral is effectively limited to a region where the approximations 
a(z£)~ cha(0) and/?A(«0— /^(0) are valid and after integration one obtains [117] 



U(z A ) 



3hca(0) 



± 



3h/3(0) 



32ir 2 eoz\ 32-k 2 Eqcz\ 



(410) 



In the non-retarded limit za -C c/u> max (u) max : maximum of all relevant atomic transition fre- 
quencies), the factors ou(z£) and limits the ^-integral in Eq. (14091 ) to a range where we 
may approximately set e~ 2 ^ ZA / c ~ 1 and neglect the second and third terms in the square brackets. 
The integral can then be performed with the aid of the definitions ( 14031 ) and ( 14041 ) 



U(z A ) = T 



48irs z A 



± 



48ne c 2 zl 



(411) 



with the first term being the Lennard-Jones potential [93]. 

The results d4091 l-( 141 lb show that a perfectly conducting plate attracts polarisable atoms and 
repels magnetisable ones. For the nonretarded limit, this behaviour can be made plausible by 
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(a) 













& 







(b) 




Fig. 18. Image dipole construction for an (a) electric (b) magnetic dipole in front of a perfectly conducting 
plate. 



noting that the CP potential is due to the interaction of the atomic dipole moment with its image 
in the plate [93]. The image of an electric dipole moment d, located at a distance za away 
from a perfectly conducting plate, is constructed by a reflection at the £i/-plane, together with 
an interchange of positive and negative charges and is hence given by d' = (— d x , — d y , d z ), cf. 
Fig[T8la). The average interaction energy of the dipole and its image is hence given by [118] 



U z d' z ) 



TTl \ - 1 (d ' d ' 

Ue[ZA> 2 47re (2zA) 3 48^4 



(412) 



in agreement with Eq. (141 II) . where the factor 1/2 accounts for the fact that the second dipole is 
induced by the first one. 

On the contrary, a magnetic dipole rh behaves like a pseudo-vector under reflection, so its 
image is given by rh' = (rh x , m y , —m z ), cf. Fig[T8"Tb). The interaction energy of the magnetic 
dipole and its image reads 



1 (rh • rh' — 3m z m' z ) (rh 2 ) 
m{ZA > = 2 4^o (2z A ) 3 = 48tco4 



(413) 



again in agreement with Eq. (141 11) . The different signs of the CP potential associated with po- 
larisable/magnetisable atoms can thus be understood from the different reflection behaviour of 
electric and magnetic dipoles. 



Half space: Let us next consider a semi-infinite magnetoelectric half space of finite permittivity 
and permeability, which is a good model for plates whose thickness is large with respect to the 
atom-surface separation. Using the reflection coefficients d55U in App. IA.4.11 the electric and 
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magnetic potentials (|407T > and (1408b take the forms [58, 96, 1 15, 1 16, 1 19, 120] 



u e( z A) = / d£^a{i£) I dK z e 

e(i£)n z + Ki z 



^—2k z Za 



H{i£)K z 
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U m (za) — 
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-2k z za 



e(i£,)n z + k 1z 



(2 KzC - \ \ ^Mhi ~ Ku 

\ £ 2 J K^) K z + Klz 



with 



(414) 



(415) 



lira ku 



■[ e (»0/i(<0-l]+« 2 , 



(416) 



where ki is the wave vector inside the half space. In the retarded limit za ^S> c/uj m i n (with w m i n 
being the minimum of all relevant atom and medium resonance frequencies) the potentials take 
the asymptotic forms [121] 



U e {z A ) = 



dv 



U m {z A ) = - 



3ftca(0) 

1 n(0)v - ^n 2 {Q) - 1 + v 2 
v 4 fj,(0)v + ^n 2 (0) - 1 + v 2 
3/1/9(0) 



1 \ e(0)v - y/n^U) - 1 + v 2 



v 4 > e(0)v + yjn 2 (0) - 1 + v 2 



(417) 



64ir 2 eocz\ Ji 



dv 



1 \ n(0)v - y/n 2 (0) -1 + v 2 



v 2 v 4 J M (0)« + ^2(0) -1 + v 2 



1 e(0)v- a/« 2 (0)- 1 + v 2 
v 4 e(0)v + y/n 2 {0)~ 1 + v 2 



(418) 



[n(0) = v^Op)], while in the nonretarded limit 71(0)2,4 <C c/u> max (w max : maximum of all 
relevant atom and medium resonance frequencies), they are well approximated by 



t/e(*A) = 
(z^) = 



16ir 2 e z A J 



327r 2 ZA 

for a dielectric half space and by 
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Fig. 19. The potential of a polarisable ground-state two-level atom (transition frequency uja, dipole matrix 
element d) in front of a magnetoelectric half space is shown as a function of the distance between the 
atom and the half space for different values of /i(0) (u>p e /ujio = 0.75, uiTe/oJio = 1.03, UTm/^io = 1, 
7e/wio = 7mMo = 0.001) [115]. 



for a purely magnetic one. One can summarise these results by stating that the CP potential is at- 
tractive for two objects of the same (electric/magnetic) nature, e.g. a polarisable atom interacting 
with a dielectric half space, while being repulsive for two objects of different nature. 

In contrast to what is suggested by the findings for a perfect mirror, the different cases fur- 
thermore lead to different power laws in the nonretarded limit, with attractive potentials being 
proportional to l/z A and repulsive ones following a much weaker 1/za behaviour. When either 
the atom or the half space simultaneously exhibit electric and magnetic properties, both attrac- 
tive and repulsive force components are present, which may be combined to form a potential 
barrier. This is illustrated in Fig. [19] where we show the potential (14141) of a polarisable ground- 
state two-level atom near a magnetoelectric half space, whose permittivity and permeability have 
been modelled by the single-resonance forms 

e(w) = i + - — — . K") = 1 + - — • < 423 > 

It is seen that repulsive force components may lead to a potential barrier at intermediate distances 
while attractive forces always dominate close to the surface due to their stronger power law. For 
a polarisable atom, repulsive force components are associated with the magnetic properties of 
the half space, so the barrier increases in height as fi(0) increases. It follows from the retarded 
limit ( 1417b that the threshold for barrier formation is /j,(0) — 1 > 3.29[e(0) — 1] for a weakly 
magnetoelectric half space and /j,(0) > 5.11e(0) for a strongly magnetoelectric one. 
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Fig. 20. The potential of a polarisable ground-state two-level atom in front of a magnetoelectric plate is 
shown as a function of the atom-plate separation for different values of the plate thickness d [(J,(0) = 5; 
whereas all other parameters are the same as in Fig.ED [115]. 



Plate of finite thickness: Finally, we consider a plate of arbitrary thickness d. Use of the 
appropriate reflection coefficients (1552b in Eqs. (1407b and ( 1408b leads to the potentials [115] 

U e (z A ) = ^| HdUMiO r dK z e~ 2 ^ ZA 

[m 2 (*£) k z — K iz\ t&nh(ni z d) 



in n l° 2 A [£ 2 {iQn 2 z - K 2 lz }t&nh(n lz d) \ 

^ e > /2e(iOK»Ku + [e 2 (<0«2 + «i.]te^(«ixd)/' 



M*0 = ^r^ 2 « r *.e-^ 

8tH Jo c 2 7 4/c 

[e 2 (i£)k 2 — «?J tanh(Ki z d) 



2s(i£,)n z Ki z + [e 2 (i£)/-c 2 + K jJ tanh(Ki z d) 



£ 2 / 2n(i£)n z n lz + [ii, 2 (i£)n 2 z + n 2 lz ] tanh(«i z d) 

For a sufficiently thick plate, d ^> z^, one may approximate tanh(/tizcj) ~ 1 to recover the 
half-space results ( 1414b and ( 1415b . In the opposite limit of a thin plate, n(0)d <C 2U> the approx- 
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imation K\ z d <C 1 results in 

hfiod 



U e (z A ) 



U m (zji) = 



8tt 2 



<i££ a(i£) / dn z e 



-2k z za 



- 1 



hfiod 



2 m) 

c 2 

,2 



dn z e 



e 2 (^)^ 2 -«? 



i 



' 12 



(426) 



(427) 



As for the half space, these potentials reduce to simple power laws for large and small atom- 
surface separations. In the retarded limit, the thin-plate potentials read 



U e (z A ) = - 



U m {z A ) 



hca(0)d 
l60ir 2 e o z A 
HP(0)d 



14e 2 (0) - 

E (0) 
'lV(O)- 

. m(o) 



6^ 2 (0) - 1" 
6e 2 (0) - 1 



160ir 2 e cz A 
while for nonretarded distances, they become 

3hd 



U e (z A ) = - 

U m (z A ) = 



QAtt 2 EqZ a a 

H/jLod 



A JO 

oo 



s(0) 

e 2 m - 1 



3 m) [e(iO-l}[HiO 



for a purely electric plate and 

hfiod 



U e (z A ) = 

U m (z A ) = 



64n 2 z 2 A 
3hd 



64tt 2 eoz\ 



du 



[/i(iO-l][3M*0 



c 2 



(428) 
(429) 

(430) 
(431) 

(432) 

(433) 
(434) 



for a purely magnetic one. A comparison with the findings of Sec. l5.2.2l shows that when moving 
from a thick to a thin plate, all power laws are increased by one inverse power as a result of the 
linear dependence on d/z A . This similarity is due to the common microscopic origins of both 
forces (Sec. EH. 

The general behaviour of the CP potential for a thin plate being very similar to that for thick 
plates, one may expect that repulsive force components may lead to potential barriers regardless 
of the plate thickness. This is confirmed in Fig. [20j where we display the potential ( 14241 ) of a 
polarisable atom near plates of various thicknesses. We have chosen sufficiently strong magnetic 
properties to ensure the existence of a potential barrier. It is very low for a thin plate, reaches a 
maximal height for some intermediate thickness and then lowers slowly towards the asymptotic 
half space value as the thickness is further increased. 



Macroscopic QED — concepts and applications 



87 




Fig. 21. An atom interacting with a magneto-electric sphere. 



5.2.3 Atom in front of a sphere 

As a second application of the general theory, let us consider an atom placed above at distance 
ta from the centre of a magneto-electric sphere of radius R, permittivity e(lo) and permeability 
fi(co), see Fig. [21] After substitution of the Green tensor ( 1575b from App. IA.4.31 the electric 
potential ( 140 It in the presence of the sphere is given by [122] 



U e (r A ) 



^ J™ ^ 3 "«) f> + 1)| [h^(kr A )]\ l s + 1(1 + 1) 



r^ (1) (fcr A )- 


2 

r l + 


- 1 d[r A ^ (1) (fer A )] 


fcr A 




fcrA dr A 



,j, (435) 

where the sphere's reflection coefficients for s- and p-polarised spherical waves read 

r i = [ z Ji ( z )] 'ji ( z i ) - 31 jz) [ziji (zi)]' 

^)[zh < t 1 \z)]'j l (z 1 )-h < i 1 \z)[z 1 j l (z 1 )]' , 
r i = ffog) [zji (z)]' 'ji(zi) - ji(z) [ziji (zi)]' 

P " m[zhf\z)]'3M-hf\z)[zx3M]' 

[z — kR, k = i£/c, Z\ = kiR, k\ = n(i£)fc, n(oS) = ^/s(lu)^l(uj), fit (z): spherical Han- 
kel functions of the first kind, ji(z): spherical Bessel functions]. Using the duality invariance 
discussed in Sec. 15.2.11 the respective magnetic potential (1402b can be obtained by making the re- 
placements a — > (3/c 2 and e <-> fj,, which as in the case of the half space amounts to an exchange 



U m {r A ) 



' dtt 3 ^ f> + 1)1 [hf\kr A )]\ l p + 1(1 + 1) 
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1 2 



(438) 
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In the limit of a large sphere, R ^> za = ta — R, the main contribution to the sums in 
Eqs. ( 1435b and ( 14381 ) is due to terms with large n and one recovers the half-space results ( 14191 )- 
(14221) . as expected. For a sufficiently small sphere, n(Q)R -C r\, the sums effectively reduce to 
terms with 1 = 1, and one finds 



U e (r A ) 
U m {rA) 



327^6 jq 



— 5- 



mthg(x) = 2e~ 2x (3 + 6x 



-o' a Jo 
5x 2 + 



9{£,r A /c) - a Q (i£)h(£r A /c) 



(439) 
(440) 



2x 3 



and h{x) = 2x 2 e (1 + 2x + x 2 ), where we 



have introduced the small sphere's polarisability [118] 

a & (uj) = 4:-ks R 3 



,3 e M - 1 



e(w) + : 
as well as its magnetisability 

4tt.R 3 - 1 



[la + 2 



(441) 



(442) 



It is worth noting that in the small-sphere limit the electric and magnetic properties of the sphere 
completely separate. Furthermore, as will be seen in Sec. l5.3.2l below. the potentials (14391 ) and 
(14401 ) of an atom interacting with a small sphere have exactly the same form as that of two atoms, 
with the sphere's polarisability and magnetisability replacing those of the second atom. In the 
retarded limit ta S> c/o; m i n , the sphere potentials further reduce to 



23fea(O)a (O) 7fa(O)/3 (O) 

Ue\ r A) — ~ — 5 o i r 



U m (r A ) 



23ftc/3(O)/3 (O) 7ft/3(O)a (O) 
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while in the retarded limit ta -C c/w max one has 
3h f°° 

U e (r A ) = - 10 _3_ 2 „ 6 / dfa(i£)a©(*0 



'0 A JO 



U m (r A ) = - 



A 



167T 3 £Qr' 



167r 3 C 3 £nr 4 



0' A JO 



(445) 
(446) 



The hierarchy of signs and power laws of these potentials is closely analogous to that found for a 
half space or a thin plate: A polarisable atom is attracted to an electric sphere and repelled from a 
magnetic one, while the findings for a magnetisable atom are exactly opposite. Again, attractive 
and repulsive potentials follow the same (l/r A ) power laws in the retarded limit, while for short 
distances attractive potentials (with their l/r A power law) dominate over the repulsive (l/r A ) 
ones. 
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5.3 Van der Waals forces 

The simultaneous interaction of two atoms with the electromagnetic field leads to the vdW force 
between them (for general literature and reviews, see [22,97, 100, 123-127]), which may sensi- 
tively depend on their magneto-electric environment. In close analogy to the single-atom case, 
velocity-independent vdW forces between ground-state atoms can be derived from an atom-field 
coupling energy. Starting point is the two-atom generalisation 



H = -Bi- + -Bs_ + fff + £int + h f + H AF + H BF (447) 

2m a 2m B 



of the Hamiltonian ( 13871 i. where the internal dynamics of the atoms is given by Hamiltonians 
of the type ( 13881 l and each atom individually interacts with the electromagnetic field via an in- 
teraction Hamiltonian of the form ( 1389b . Applying the Born-Oppenheimer approximation by 
integrating the internal atomic dynamics for given centre-of-mass positions f a, f b and momenta 
PA, Pb, one obtains the effective Hamiltonian 

~ 2 ~ 2 

H cS = + + E + AE. (448) 

Zttia Ivtla 

Here, E denotes the energy of the two uncoupled atoms and the field, and the energy shift 

AE = AE + AE(r A ) + AE(r B ) + AE{r A ,r B ), (449) 

can be separated into a position-independent part (which contains the Lamb shifts of both atoms), 
two parts depending only on the positions of one of the atoms and a genuine two-atom part. The 
Hamiltonian ( 14481 ) generates the following equations of motion for atom A: 



■ 1 

m A r A = tt 
in 

Fa = m A r a = \ 
in 



PA, (450) 



m A rA,H t 



cir 



= -V A U(r A ) - V A U(r A , t b ) (451) 



(similarly for atom B). The atom is thus subject to two forces, the CP force, which as discussed 
in the previous Sec. l5.2l can be derived from the CP potential d394t , and the vdW force which is 
due to the additional atom. The associated vdW potential is given by the two-atom part of the 
energy shift 

U(r A ,r B ) =AE{r A ,r B ) (452) 

and it describes not only the direct, free-space interaction of the two atoms but also accounts 
for modifications of this interactions due to the presence of magnetoelectrics. As a conse- 
quence, Newton's third law F AB — — F BA does not necessarily hold for the vdW force Fab = 
— V A U(r A , fa), due to the contribution of the bodies to the momentum balance. 



5.3.1 Perturbation theory 



In close analogy to the single-atom case, the vdW potential can for sufficiently weak atom-field 
coupling be obtained from a perturbative calculation of the energy shift. With each atom being 
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0aJb)\U) 



Tab. 4. Intermediate states contributing to the two-atom vdW interaction, where we have used the short- 
hand notations = 11x^(^^,0}^)), = \h li i ll {*y,,U f t)l\ viv (T v ,W u )). 



linearly coupled to the electromagnetic field via an interaction ( 1395b . two-atom contributions 
start to appear in the fourth-order energy shift 



AE = £ 



(G\H AF + H BF \III)(III\H AF + H BF \II) 



1,11,111*0 ( Eg ~ ElIl) 



(II\H AF +H BF \I)(I\H AF + H BF \0) 
{Eq — Ejj){Eq — Ei) 



where |G) = |0a)|0b)|{0}) is now the (uncoupled) ground state of + + H F . In order 
to give a nonvanishing contribution to the energy shift, the intermediate states |J) and \III) 
must be such that one of the atoms and a single quantum of the fundamental fields are excited, 
while three possibilities exist for the intermediate states | II) : Either both atoms are in the ground 
state and two field quanta are excited [|1a(i", oj)1\> (r', lj')) = ^f|(r, cj)f|, (r', a/)|{0})], or 
both atoms are excited and the field is in its ground state or both atoms and two field quanta 
are excited. When invoking the additional requirement that each atom must undergo exactly 
two transitions, there is a total of ten possible combinations of intermediate states, which are 
listed in Tab. [4] The fourth-order energy shift ( 14531 ) thus involves both transitions from zero- 
to single-quantum excitations and those between single- and two-quantum excitations of the 
electromagnetic field. The former are given by Eqs. ( 1397b and (13981 ), the latter can, upon recalling 
the field expansions ( 1225b and ( 12261 ), be found to be 
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(k A \(l\ lll (r 1 ,uj 1 )\m A ■ B(r A )|7 A2l2 (r 2 , w 2 )^A 3 i 3 (r 3 ,w 3 )) |0a) 



£(13) 

v/2 



{m^°-VAxG A2 (r A ,r 2 ,W2)}. 5 (12) {m^V A xGj 3 (r A ,r3,w 3 )} 



1^2 



(12) 



with 



(455) 



(456) 



Comparison with Eqs. ( 1397b and ( 13981 ) reveals that the matrix element for the one- to two-photon 
transition is equivalent to the two possible combinations of one photon being created and one 
propagating freely. This is schematically depicted in Fig. [22 a). 

The various contributions to the vdW potential can be calculated by substituting the inter- 
mediate states from Tab. Eland the matrix elements ( 1397b . ( 13981 ). ( 14541 ) and ( 14551 ) into Eq. ( 14531 ). 
Let us begin with the intermediate-state combination (1) from Tab.|4] the calculation of which is 
schematically represented in Fig. l22T b). After expanding the product in Eq. ( 1453b , evaluating the 
delta functions ( 1456b and making use of the integral relation ( 1230b . one obtains 



k,l 



du> 



|w a w /2 dSf ■ImG(r A ,r B ,w) ■ dgdjf ■ ImG(r A ,r B ,a/) d ! 



- 0JU1 



df .JmG(T A ,T B ,Lj) x V'-mgJ [m™ ■ V x Im G(r B , r A , J) ■ A^] 
ojlo 1 [m° A k • V x lmG{r A ,r B ,uj) ■ dg] [dg • ImG(r B ,r A ,w') x V' • m^ fc 
m° A k ■ V x Im G(r A , r B , u) x V' • mgl [m™-Vx Im G(r B , r A , w') V ■ mf 



(457) 



Note that the single- to two-photon transition ( 1454b [or d455H with its two possible processes 
[Fig. l22T a)l enters the energy shift quadratic ally, so after expanding the product, four different 
terms arise, see Fig. l22T b). They can be grouped into pairs of equal terms, so one ends up 
with only two distinct contributions. After assuming real dipole matrix elements, these two only 
differ in their frequency denominators D{i a ) an d ^(16) which are given in Tab. [5] The energy 
shift (1457b contains electric-electric contributions (where both atoms undergo purely electric 
transitions), magnetic-magnetic contributions (where both atoms undergo magnetic transitions) 
and mixed electric-magnetic and magnetic-electric ones. As in the single-atom case, we have 
assumed nonchiral atoms and discarded all those contributions where at least one atom undergoes 
an electric and a magnetic transition. 

The contributions to the energy shift which are associated with the remaining intermediate- 
state combinations listed in Tab. |4]can be calculated in a similar way; all contributions are de- 
picted in Fig.|23](cf. also Ref. [22]). The contributions (4)-(6), (9) and (10) contain two-photon 
states, so they correspond to two distinct terms, similar to the contribution (1) studied above. 
However, for the contributions (4), (5), (9) and (10), one of those two terms separates into two 
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Fig. 22. Schematic representation of three-photon interactions (a) and the contribution (1) to the fourth- 
order energy shift (b). Solid lines represent atomic states and dashed lines stand for photons. We do not 
distinguish electric and magnetic interactions. 



single-atom processes [diagrams (4b), (5b), (9b) and (10b) of Fig. [23) and does hence not con- 
tribute to the vdW potential. For real dipole matrix elements, all genuine two-atom contributions 
only differ from Eq. (1457b by the frequency denominators Dr^ and possibly also by the sign of 
the terms in the third and fourth lines of the equation; signs and denominators are listed in Tab. [5] 
The total vdW potential U (ta,tb) — AE(ta, rg) = J^i can be obtained by summing 
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Tab. 5. Signs and frequency denominators associated with the intermediate-state combinations given in 
Tab.Hl 
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(458) 



which holds because the remaining parts f(tu, lu') of the integrands in Eq. ((457) are symmetric 
with respect to an interchange of lu and u>' . The ^/-integral can be performed with the help of 



1 



1 



LU + LU 1 LU — LU' 



7 ) Im G(r B ,r A ,Lo') 



--Lu[G{r Bl r A ,Lu) + G*(r B ,r Al iu)]. (459) 



Transforming the ^-integral to run along the positive imaginary axis, one finds [58, 100, 128-130] 

U(v A ,v B ) = U ee (r A ,r B ) + U em (r A ,r B ) + U me {r A ,r B ) + U mm (r A ,r B ), (460) 

U ee (r A ,r B ) =- M J°° da 4 tr[ a A(^) • G(r A ,r B ,i0 ■ a B (#) ■ G(r B ,r A ,iQ] 

= -^J Q da 4 a A (i0a B (i0tr[G(r A ,r B ,it;) ■ G(r B ,r A ,it)], (461) 
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Fig. 23. Schematic representation of all two-atom contributions to the fourth-order energy shift. 
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where the atomic polarisabilities and magnetisabilities are given by Eqs. (1403b and (1404b and 
the respective second lines of these equalities hold for isotropic atoms. In close similarity to the 
single-atom CP potential d400b — ( f402b . the two-atom vdW potential can hence be expressed in 
terms of the atomic response functions and the Green tensor of the electromagnetic field, where 
the latter connects the positions of the two atoms. 

In order to treat the interaction of atoms which are embedded in a magnetoelectric, we have 
to generalise the above results by taking into account local-field effects. In close analogy to 
the single-atom case, this can be achieved via the real-cavity model by substituting the respec- 
tive local-field corrected Green tensors (15961 l. ( 1598b — ( f600b from App. lA.6l into Eqs. ( 146 1 b — (1464b 
resulting in [58, 1 14] 
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Fig. 24. Interaction of two atoms embedded in a bulk magneto-dielectric medium. 



where e a (lu) — s{ v A, w) and fi A (uj) — I^{^a, (and similarly for atom E). 

The behaviour of the vdW potential under a duality transformation a <-> (3 /c 2 , e «-* [i follows 
immediately from the respective transformation laws of the Green tensor derived in Apps. IA.2I 
and |A.6l Using Eqs. d524t -( |527l >, one sees that the free-space potentials ( |461t -( |464t transform 
into one another according to U ee <-» U mm , U em «-* U me , so that the total vdW potential ( 1460b 
is duality invariant [49]. The same is true for embedded atoms when including local-field cor- 
rections, as Eqs. ( 1603b — ( f606b show. The duality invariance can be exploited when applying the 
general potentials to specific geometries; e.g., after calculation of U ee for a certain magnetoelec- 
tric body, U mm can be obtained by replacing a — > (3 /c 2 and exchanging e <-> /i. 



5.3.2 Two atoms inside a bulk medium 



Let us first consider the vdW potential of two isotropic atoms that are embedded in an infinite 
homogeneous bulk medium of permittivity e(oj) and permeability (i(oj) (Fig.l24l>. Using the bulk 
Green tensor d533K the local-field corrected potentials ( |465t and ( |466t take the forms [58, 1 14, 
131,132] 



UeeirA ' rB)= - 16ni r AB J [ 2£ 8 ( 1 4+ ) l]4 9inmr AB /c] 

U em (r A ,r B ) = TT J zri / d^ 2 a A {i^)/3 B (i^) 



(469) 



^r\ B J s s ^ s; [2e(i0 + 1] 2 [2 M (*£) + I] 2 

x h[n(it)Zr A B/c], (470) 



[t A b = \ V A — v b\\ n(u>) = y/e(u))/j,(u)), refractive index of the medium] with g (x) = e~ 2x (3 + 
6x + 5x 2 + 2x 3 + x A ) and h(x) = e~ 2x (l + 2x + x 2 ). Making use of duality invariance, we 
can apply the replacements a <-» (3/c 2 , e <-* u, to directly infer the remaining two potentials 
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[58,114,131,132] 



x h[n(i^r AB /c], (471) 

(472) 



These results and their retarded and nonretarded limits given below generalise the well-known 
free-space potentials [117, 133-136]; in particular, U ee reduces to the famous Casimir-Polder 
potential of two polarisable ground state atoms in free space for e = fi = 1 [94]. In free space, 
the potential between two electric or two magnetic atoms is attractive and that between an electric 
and a magnetic one is repulsive, in agreement with the general heuristic rule that dispersion forces 
between objects of the same electric/magnetic nature are attractive and those between objects of 
opposite nature are repulsive. 

An inspection of Eqs. ( I469l )-( l472l i reveals that a bulk magneto-electric medium can influence 
the strengths of the various two-atom interactions, but cannot change their signs. In order to 
discuss the effect of the medium in more detail, it is useful to consider the limits of large and 
small interatomic separations. In the retarded limit r AB 3> c/w m i n , the potentials are well 
approximated by 



U ee (r A ,r B ) 
U em (r A ,T B ) 
U me (r A ,r B ) 
U mm (r A ,r B ) 



23hca A (0)a B (0) 81e 2 (0) 



647T 3 £n^ 
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647r 3 £ ^ 



' AB n(0)[2 £ (0) + l] 2 [2 M (0) + ip 
7hcno/3 A (0)a B {0) 81^(0)e(0) 



AB 



n(0)[2 M (0) + l] 2 [2e(0) + l] ;; 



647T 3 £o? 



23hc^ (3 A (0)p B {0) 81/i 2 (0) 



64ttV 



AB 



n(0)[2/x(0) + I] 4 



(473) 
(474) 
(475) 
(476) 



In this case, the influence of the medium on all four types of potentials is very similar: The 
coupling of each atom to the field is screened by a factor 3e(0)/[2e(0) + l] 2 for polarisable 
atoms, and a factor 3/i(0)/[2/j(0) + l] 2 for magnetisable atoms; in addition, the reduced speed 
of light in the medium leads to a further reduction of the potential by a factor n(0). In the 
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Fig. 25. Interaction of two atoms in the presence of a magnetoelectric half space. 



nonretarded limit tab <C c/u mw , the medium-assisted potentials simplify to 



U ee (r A ,r B ) = 
U em (r A ,r B ) = 
U me (r A ,r B ) = 
U mm (r A ,r B ) = 



16^ £ 2 r% J 



167T 3 '/ 



AB JO 



~ 16tt 3 / 6 



8l£ 2 (^) 

[2e(i0 + If ' 

Mig) + i] 2 [2 M (iO + 1] 2 ' 

8l e 2 (iQ/x 2 (iQ 
[2 £ (iC) + l] 2 [2 M (iO + l] 2 ' 



deC 2 /3A(^)aB«) 



(477) 
(478) 
(479) 
(480) 



The potentials [7 ee and U mm are thus again reduced, where in the nonretarded limit, U ee is only 
influenced by the electric properties of the medium and U mm only by the magnetic ones. On 
the contrary, the mixed potentials U em and U me are enhanced by a factor of up to 81/16 in the 
nonretarded limit. 



5.3.3 Two atoms near a half space 



A modification of the vdW interaction does not only occur for atoms embedded in a medium, 
but can also be induced by a magnetoelectric body placed near to atoms in free space. To see 
this, let us consider the body-assisted vdW interaction of two atoms placed above (z A , z B > 0) a 
semi-infinite half space of permittivity e(uj) and permeability (i(oj), see Fig. [25] For simplicity, 
we assume both atoms to be isotropic and nonmagnetic, so that U = U ee as given by Eq. (146 U . 
Separating the Green tensor into its bulk (free-space) and and scattering parts according to G = 
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G?(°) + G^ s \ the potential reads [130] 

ha 2 r°° 

U(r A ,r B )=-^ / dU 4 a A ^)a B (i0tv[G^(r A ,r B ,z0 ■ G^(r B ,r A ,zO] 
Z7r Jo 

-Mr dZ?a A (iQa B mtr[GW(T A , r B ,i£) • G (S) (r B , v A , i$ 
27T Jo 

+ G {s \r A ,r B ,iO ■ G {0 \r B ,r A ,iO] 

/ de^(ie)a B (»e)tr[G( s )(r A ,r B) tO-G (5) (r B) r A ,iO]. 
z!7r Jo 

(481) 

The first term is the free-space potential Uq{t a , t b ) which is due to a direct exchange of two 
photons between the two atoms according to one of the various processes depicted in Fig. [23] 
The second and third terms represent the body-induced modification of the potential; they are 
due to processes where one or both of the two exchanged photons are scattered of the surface of 
the half space before being absorbed. The free-space part of the potential is just a special case of 
the bulk-medium potential ( 14691 ) calculated in Sec. 15.3.21 while the body-induced change can be 
found by employing the free-space Green tensor ( 15361 ) from App. lA.3l together with the scattering 
tensor of the half space from App. IA.4.11 As the resulting expressions are rather involved due 
to the large number of geometric parameters, we only give analytical formulae for some simple 
special cases. 



Perfect mirror: We first consider a perfect mirror, with the reflection coefficients being given 
by r s = —r p = — 1 for a perfectly conducting plate and by r s = —r p = 1 for an infinitely 
permeable one. In the retarded limit z A , z B , r AB ^> c/uj m i n , the vdW potential in the special 
case x AB <C z A + z B — z + can be given in closed form [130, 137, 138] 



U(r A ,r B ) 



23hca A (0)a B {0) 



64tt 3 £ ! 



1 



AB 



32 



L AB 



6r AB 



23 r\ B z + (r AB + z + f 



1 



(482) 



(r AB = \r A — r B |, x AB = \x A — x B \) where the first term in square brackets is the free-space 
potential due to the direct exchange of two photons and the second and third terms represent the 
mirror-induced modification of the potential due to scattering of one or two photons of the mirror 
surface (where the different signs refer to the two cases of a conducting or permeable mirror). 
An infinitely permeable plate thus always leads to an enhancement of the retarded interaction 
of two polarisable atoms whereas two cases need to be distinguished for a perfectly conducting 
one: When the atoms are aligned parallel to the plate (z AB = \z A — z B \ = 0), Eq. (14821 ) shows 
that the potential is always reduced due to the presence of the plate, while for a perpendicular 
alignment (x AB — 0), the potential of reduced if z B /z A < 4.90 (atom A being closer to the 
plate than atom B). 

In the nonretarded limit z A , z B ,r AB <§; c/w max , the vdW potential in the presence of a 
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perfect mirror is found to be [130, 137, 138] 

U(r A ,r B ) = - 32 / c^aA(^) a s(^) 
io 71 " Jo 



1 _ — 2z AB z? + x AB (z 2 _ + z 2 AB ) 1 



V AB AB + ~ ^ A B \ * - ■ ■ './_;' 



(483) 



(r + = y + z\_ ). For a parallel alignment of the two atoms, a perfectly conducting plate 
thus reduces the vdW potential of two polarisable atoms while a permeable plate leads to an 
enhancement. In particular in the on-surface limit z + — > 0, the reduction and enhancement 
factors with respect to the free-space potential are given by 2/3 and 10/3. On the contrary, for a 
vertical alignment the potential is reduced by a conducting plate while for a permeable plate it is 
reduced if z B j za < 14.82 (atom A being closer to the plate than atom B). 

The enhancement and reduction of the nonretarded vdW interaction due to a perfect mirror 
can be understood from the interaction of the fluctuating atomic dipole moments A and B and 
their images A' and B' in the plate, with 

i^dipoic = Vab + V AB > + Vba> (484) 

being the corresponding interaction Hamiltonian. Here, Vab denotes the direct interaction be- 
tween dipole A and dipole B, while Vab 1 and V B A' denote the indirect interaction between 
each dipole and the image induced by the other one in the plate. The leading contribution to the 
two-atom energy shift is of second order in i?dipoic, 



, r , , (0 a I (0B\H int \l B )\k A ) (k A | (l B \H int \0 B ) \0 A ) ... 
AE(r A ,r B ) = - ^ n(ljk o 7 . ,m ■ (485 > 



The three terms in Eq. ( 1483b can be identified as different contributions to this energy shift: The 
first term is due to contributions that are quadratic in the direct interaction Vab and is always 
attractive due to the global minus sign in Eq. ((485}. The third term which is associated with 
quadratic contributions of indirect dipole-image interactions Vab 1 or Vba' i s attractive for the 
same reason and thus always acts as an enhancement of the free-space potential. The second term 
in Eq. ( 14831 ) is due to mixed direct-indirect interactions, its sign is negative for attractive dipole- 
image interactions (thus tending to enhance the free-space potential) and positive for repulsive 
dipole-image interactions (leading to a reduction). The overall effect of the perfect mirror thus 
depends on the attraction or repulsion of the dipole-image interaction and, in the case of repul- 
sion, on the relative strength of the mixed contributions compared to the purely indirect one. 

We begin our investigations with the case of two (polarisable) atoms placed aligned parallel to 
a perfectly conducting plane. The image-dipole construction [Fig. l26t a)l reveals that the indirect 
interactions Vab 1 and Vba' are repulsive, explaining why the second term in Eq. ( I483l l acts 
to reduce the interaction. Equation ( 14831 ) shows that it is always dominant over the enhancing 
purely indirect interaction, leading to an overall reduction of the free-space potential. The image- 
dipole construction further reveals that the quantity r + entering the mirror-induced modification 
of the potential is simply the distance between an atom and the image of the other atom. Let us 
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(a) (b) 




Fig. 26. Image dipole construction for two (a,b) electric and (c,d) magnetic dipoles in front of a perfectly 
conducting plate. 

consider next the case of perpendicular alignment of the two atoms. As seen from the image- 
dipole construction [Fig. l2oT b)l. the indirect interactions are attractive in this case, so that the 
mirror enhances the potential, in agreement with Eq. (1483b . 

The case of two polarisable atoms near an infinitely permeable plate can be addressed by 
studying the interaction of two magnetisable atoms near a perfectly conducting plate instead, 
since the two problems are equivalent by virtue of duality. Noting that magnetic dipoles behave as 
pseudovectors under spatial reflection, one finds that for parallel alignment of the two atoms the 
dipole-image interaction is attractive [Fig.|26lc)l like the dipole-dipole interaction, so all mirror- 
induced terms enhance the free-space potential as predicted by Eq. (148 31 l. For perpendicular 
alignment the indirect dipole-image interaction is seen to be repulsive [Fig. l26*l' d)l. so the mixed 
direct-indirect interaction tends to reduce the interaction [cf. the second term in Eq. ( I483H 
while the purely indirect one acts towards an enhancement. For small atom-atom separations, 
the former clearly dominates due to the enhanced strength of the direct interaction, leading to a 
reduction of the overall potential while for large atom-atom separations, an enhancement due to 
the influence of the purely indirect interaction may be expected. The two cases are obviously 
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Fig. 27. The vdW potential for two identical two-level atoms (transition frequency lja) aligned parallel 
(zab = 0) to the surface of an (a) purely dielectric half space (up e /uA = 3, ujtb/^a = 1, Je/wA = 
0.001) (b) purely magnetic half space {up m /uA = 3, ujTm/uA = 1, 'Ym/wA = 0.001) is shown as 
a function of the interatomic separation tab- The potential is normalised with respect to the free-space 
potential Uo{tab)- The atom-half-space separations are za = zb = 0.01c/wa (solid line), 0.2c/o;a 
(dashed line), and c/uja (dotted line) [130]. 



separated by the inequality given below Eq. ( 14831 ). 

Magnetodielectric half space: For a semi-infinite half space of finite permittivity e(uj) and 
permeability fi(ui), the vdW potential ( 148 It is found by substitution of the free-space Green 
tensor ( 1536b and the scattering Green tensor from App. |A.4T| together with the reflection coeffi- 
cients ( 15511 ) of the half space. The arising integrals over frequency and wave-vector components 
have to be solved numerically in order to obtain the half-space assisted potential for arbitrary 
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Fig. 28. The vdW potential for two two-level atoms aligned perpendicular (xab = 0) to an (a) purely 
dielectric half space and (b) purely magnetic half space is shown as a function of the interatomic separation 
tab- The distance between atom A (which is closer to the surface of the half space than atom B) and 
the surface is za = 0.01c/wio (solid line), 0.2c/ljio (dashed line), and c/ujio (dotted line). All other 
parameters are the same as in Fig.|27||T301. 

interatomic and atom-half space distances. The results are shown in Figs. l27l and |281 for two 
identical two-level atoms with the permittivity and permeability of the half space being given 
by the single-resonance models d423K where we display the relative modification of the vdW 
potential with respect to its free-space value Uq as found from Eq. ( |469l l. 

Figure [27] shows the case of the two atoms being aligned parallel to the surface of the half 
space (zab — 0). In agreement with the findings for a perfect mirror, Eqs. ( 14821 ) and ( 14831 ), 
the potential is always reduced by a purely dielectric half space [Fig. 1271^ )1 and enhanced by 
a purely magnetic one [Fig.lZTl'b)!. The relative reduction/enhancement becomes noticeable as 
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Fig. 29. Two atoms interacting with a magneto-electric sphere. 



soon as the interatomic separation becomes comparable to the atom-surface distance, it saturates 
for large interatomic separations. The figure further shows that the relative reduction for a di- 
electric half space has a pronounced minimum, the enhancement due to a magnetic one increases 
monotonically with interatomic distance. 

In Fig. [28] we study the case of perpendicular alignment (xab — 0). In agreement with 
the perfect-mirror results (14821 ) and ( 14831 ). a purely dielectric half space is found to always en- 
hance the potential. The enhancement sets in when interatomic and atom-surface distances be- 
come equal, reaches a maximum in some cases (i.e., whenever the atom-surface separation of 
the closer atom A is sufficiently small) and saturates for large interatomic distances. At first 
glance, the findings for a purely magnetic half-space seem to indicate a global enhancement 
which monotonically increases with the interatomic separation, thus being in contradiction with 
the perfect-mirror result (14831 ). However, a closer look [cf. inset in Fig. [28lb)l reveals that for 
very small interatomic and atom-surface separations a reduction can indeed be found, as pre- 
dicted by Eq. d483l 



In order to demonstrate the effect of different geometries on the vdW potential, consider next to 
polarisable atoms placed at distances ta and tb from the centre of a magneto-electric sphere of 
radius R, permittivity e{u>) and permeability fj,(us), with the separation angle of the two atoms 
being denoted by Bab, see Fig. [29] The sphere-assisted vdW potential is given by Eq. ( 1481b 
where the scattering Green tensor G^- ' for the sphere can be found in App. IA.431 For a small 
sphere (R <C ta, t b) in the nonretarded limit (ta, tb,tab <C c/uj max ), one can show that the 
sphere-induced part Uq = U — Uq of the potential reduces to [139] 



[&A = ?aI t a, e B = ?bI t b, e AB = { r A — r B)/ r AB] where a Q (uj) is the polarisability of the 
sphere as given by Eq. ( 1441b . When replacing the sphere by a third atom C (a — + cue), our 
results coincides with the Axilrod-Teller potential of three atoms [140]. The sign of the three- 
body potential ( 14861 ) depends on the geometric arrangement of the three objects. For instance, an 
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Fig. 30. The vdW potential for two identical two-level atoms placed at equal distance (ta =Tb) from the 
centre of an (a) purely dielectric sphere (b) purely magnetic sphere is shown as a function of the interatomic 
angle 9ab- The potential is normalised with respect to the free-space potential Uo{6ab)- The sphere 
radius is R = c/lua and the distances of the atoms from the centre of the sphere are ta — rs = 1.03c/uja 
(solid line), 1.3c/wa (dashed line), and 2c/ua (dotted line). The medium parameters are the same as in 
Fig.[27][139]. 



attractive potential is found when they are placed in a straight line while the potential is repulsive 
when they form an equilateral triangle. 

In Fig. [30] we show the total vdW potential of two identical two-level atoms placed at equal 
distance from a purely dielectric or purely magnetic sphere whose magnetoelectric response 
is given by Eqs. ( 14231 >. For small separation angles and atom-sphere separations, the sphere- 
induced modification of the vdW potential is very similar to our findings for two atoms placed in 
parallel alignment with respect to a half space: A dielectric sphere reduces the potential, while 
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a magnetic one leads to an enhancement. However, the figure also shows that the behaviour is 
reversed as the angular separation grows. For an electric sphere, this enhancement of the potential 
for large angular separations can be understood from the fact that the Axilrod-Teller potential is 
attractive when the atoms and the sphere are situated on a straight line. 



5.4 Relation between dispersion forces 

The three types of dispersion forces considered in the previous three sections have a common 
origin and are thus closely related. To see this, let us start from the Casimir force ( I382l i on a 
dielectric body of volume Vi which is situated in free space and interacts with a second dielectric 
body of volume V%. Assuming the first body to consist of a dilute gas of atoms [number density 
r\\ (r), polarisability a(oj)] such that the linearised Clausius-Mosotti law [118] 

r)(r)a(oj) 

e(r,u) - 1 = /W V ' (487) 



holds, one can make use of the linear Born expansion given in App. lA.5l to show that [105] 

d 3 rj)(r)VC/(r) (488) 



'Vi 

where U(r) is the Casimir-Polder potential ( I40U . To leading order in the atomic polarisability, 
the Casimir force on the body is simply the sum of the CP forces on the atoms contained inside 
it. One can repeat the exercise for the second body to find 



U{v A )= / d 6 r V {r)U{r A ,r), (489) 
Jv 2 

so the CP potential U(r) of each atom in body 1 with body 2 is due to its vdW interactions ( 146 It 
with the atoms in body 2. Combining these results, one has [105] 



F = - [ d 3 rr](r) [ dVr((r')VC/(r,r'). 



(490) 

'Vi Jv 2 

Hence, to leading order the Casimir force between the two bodies is a sum over all possible vdW 
forces between atoms in body 1 and atoms in body 2. For weakly dielectric bodies, both Casimir 
and CP forces may thus be regarded as a consequence of two-atom vdW forces. These results 
naturally generalise to weakly magnetic bodies and magnetic atoms. 

For bodies with a stronger magnetoelectric response, this simple additivity breaks down due 
to the influence of many-atom interactions [141-144]. For instance, using the exact Clausius- 
Mosotti law [118] 

e(r,w)-l = - (491) 

1 - r?(r)a(w)/(3e ) 



together with the full Born expansion, one can show that the CP potential ( 14891 ) generalises 
to [143,144] 



U(r A ) = J2 /"'/ ^i^i) ' ' ' d3r K ti(t k )U(t a ,tu . . . ,r K ) 



(492) 
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Tab. 6. Asymptotic power laws for the forces between (a) two atoms, (b) an atom and a small sphere, (c) 
an atom and a thin ring, (d) an atom an a thin plate, (e) an atom and a half space and (f) for the force per 
unit area between two half spaces. In the table heading, e stands for an electric object and m for a magnetic 
one. The signs + and — denote repulsive and attractive forces, respectively. 



where U{r\, . . . , r/v) denotes iV-atom vdW potentials [145, 146]. The full CP potential of a 
single atom in the presence of a body is just due to the whole hierarchy of its 2,3,... iV-atom 
interactions with the body atoms. 

To illustrate the close relation between dispersion forces, let us look at a few simple examples, 
most of which have already been studied throughout this Sec. [5] In Tab. |6l we list signs and 
leading power laws of the dispersion forces between two atoms (Sec. 15.3.21 ). an atom interacting 
with a small sphere (Sec. I5.2.31 l. a thin ring [144], a thin plate (Sec. I5.2."2b and a semi-infinite 
half space (Sec. l5.2.2l i. and that between two half spaces [147, 148]. For weakly magnetoelectric 
objects, the rows (b)-(f) of the table follow by pairwise summation of the vdW forces displayed 
in row (a): Summation over the compact volumes of a small sphere (b) or a thin ring (c) does not 
change the respective power laws, while summation over an infinite volume lowers the leading 
inverse power according to the number of infinite dimensions. So, the leading inverse powers 
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are lowered by two and three for the interaction of an atom with a thin plate of infinite lateral 
extension (d) and a half space (e), respectively. The power laws for the force between two half 
spaces (f) can then be obtained from the atom-half-space force (e) by integrating over the three 
infinite dimensions where integration over z lowers the leading inverse powers by one while the 
trivial integrations over x and y yield an infinite force, i.e., a finite force per unit area. Note that 
all of these power laws remain valid for objects with stronger magnetoelectric bodies, so we can 
conclude that many-atom interactions do not change the asymptotic power laws. 

The common features of all dispersion forces listed in the table are as follows: Forces be- 
tween objects of the same (electric/magnetic) nature are attractive while those between objects 
of opposite nature are repulsive. The combinations e ^> e and m <-> m obviously lead to the 
same behaviour of the force as an immediate consequence of duality invariance, the same hold 
for the combinations e <-> m and m <-> e. Attractive and repulsive forces generally follow the 
same power law in the retarded limit, while in the nonretarded limit attractive forces are stronger 
by two inverse powers in the object separation. 



5.5 Thermal effects 



The Casimir-Polder force acting on a stationary atom at finite temperatures can be derived from 
the Lorentz force ( 13771 which for a stationary nonmagnetic atom in electric dipole approximation 
takes the simple form 



F(r A ,t) 



Vd • E(r) 



(493) 



This expression can be evaluated using the statistical averages of the relevant electromagnetic 
field operators given in Eqs. 124ZI — A2451 . We use the multipolar-coupling Hamiltonian ( |294t in 
electric-dipole approximation to solve the coupled dynamics of the electromagnetic field and the 
atom which, upon using the solutions to the equations of motion for the internal atomic dynamics 
(Sec. reads as [149] 



F(v A ,t) = / dwuj 2 Vd nr ImG( s '(r,r A , W ).d s 



dr(A nn (T)) 



{n th He [l( " + ^ fcMr " +rfc)/2l( '- T) + [nthM + 1] e^ u ~ u »»)-(r»+r»)/2](t-r) J 
c.c. (494) 



After evaluating the t integral in Markov approximation and the lu integral using countour- 
integral techniques, the thermal Casimir-Polder force on an atom in an incoherent superposition 
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of energy eigenstates is given by Ffr^, t) = J2 n a nn(t)^n[ r A) with force components 



F„(r A ) = -^T ]T (l - ±6 m )& 

x Vtr{[a n (i£ N ) + ot n (-i£ N )} ■ G {s \r A ,r,i£ N )}\ r=rA 

+ i"o y^{Q(^nfc)^nfc[nth(Onfc) + l]Vd„ fe • G (s) (r, r A , O nfe ) • d fc „| r=r 

- 6(a; fcn )n^ nth (o* n )Vd nfe • G< s )(r, v A ,n* kn ) ■ d kn \ r=rA + c.c.} 
[O n fe = o) rl fc + i(r„ + Tk)/2] and atomic polarisability 

1 Y^r d nfe ® dfe„ , dfc n ®d n fe~ 



-a 



-a* 



(495) 



(496) 



One important feature of the result (I495l l is the appearance of resonant force contributions 
proportional to nth(w ra fc) and nth(w„/c) + 1 which are due to absorption and emission of photons, 
respectively. Even for ground-state atoms there exists a resonant force contribution 



F cs (r A ) = -i^fco«thK-o)|d fe| 2 V A trReG (5) (r A ,r A ,w fco ) 



(497) 



associated with absorption of thermal photons at the atomic transition frequency tvko- This result 
is in contrast to the frequently used Lifshitz result [150] 



TiLifshitz 



oo 



(t a ) = -MbT > (1 - ^S m )^ N a(^ N )V A trG^ s >{r A , r A , i£ N ) 



N=0 



(498) 



which only contains the non-resonant force contributions. Clearly, because of the thermalisation 
of the atom, these resonant forces can only be observed on time scales that are short compared 
to the inverse ground-state heating rates T^ 1 . As their magnitude scales with the mean thermal 
number n t h(<jJko), the atomic transition frequency must not be too large in comparison with the 
ambient temperature, hiouo S ^bT. Ideal candidates to observe these resonant force contribu- 
tions would therefore be polar molecules whose vibrational and rotational frequencies are small 
enough to yield large thermal photon numbers at room temperature, and whose heating times can 
reach several seconds (see also Sec. [43} [78]. 

The atom thermalises after times much longer than T^ 1 and reaches its stationary (thermal) 
state 



exp 



J2 n M> n A nn /(k B T) 



tr exp 



J2 n h£j n A nn /(k B T) 



5u> r , 



(499) 



in the long-time limit. In this limit all resonant force contributions cancel, and the Casimir- 
Polder force can be written in the form ( 14981 ) only if the atomic polarisability has been identified 
with the thermal polarisability 



VT,nnOt n {u) ■ 



(500) 
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Without this identification, the equilibrium force ( 1498b is generally larger than the one obtained 
from Eq. ( 14951 . The microscopic degrees of freedom of an atomic system thus give rise to a 
rich additional structure of CP forces, which is not explicitly present in the macroscopic Casimir 
force. 
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6 Cavity QED effects 

In this final section, we turn our attention to the application of macroscopic QED to strong atom- 
field coupling scenarios. This is of particular interest in cases when the atom or molecule is 
confined in resonator-like structures in which the photon round-trip time is short compared to the 
atomic relaxation time. 



6.1 Quantum-state extraction from a high-Q cavity 

The interaction of a two-level atom with light in a resonant cavity is of fundamental importance 
to our understanding of quantum optics, and provides a playground to test ideas of quantum infor- 
mation processing such as the generation of entangled quantum states of two or more atoms [151] 
or the realization of single-photon sources. In a previous section (Sec. 12.2.41 ) we described the 
idealised situation of perfectly coherent atom-light interaction that led to the Jaynes-Cummings 
model. Within this model it is possible to completely transfer the excitation of a two-model atom 
to the cavity field and back coherently. The idea to use a high-Q cavity to deterministically gen- 
erate single photons in well-defined quantum states is at the heart of cavity QED [152, 153]. Here 
we will examine the closely related problem of how to extract a photon from a cavity within the 
framework of macroscopic QED which provides us with the means to treat leaky cavities in a 
quantum-theoretically consistent way [51]. 

For simplicity, we will focus on one-dimensional cavities such as in Fig. [5T| where region 
denotes a perfect mirror, region 2 the semi-transparent mirror, region 1 the cavity interspace, and 
region 3 the free space surrounding the system. In order to calculate the response of the atom- 
cavity system, we assume that the atom is initially being prepared in its excited state |e) and the 
cavity field is in its ground state |{0}). Then, the state vector at a later time can be expanded into 

oo 

\m) = O»(t)|e)|{0}) + JdzJ dujC g (z,uj,t)e- iut \g)\l(z,u;)) . (501) 

o 

The Schrodinger equation with the electric-dipole Hamiltonian (12971 ) yields coupled equations 




3=0 1 2 3 



Fig. 3 1 . The semi-transparent mirror of the cavity (region 2) is modeled by a dielectric plate, and the atom 
inside the cavity (region 1) can be embedded in some dielectric medium. Figure taken from Ref. [51]. 
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of motion for the unknown coefficients C e (t) and C g (z, u>, t) as 



C e (t) = - 



— , diu— dz y/lme(z, u)G(z A , z, uj)C g (z, u, t) e -<(«-«A)t > 

o 

(502) 
(503) 



- y/Jme(z, v)G* (z A , z, u)C e {t)e^- UA)t 



Cg(z,U),t) = , 

y/7Tn£ A (r 

[A: cross-sectional area of the cavity]. Inserting the formal solution of Eq. (15031 ) into Eq. ( 1502b 
yields an integro-differential equation 



C e (t) = / dt'K(t-t')C e {t') 



with the integral kernel 

oo 



(504) 



(505) 



The spectral response is carried by the Green function G(z, z' , u>) which has to be determined 
by adjusting the boundary conditions at the interfaces between the regions of piecewise constant 
permittivity. Because of multiple scattering inside the cavity region 1, the Green function has 
poles at the complex frequencies f2& = oj k — iF k /2 where the condition 



£>i(n fc ) - 1 + n 3 (Q k )e 2ie ^ n ^ c = 



(506) 



[7*13: generalised Fresnel reflection coefficient (I552H is met. We assume that the line widths T k 
are much smaller than the line separations T k <C (wk+i ~ Wfe_i)/2 so that the integration over 
frequency in the kernel can be restricted to subintervals A k = [(uJk-i + u k )/2, (u> k + ui k +i)/2]. 
Near a cavity resonance, the time integral in Eq. ( 1504b can then be performed in the Markov 
approximation to provide the cavity-induced shift of the atomic transition frequency as [51] 



5lo = 



UJ A UJ k ' 



\n k , 



2 _ u A r k > ln / ujj^ 

47T V U) A 



with 



4MI 2 . 2 

ak = 1 — Ti sm 
neoAl 



w k z A 



(507) 



(508) 



[u> A — lu a — 5uj]. Within these approximations, the kernel function K(t) for the excited-state 
probability amplitude C e (t) = C e (t)e~ lSu)t is then 



(509) 
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Fig. 32. The efficiency of single-photon Fock state preparation rj(t) (solid curve) and the excited-state 
probability \C e {t)\ 2 (dashed curve) for uik = 2 x W 8 Fk, (*k = 5 x lCF r Ffc, u>k — &A = 0.1Ffc,andthe 
radiative damping rate 7fc = 0.9Ffc. Taken from Ref. [51]. 

The solution for C e (t) is the prerequisite for the calculation of the photon extraction ef- 
ficiency r](t) from the cavity [51], which is defined as the probability of having the outgo- 
ing field prepared in a single-photon Fock state. The result for a typical situation is shown in 
Fig. [32] For very long times, a sufficiently high-Q cavity and almost exact resonance, the ex- 
traction efficiency approaches r/(t) — ► 7fcrad/Tfc. Here, the 7fc ra d is the radiative damping rate 
7 fcrad = c/(2l)\Tk\ 2 and Tk the transmission coefficient Tk = ti3(f2fc)e Mfci / c . The decay rate 
Tk of a cavity resonance is the sum of the radiative damping rate 7/crad and the absorptive damp- 
ing rate jkubs = c/ (2l)\Ak\ 2 where Ak is the absorption coefficient at the cavity resonance 
frequency uik- Hence, the extraction efficiency in the stationary limit is [154, 155] 

v ( t ) -» ^ . (510) 

7fcrad + 7fcabs 

For intermediate times, the efficiency shows an oscillating (tidal) behaviour. Thus, during each 
emission/re-absorption cycle within the cavity, the single-photon wave function builds up outside 
the cavity until the steady-state solution (empty cavity) is reached. 

6.2 Spherical microcavities 

Spherical microresonators possess a rich structure of field resonances to which an atom can be 
coupled when placed inside (whispering gallery modes) or outside (surface-guided modes) the 
resonator. 

6.2.1 Atom inside a spherical microcavity 

In cases in which the atom is embedded in a spherical microcavity, we can distinguish between 
large and small cavities, weak and strong coupling. In Sec. 14.11 we have discussed the case in 
which an atom was placed inside a very small cavity. This led us to the notion of local-field 
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Fig. 33. Scheme of the spherical cavity with outer radius Ri and inner radius ife- 




Fig. 34. Spontaneous decay rate T/To as a function of frequency (left panel). Excited-state probability 
I Ce (i ) | 2 as a function of decay time (right panel). For parameters, see text. Figures taken from Ref. [156]. 



corrections. We now have a brief look at a cavity whose radius is much larger than the relevant 
atomic transition frequency. In view of the three-layered structure depicted in Fig. [33] we thus 
require that R,2UJa/c ^> 1. Then, the rate of spontaneous decay can be calculated from Eq. (13421 > 
in the limit of thick cavity walls as [156] 

r p r> \ n ( UJ A) - it£m(R 2 uJA/c) ] 

r~r Re . -, , — Y~E> T\ ■ ( 511 ) 

The left panel in Fig.[34]depicts the normalised decay rate r/r as a function of frequency. We 
assume a single-resonance Drude-Lorentz model ( 1423b with ujp = 0.5ujt and 7 — 0.01o->t- The 
microresonator has parameters R2 = 30 Ay and Bi — R2 = Xt- One observes narrow -band 
enhancement (T/Tq ^> 1) alternating with broadband suppression (T/ro 1) of spontaneous 
decay. The narrow peaks are located at the cavity resonances. 

When the atomic transition frequency approaches one of the cavity resonances, the atom-field 
coupling becomes stronger. At resonance, when uja = where ujc 1S one of the resonance 
frequencies of the microcavity, the time evolution for the excited-state occupation probability is 
the one of a damped oscillator with damping constant Sue, the width of the cavity resonance. For 
small material absorption, 7 wj, ujp, uj p /ojt, the cavity resonance line can be approximated 
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by a Lorentzian with width 

s ^'n§kry (512) 

In the strong-coupling regime when the condition T(uic) ^> Sloc is fulfilled, the excited-state 
probability 



|C.(t)| 9 = e -*"°W J^^£ t ) ( 513) 



shows damped Rabi oscillations with Rabi frequency 57 = ^J2T(loc)8loc- The right panel 
in Fig. [34] shows the temporal evolution of the excited-state occupation probability for loa = 
1.046448w T , r A T /(2c) = 1CT 6 , and 7 = l(r 4 w T (solid line), 7 = 5 x l(r 4 w T (dashed line), 
and 7 = 10~ 3 ojt (dotted line). The other parameters are the same as those used for the curves in 
the left panel. For comparison, the exponential decay in free space is also shown (dashed-dotted 
line). 

6.2.2 Atom outside a spherical microcavity 

Atoms need not be located inside a closed dielectric or metallic structure in order to be coupled 
strongly to the electromagnetic field. Instead of whispering gallery modes that are due to total 
internal reflection inside a sphere, surface-guided modes can be excited on the outside surface of 
a spherical microcavity [157]. Assuming again a Drude-Lorentz model for the permittivity of 
the dielectric sphere that features a band gap between the transverse resonance frequency lot and 
the longitudinal resonance frequency ujl = \/w T + Wp. It is clear that for frequencies below 
lot (j^A < wr) whispering gallery modes are excited, whereas for frequencies within the band 
gap {lot < < ^l) surface-guided modes are excited. 

Figure[35]shows the rate of spontaneous decay of a radially oriented dipole very close (za = 
ta — R = O.lA^t) to a microsphere of radius R = 2Xa- In this near-field limit, the decay rate is 
approximately given by Eq. ( 13401 ) for a planar interface, 



T 1 - 31me(u A ) ( c 
TY ~ 4\ £ (lo a ) + 1\ 2 \loaza 



0(1). (514) 



The strong enhancement in the band-gap region (lot < < wl) is due to surface-guided 
modes. 

These surface excitations can be used to induce entanglement between two atoms located 
diametrically opposite around the microsphere [157, 159]. If the two atoms A and B are initially 
prepared in their respective excited and ground states, their combined initial state can be written 
in terms of the (entangled) coherent superpositions |±) = (|e^, gs) ± \gA, e s))/ which are 
initially equally excited. Associated with these superpositions are decay rates 



TV ~ 2^ Re | ik A r A ? K {kArA) 



1=1 



[ji{k A r A ) + rlfiLOA^ikArA)] [l T (-1)'] ) (515) 
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Fig. 35. Left panel: Decay rate as a function of the atomic transition frequency for a radially oscillating 
dipole near a microsphere of radius R — 2Xa- The chosen parameters are lop = 0.5i^t, 7 = W~ 4 lot, 
Ar — O.IXa- Figure taken from Ref. [158]. Right panel: Frequency dependence of the decay rates r± 
for two radially oriented dipoles at Ar^ = Ars = 0.02Aa- The other parameters are uop = 0.5u>t, 
7 = 10~ 6 oj T ,i? = 10A T [159]. 



[fcji = lua/c; ji(z) and /ij ^(z), spherical Bessel and Hankel functions; r^~p(u>A), generalised 
reflection coefficient given in Sec. lA.4.3l . If the atomic transition frequency uia coincides with 
one of the microsphere resonances I, the single-atom decay rate can be approximated by [158] 



r 



-Z(/ + l)(2/ + l)Re 



h?\k A r A ) 

k A TA 



(hi 
r P,i 



(516) 



and therefore the decay rates of the superposition states reduce to T± ~ T[l =p (— 1) This 
means that if I is even (odd), the state |— ) (|+)) decays much faster than the state |+) (|— )). The 
frequency dependence of T± is shown on the right panel in Fig. [35] Hence, depending on I one of 
the superposition states |±) is superradiant while the other is subradiant. After the superradiant 
combination has decayed, the two atoms reside in the subradiant entangled state |+) or |— ). 
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7 Outlook 

The theory of macroscopic quantum electrodynamics is a powerful tool that provides the link 
between isolated atomic ensembles and absorbing solid-state surfaces. It extends the well- 
established field of free-space quantum optics to new horizons. It is remarkable that, at least 
within the framework of causal linear response, one finds a Hamiltonian description of the 
medium-assisted electromagnetic field that does not have to fall back onto a master equation 
description. Instead, unitary evolution equations can be found that take into account absorp- 
tion processes as well. Examples of such unitary evolutions include the absorbing beam splitter 
(Sec. 13.2b and internal atomic dynamics (Sec. @J. In many cases, macroscopic quantum elec- 
trodynamics can provide more information than perturbative calculations, as exemplified by the 
theory of thermal Casimir-Polder forces (Sec. 15.5b . 

Obvious extensions towards nonlinear interaction Hamiltonians (Sec. 13.41 ) will be further 
studied to bridge the gap between classical descriptions of nonlinear processes to include both 
linear and nonlinear absorption, and to establish the connection with microscopic models of 
effective nonlinear interactions. Such a theory is likely to provide detailed information on the 
quantum state of light that emerges from, e.g. spontaneous parametric down-conversion crystals, 
and on the effect of higher-order nonlinearities. 

As it is evident from our collection of examples, for the time being we have always assumed 
that all macroscopic bodies as well as all atomic systems are at rest. Such a restriction is clearly 
not necessary. Centre-of-mass motion of atomic systems can be accounted for by either imposing 
a classical trajectory or more consistently by solving Newton's equations of motion in addition 
to Maxwell's equations. To assume a classical trajectory for an atom or a molecule means to 
neglect the backaction onto the motion which in many cases can be justified. 

For example, from our discussion in Sec.|5]it is clear that an atom moving parallel to a planar 
dielectric surface will not only experience a Casimir-Polder force in the direction perpendicular 
to its motion, but also in the direction of its motion. This effect is sometimes called quantum 
friction and can be understood in two ways. In our intuitive picture of Casimir-Polder forces 
being caused by forces between the atomic dipole and image dipoles in the magnetoelectric 
material, one can view quantum friction as a drag effect acting on the image dipoles due to finite 
conductivity or even resonances. Another way to view quantum friction is as an interaction with 
the surface plasmon at the magnetoelectric interface. 

Taking the quantum-mechanical character of atomic position and momentum into account 
will lead one to a theory of motional heating that is important in particular in ion traps. The 
opposite effect of cooling to the motional ground state, enhanced by coupling to a resonant 
microcavity, can also be understood within the framework of macroscopic quantum electrody- 
namics. 

Eventually, understanding how motion affects dispersion forces acting on isolated atomic 
systems will lead to a theory of a quantum theory of moving dielectric materials which thus 
far only exists for strictly non-absorbing materials. This in of fundamental interest as it makes 
contact with the low-energy limit of effective action in quantum field theory. 
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A Dyadic Green functions 

The notorious problem one perpetually faces in macroscopic QED is to find the Green tensor or 
dyadic Green function (DGF) associated with the classical electromagnetic scattering problem. 
This is not unusual as the DGF generalizes the usual mode expansion in free space to the solution 
of a more general boundary value problem. In what follows we will restrict ourselves to the 
special case of spatially local magnetoelectric materials. Extensions to materials with spatial 
dispersion can be found in Ref. [38]. 

We will need to look for the fundamental solution to the Helmholtz equation 



uj 2 

V x K(r,cj)V x T e(r,uj) 



E(r, to) = i(J, u)}(r, uj) (517) 



with the dielectric permittivity e(r, uj) and inverse magnetic permeability n(r, u>) = fj, 1 (r, uj) 
being arbitrary. The solution to Eq. d517t can be written in terms of the DGF as 

E(r, uj) = ifiouj J dV G(r, r', to) • j(r', to) (518) 

where the dyadic Green function satisfies 



uj 2 



V x k(t,cj)V x — -e{v,uj) 



c 



G{t,t\u>) = 8{t-t'). (519) 



Together with the boundary condition that G(r, T 1 ,u>) vanishes as |r — r'| — > oo, Eq. d519t 
has a unique solution provided that the strict inequalities Imefr, w) > and Im/i(r, uj) > 
hold. Physically, these requirements mean that all dielectric materials have to be passive, i.e. 
absorbing, media. 

A.l General properties 

The dyadic Green function inherits the analytical properties of the permittivity and permeability. 
That is, it is meromorphic in the upper complex frequency half -plane, and the Schwarz reflection 
principle 

G*(r,r',uj) = G(r,r',-uj*) (520) 
holds. The Onsager reciprocity theorem requires that the DGF satisfies the additional relation 

G(r',r,uj) = G T (r,r',Lo). (521) 
Particularly useful is the following integral relation involving products of Green functions, 

,2 



J d 3 s Im e(s, uj)G(r, s, uj) ■ G*(s,r',uj) 



+Im«(s,ai) |G(r,s,w) x vj • [V x G*{s,r',uj)] j = ImG(r,r',w) . (522) 

This relation is essentially the linear fluctuation-dissipation theorem, written in terms of the 
dyadic Green function. 
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A.2 Duality relations 

Another useful property of the DGF is its behaviour under a duality transformation. The dual 
DGF is defined by the equation 



V x K*(r,cj)V x T e*(r,Lu) 



G*(r,r',cj) = 5(r-r') 



(523) 



with e* = /i = 1/k, k* = l//i* = 1/e, i.e., it is the solution to the Helmholtz equation with 
e and /1 exchanged. By applying the duality transformation to the field expansions ( 1232b — ( f235l >. 
using the transformation rules listed in Tab. [2] and comparing coefficients, one can easily verify 
the following relations between the original DGF and its dual [49]: 



^G*(r,r» = 

V x G*(r,r',w) x V' = 

V x G*(r,r',w) = 

G*(r,r»xV' = 



-re(r,w)V x G(r,r',w) x V'k(v',uj) 
-«(r, uj)8(y — r'), 



(524) 



-e(r, uj) G(r, r', u)e{r',oj) + e(r, u)6(r - r'), (525) 



-e(r, w)G(r,r',w) x Vfi;(r',w), 
— k(t,lu)V x G(r, r', w)e(r', to) , 



(526) 
(527) 



In the next subsections, we will list the explicit forms of dyadic Green functions for important 
highly symmetric geometric arrangements of magnetoelectric bodies. 



A.3 Bulk material 

The simplest situation one can envisage is one of a homogeneous, isotropic dielectric medium. 
In this case, the dielectric permittivity e(lo) and magnetic permeability do not depend on 
the spatial position. This means that the magnetoelectric material and hence the dyadic Green 
function that describes it must be translationally invariant. Thus, the DGF can only depend on 
the difference between the coordinates q = r — r', G(r, r', lo) = G(g, u). We can therefore use 
Fourier transform techniques to solve Eq. d519t . Defining 

G(k,w) = J ^^G(0,u,)e- ik -s, (528) 

the Fourier transformed Helmholtz equation is 

uj 2 

-k x k x G(k, lo) - ^= r e(w)/x(w)G(k, uj) = u,(oj)I . (529) 

c z 

This is now an algebraic equation that can be easily solved. To do so, we decompose the unit 
tensor into its transverse and longitudinal parts with respect to the wave vector k, 

/ k®k\ k(g)k 
I= i 1 -— + — ■ (530) 
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Note that 

k x kx 



k 2 I 



k(g>k 
k 2 



(531) 



is transverse. Hence, the Green tensor is, after inserting into the inverse Fourier transform, 



G(q,lu) 



d 3 k 
(2^2 



MM 



k 2 c 2 



k® k 

k 2 



k® k 



u) 2 s(tu)fj,(uj) k 2 

(532) 



which can be transformed back into configuration space using contour integral methods. Noting 
that the dielectric permittivity s{lo) and the magnetic permeability u,(u)) have no poles or zeros 
in the upper complex frequency half -plane, the result is 



G(q,u>) 



iq{u>)g 



fj,((j) 



4irq 2 (ui)g 



(533) 



[q 2 (oj) — ^t£(ijj)ii(uj)]. After evaluating the derivatives one arrives at a decomposition into 
transverse and longitudinal parts as 



G" (<?><") = 
G ± (q,u) = 



Airq 2 
Airq 2 
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In particular, in free space where e(u>) = 1, we find that 

„2 



f 



\ug 



\ug 



(534) 



(535) 



(536) 



with f(x) = x + ix 2 — x 3 and g(x) = x + 3ix 2 — 3x 3 . 



A.4 Layered media: planar, cylindrical, spherical 

A more challenging situation arises if one considers magnetoelectric materials that are structured 
to form boundaries of certain symmetry. What we envisage here are examples of stratified ma- 
terials in planar, cylindrical and spherical geometries. The idea is to expand the dyadic Green 
function in terms of vector wave functions associated with the symmetry of the problem. The 
success of this procedure is directly related to the problem of separability of the Helmholtz op- 
erator V x /i _1 (r, w)Vx — %£(r,w). There are only a few coordinate systems in which this 
operator is separable [160], namely the cartesian, cylindrical, spherical and spheroidal coordi- 
nate systems. Hence, only there can the solution to the three-dimensional Helmholtz operator be 
reduced to seeking the scalar solution of a one-dimensional wave equation. 
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coordinate system 


pilot vector 


scalar eigenfunctions 


cartesian (x, y, z) 






cylindrical (r, ip, z) 


e 2 


ip n (k ry k Z ) r) = J n (k r r)e in ^' e tlizZ 


spherical (r, 6, <p) 




ipl, m (k,r) =j l (kr)P{ n {cosey m <<> 



Tab. 7. Pilot vectors and scalar eigenfunctions that are regular at the origin for cartesian, cylindrical and 
spherical coordinate systems. The functions J n (kr) and ji(kr) denote cylindrical and spherical Bessel 
functions, respectively. The P ; m (cos 9) are associated Legendre polynomials. 

The scalar Helmholtz equation takes the form 

(A + k 2 ) V'(k, r) = (537) 

[k 2 = ^e(r, u)]. Table|7]lists the scalar eigenfunctions that are regular at the origin for carte- 
sian, cylindrical and spherical coordinate systems. 

The scalar wave function ^(k, r) is then used to construct the irrotational (L) and solenoidal 
(M and N) eigenfunctions with respect to a pilot vector c as [161] 



L(k,r) = W(k,r), (538) 

M(k,r) = V x ft/>(k,r)c] , (539) 

N(k,r) = yV x V x [V>(k,r)c] . (540) 
k 

Because of the orthogonality of the scalar eigenfunctions, 

/ d 3 r?/>(k,r)V>(-k',r) = (27r) 3 <5(k - k') , (541) 



the vector wave function M, N and L are mutually orthogonal. 

Due to the Helmholtz theorem, these three types of eigenfunctions form a complete set of 
basis functions for the Helmholtz operator. The dyadic Green tensor, written in terms of these 
eigenfunctions, takes the form of a Fourier integral 

G(r, r', w) = J d 3 k [M(k, r) <g a(k) + N(k, r) b(k) + L(k, r) <g> c(k)] (542) 

with yet unknown expansion coefficients a(k), b(k) and c(k) which can be obtained using 
the mutual orthogonality of the vector wave functions. The interested reader is referred to the 
excellent textbook [161] for further details. 

For us, it suffices to know that these vector wave functions can be used to determine the (un- 
bounded) dyadic Green function in various coordinate systems. Our interest is now focussed onto 
some physically relevant situations in which two (or more) dielectric materials form a planarly, 
cylindrically or spherically layered structure. More specifically, if one were to consider only two 
adjoining materials, one would be looking at two half-spaces joined by a planar interface, an 
infinitely long wire or a dielectric sphere, respectively (see Fig.l36l>. 
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Fig. 36. Inhomogeneous dielectrics for which the dyadic Green function is analytically known. 

The unbounded DGFs for each individual dielectric material have to be joined at the common 
interface(s) across which the transverse field components of electric and magnetic fields have to 
be continuous. Let n be a unit vector normal to the interface r = r$ between two dielectrics. 
Then Maxwell's equations imply that at the interface 

n x G {fs) = n x G [u+1)s] , (543) 
— n x V x G (/s) = — !— n x V x G [(/+1)s] , (544) 

where the superscripts denote that the source point r' connects to the /ield point r via a Green 
function G^ s \ The dyadic Green function can then be decomposed into two parts, an un- 
bounded DGF C?( '(r, r', uj) that represents direct propagation from r' to r in an unbounded 
medium, and a scattering part G^ s \v, r', lu) that describes the contributions from multiply re- 
flected and transmitted waves, 

G (/S) (r,r» = G (0 ' /S) (r,r' ,u)8 fs + G (SJs) (r,r' ,u) . (545) 

Mathematically speaking, the scattering part of the DGF has to be included in order to fix the 
boundary conditions at the medium interfaces, whereas the DGF of the unbounded medium is 
responsible for the correct boundary conditions at infinity. 

A.4.1 Planarly layered media 

There are two distinct methods for finding the dyadic Green function for planarly layered media. 
Here we will describe the Weyl expansion method [41, 162] that gives very compact expressions 
for the DGF. The second method using vector wave functions [163] is somewhat less transparent 
and will be reserved for cylindrically and spherically layered media. 

The Weyl expansion is based on the translational invariance of the dyadic Green function 
with respect to the directions parallel to the planar interface, 

G(r,r» = J ^e ik r(P-p')^ ul ) G (k ll ,z,z',u } ) (546) 
[r = (p, z)]. The matrix components for the scattering part of the DGF (we drop the superscript 
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S for notational convenience) associated with reflection from the interface are given by [41, 162] 



G^{k {] ,Z,z',Lj) 
G { £ ] {\L h Z,z\Lu) 



2fci 



2fci 



Jk lz (\z\ + \z'\) 



ik lz (\g\ + \z'\) 



2fci 



■ r j( w )-BjfT- r »( w )y 



r p (w) 



2/r 



!_ e ifcl.(|2|+|2'|) 



-sgn(z')r p (w) J- 



(547) 

(548) 
(549) 
(550) 



where fcf = (w 2 /c 2 )£i(a;)/ii(a;) and fc 2 z = kf — A; 2 . The functions r s (u) and r p (w) are the 
Fresnel coefficients of s- and p-polarised waves, 



2r 



/i 2 (w)fci z - [i\(u)k 2z ' 



r p (w) 



e 2 (^)fci z - ei(o;)fc 2z 



(551) 



The remaining vector components can be deduced by employing the replacement rule Gj,^ = 

Got ^ (Mb fey) ar, d the reciprocity theorem G(r, r',u/) = G T (r',r,ix)) which translates into 
G( 11 )(k||,z,z',c) = G( 11 ) T (-k| h z',z, W ). 

The Fresnel reflection coefficients ( 15511 ) describe a single planar interface. For more than one 
interface, say a planar layer of thickness d, the reflection coefficients can be combined (for both 
polarisations) as 



ri2 



ri2 + r 23 e 



2ik,2z d 



1 - f2ir 23 e 



(552) 



where ry is the Fresnel coefficient for the interface between layers i and j. 

The matrix components of the scattering DGF associated with transmission through the pla- 
nar interface are 



GL 2) (k||,z,z',u; 
G^(k {] ,z,z',cu 

G& 2 >(k||,z,*> 
GW(k [b z,z',u 



ikiz\z\+ik 2 z\z'\ 



2fci 



. , ,k lz k 2z k 2 x , , , fcy 



2/r 



* gifcl*|z|+ifc2*|z'| 



2k lz 



2fci 



3 ifei z |z|+ife 2i |z'| 



,ifclz|z|+ifc2z|z'| 



? c iki z \z\+ik2z\z'\ 

2k lz 



sgn(z%(w) 
sgn(z')t p (uj) 

t p (uj) 



k\ z k x 



k\k 2 
k 2z k x 



k\k 2 



k\k 2 



(553) 

, (554) 
(555) 
(556) 

(557) 
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The transmission coefficients for s- and p-polarised waves are 

t s (u) 



M2M 



[1 + r a ((j)] , t p (u>) 



£2M 



[l + r„(w)] , 



(558) 



(12) (12) 

The remaining vector components can again by symmetry arguments such as G yx = G xy and 

Gyz — Gxz {J$x ~ {kx ^y) ^nd Cyy — Gxx {J^x 4r ^ &y)* 



A.4.2 Cylindrically layered media 

For the two examples that follow here, we use the expansion of the DGF in terms of vector wave 
functions. We first define the cylindrical vector wave functions [164, 165] that depend on the 
radial (77) and axial (ft) components of the wave number [fc 2 = e(Lo)/j,(Lo)ui 2 /c 2 = i] 2 + ft 2 ] as 



M. 



(ft) = 



n sin dZ n {r\g) cos 

T-Z n (rjQ) ntpe e . rupe^ 

g cos dg sin 



ihz 



Ah) = 



1 



dZ n (rjg) cos n sin 2 cos 

ih . rnpe e ^ih-Z n (j]g) nipe v + rj Z n (r]g) . nipe 2 

dg sin g cos sin 



(559) 

e lhz . 
(560) 



The trigonometric functions have to be chosen appropriately for the even and odd types of func- 
tions. In terms of those, the dyadic Green function for a cylindrical wire can be written in the 
form [164, 165] 



CO 

/CO 

n=0 



(2 - S n0 ) 



r% M Mf nm (ft) M yj m (-ft) + (ft) N ^ (-ft) 

+ 4N?!„ (ft) ® M'i^ (-^) + rS^Mi 1 ^ (ft) ® N'i^ (-*) 



r'(i) 



'(i) 



(561) 



OO 

/OO 

n— 



(2 - 5 n0 ) 



'(1) 

21 tvt. fu\ o.n/r'W f_M 1 j-21 A/f „ (M^lv'W 



t 21 

'A/M 



;m)2 (ft) ® Mi^ (-ft) + ^JvNe^ 



Wsn'.^ho 



+ 4VN>, 2 (ft) ® Mi^ (-ft) + ^M >7?2 (ft) (-ft) 



(562) 



00 

/CO 

n=0 



(2 - <S n0 ) 



% 2 



r(i) 



W ® M'e„„ 2 (-ft) + 4VNll, (ft) ® N'e„„ 2 (-ft) 



z 12 lvr 1 - 



+ t 12 M N^ (ft) ® M'. n)Ja (-ft) + t 12 N Mi l > (ft) ® N'. (-ft) , (563) 



(i) 



;nr; 2 1 
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G( 22 >(r,r» 



»M2 

8tt 



dh 



(2 — <S n p) 

n=0 



% 2 



nr)2 V 



(ft) ® M'e^-ft) + 4N !n , 2 (/,) ® N' !m)2 (-/i) 



+ r 



21 No^(h) ® M'ef-ft) + r^Mo (ft) 



NM 1 ^ e nr/ 2 < 



N'e n „ ,(-ft) 



(564) 



In this compact notation, a superscript means that the function Z n in the respective vector 
wave function is a Hankel function of the first kind, Hn \ Without that superscript, it is under- 
stood that the Bessel function J n is used. 

The reflection and transmission coefficients are determined by a 4 x 4-matrix equation 

-l 



T {H,v) 



■n 



(H,V) 



with 



Mi' 1 (r)j R) 
dR 


( ± C 3 H^( Vj R) 



i? 



OR 


dR 




OR 




dR 




PjJ n {rijR) 

dR 




± 



(ft 



R 

PjJ n (r)jR) 

T j dJ„(y j R) 
dR 




OR 


T jPj J n {r)jR) ) 



and the abbreviations 







ihn 



Pi 



With these preparations, the reflection and transmission coefficients can be derived as 

. -l 

r ii 
-11 

NM,MN 
.21 

MM,NN 
f 21 
NM.MN 



f 12 

MM,NN 
f 12 

NM.MN 



22 

r MM,NN 
22 

r NM,MN 




(565) 



(566) 



(567) 



(568) 

(569) 
(570) 
(571) 
(572) 



Note that, in contrast to planarly or spherically layered media, s- and p-polarised waves, repre- 
sented by the M and N vector wave functions, becomes mixed for all angular momenta n > 0. 
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A.4.3 Spherically layered media 



Similarly to the cylindrically layered case, we use the expansion of the dyadic Green function in 
terms of vector wave functions. The spherical vector wave functions are given by [166] 



,., m STvmt ^\ sm f 7 \ dPJ 71 (cos 0) cos 

M« ml fc = T^-^zi(kr )P l "(cose rrvp e - z t {kr) — - . mipe^, 

sin 6 cos rfB sin 



N. 



ml 



(k) = 



cos 



zAkr) P, m ( cos 9) mipe r 
kr sin 



+ 



£[(fer) 
kr 



dPI n ( cos 0) cos to „. , sin 

— ^ . m^eeT^-pr^ cose mtpe v 

dkd sin sm0 cos 



With the help of those, the dyadic Green function for a dielectric sphere is [166] 

oo I 



i[i\k 



2/ + 1 (Z - m 



G(21) (-» ^ ^EEE(-m|^§^ 



e,o i=l m=0 

r'M 



if M e m , (fc 2 ) ® M (fei ) + t^N >z (k 2 ) 



G< 12 )(r,r» 



l^ 2 k 2 
47T 



2/ + 1 (Z - m 



Z(Z + 1) (Z + m 



e.o 1=1 m— 

oo / 



i2i\/r( 1 ) cfc^^A/r'. rio _i_ +i 2 n1^ 

*M2^2 \^ \ ^ x ^ 2/ + 1 — m 



G (22) (r,r» = »EEE( 2 -^) 



Z(Z + 1) (? + m 



e.o /— 1 ??i—0 

x rf M. ml (i 2 ) ® M'e ml (k 2 ) + rf N >; (fc 2 ) ® N' >; (fc 2 ) 
For the following, it is convenient to define the Ricatti functions 



In terms of those, the Mie scattering coefficients read 

(ii) yihviih fyjiikiR) - tt2kir)i(k\R)ji(k2R) 



x—kiR 



„(n) 



fi 1 k 2 r 1 i(k 2 R)h ( l 1) (hR) - ^kitiihR^ifaR) 
Hik 2 rn(k 1 R)j l (k 2 R) - ^ 2 k 1 rji(k 2 R)ji(kiR) 

IHfobfafyjifaR) - M2fci-7/(fc 2 i?)^ (1) (kiR) 



(573) 



(574) 



(575) 

(576) 

(577) 

(578) 
(579) 

(580) 
(581) 
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r(22 ) = tHk^jkzR^jhR) - ^k^jhR^jhR) 



r(22 ) = viki&jhR^jhR) - /x 2 fci6(fc 2 fl)/i, (1) (fcifl) ( . 81) 
P ^ik 2 ^i(k 1 R)ji(k 2 R) - fi 2 k 1 r]i(k 2 R)h\ 1 \k 1 R) 



(12) = vMm^R^ihR) - ti{k 2 R)ji{k 2 R)] 

t(1 a) = ^ik 2 [^i{k 2 R)ji{k 2 R) - m{k2R)h { p {k 2 R)\ 
P (iM(k!R)ji(k 2 R) - wkiriifafyhP (kiR) 



,(21 



(584) 
(585) 

(586) 



t (2i) = KhfafafyMhR) - mihR^jhR)] 

P vMifaR^ifaR) - wkwifafyhP (hR) ' 

Note that the transmission coefficients t SjP can be further simplified by using the Wronski deter- 
minant between spherical Bessel and Hankel functions. 

A.5 Born series expansion 

In most situations, the geometric arrangement of dielectric bodies is not symmetrical enough 
to yield a separable Helmholtz operator in which case an expansion into vector wave functions 
would be feasible (Sec. |A.4| i. Instead, we will describe an iterative method known from quantum 
mechanics and quantum field theory as Born (or Dyson) series expansion of the dyadic Green 
function. This method applies to an arbitrary arrangement of dielectric bodies but, in general, 
converges quickly only if the dielectric contrast between bodies and the surrounding material is 
sufficiently small. 

Let us assume that from an arrangement of dielectric bodies, described by a dielectric per- 
mittivity e(r, lu), one can separate a part whose DGF (r, r', lj) is analytically known (e.g. 
vacuum, bulk material or layered media) and is described by a permittivity So(r, w) (see, for ex- 
ample, Fig. [3Tb. The dyadic Green functions are solutions to the respective Helmholtz equations 

V x V x G(r,r',w) - — e(r, w)G(r, r', w) = 6(r - r') , (588) 

& 

V x V x G (0) (r,r',w) - — £ (r,w)G (0) (r,r',tj) = 5(r - r') . (589) 

cr 

Subtracting both equations from one another, we find that the difference between both DGFs 
[scattering Green function G^ s '(r, r', oj)] solves the inhomogeneous Helmholtz equation 

LJ 2 



V x V x G (S) (r,r',w)- — e (r, w)G (s) (r,r',w 

& 

= ^Se(r,u) [G (0) (r,r» + G (5) (r,r»' 



(590) 
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e(r,co> 




^e (r,co)— 



Fig. 37. Dielectric body of arbitrary shape in front of a dielectric wall. The whole arrangement of bodies is 
described by a permittivity e(r,ui). The DGF for the half-space alone [with permittivity eo (r, u)] is known 
and can be used as starting point for the Bom series. 



where 5e(r, to) = e(r,u>) — eo(r,u>) is the perturbation from the permittivity £o(r, a;). Equa- 
tion d590t bears some resemblance to the Helmholtz equation (1517t for the electric field, with 
the current density i(ioUjj(r,ui) being replaced by the rhs of Eq. d590l l. Its formal solution is 
therefore 



(591) 



G^(r,r» = J d 3 sG ( - Q \r,s,u;)~Se(s,u) ^G^ (s, r', lu) + (s, r', u 

Because the unknown scattering DGF G^ (r, r', lS) appears on both sides of Eq. d59U . making 
it a Fredholm integral equation of the second kind, one can solve it iteratively as [144] 

G^(r,r» = ^ / d 3 s'G^(r,s',iu)-6s(s',cu)G^(s',r',Lu) 

J I d 3 s W'Gl )(r, S ».fe( S »Gl°)( S ', S » • Ss(s",w)G^(s",t',uj) 



2\ 2 

CO 



+ ... (592) 

Equation (1592b is known as the Born (or Dyson) series expansion of the scattering dyadic 
Green function G^fr, r', lu). It has a simple diagrammatic representation as shown in Fig. [38] 
The field propagates from the source point r' to the observation point r via the intermediate 
positions s^. Each of the arrows represents the dyadic Green function G^ (which we assumed 
to be analytically known), and each intermediate position (or vertex) carries a weight factor 
^fe(s^\o;). The intermediate positions are integrated over, and the number of those points 
increases with the order of the iteration. 

As the Born series ( 1592b is a perturbation series in fe(r, to), it is clear that fast convergence 
is only guaranteed for small enough permittivity (or equivalently, refractive index) contrast. It 
should be noted, however, that the series does always eventually converge to its unique solution 
G( s )(r,r',w). 
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r' r' r' 



Fig. 38. Diagrammatic representation of the Born series, Eq. l |592| l. Each arrow corresponds to the propaga- 
tion according to the Green tensor G^°\ each intermediate vertex at position s' 1 ' contributes with a weight 
factor i ^ T Se{s i - 1 \uj). The integration is over all intermediate positions s' ! '. 



A.6 Local-field corrected Green tensors 



While the DGF G(r,r',ui) connects the macroscopic electric field E(r) with a macroscopic 
source j(r') [recall Eq. ( I518H . we are often interested in the coupling of the local electromagnetic 
field to microscopic sources such as atoms. For atoms which are embedded in a magnetoelectric 
medium, this difference between microscopic and macroscopic quantities leads to local-field 
corrections which can be implemented via the real-cavity model: One assumes that both field 
point ri and source point r 2 are not situated directly in the medium but surrounded by small 
free-space cavities of radius R cav (cf. Fig. [39t. Modifying the original permittivity e(r, uj) and 
(inverse) permeability n(r, u>) describing the present media to 



eioc(r,w), «ioc(r,w) 



1 if |r - 

e(r, cj), k(t, uj) else, 



ri| < i?cav or |r - r 2 | < R c 



(593) 



the required local-field corrected DGF G\ oc can be found as the solution to 

, ,2 



V x Ki oc (r,w)V x j£i oc (r,t 



G(r,r',uj) = 6{r-r'). 



(594) 



The real-cavity model thus accounts for the fact that an atom occupies some space within the 
host medium, where the cavity radius is of the order of one half the lattice constant of the latter. 

For magnetoelectrics (but not for metals), the condition | yfejl\ R ca ,vUJ /c<C 1 is typically valid 
for the relevant frequencies of the electromagnetic field. In this case, the effect of the introduced 
cavities can be studied in a perturbative way by means of the spherical DGFs given in Sec. lA.4.3l 
One finds that the local-field corrected single-point DGF is given by [71, 167] 



MS), \ w /3(ei 



+ i 



- 1) c 3 
9einf 



(2ei 



1) 



9[e?(5A*i - 1) - 3ei - 1] 



5(2ei 
3ei 



I) 2 

2 



ujR c 



2ei 



1 



G (s) (ri,ri,o;) (595) 



where G is the uncorrected DGF, e$ = e(r i; uj) denotes the permittvity of the unperturbed host 
medium at the source and field points (similarly for /i) and = ^jEi\ii is the respective refractive 
index. The first term in Eq. ( 15951 ) represents contributions to the electric field which are reflected 
at the interior of the cavity and never reach the surrounding medium [type (i) in Fig. [39la)l. 
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Fig. 39. Real-cavity model for applying local-field corrections to (a) single- and (b) two-point DGFs. Three 
typical contributions (i)-(iii) to the corrected DGFs are schematically indicated. 



The second term represents contributions which are transmitted into the host medium and, after 
possible transmissions and reflections, transmitted back into the cavity [type (ii)]; each of the 
transmissions through the cavity surface gives rise to one factor in round brackets. Processes 
where the field is backrefiected from the outside of the cavity [type (iii)] have been neglected 
since they are of higher order in the small parameter \y/eJI\R cav uj / c. Similarly, the local-field 
corrected two-point DGF is found to be [71] 

3s 3co 

Gioc(ri,r 2 ,cj) = — - — G(ri,r 2 ,w) forri^r 2 . (596) 

lE\ + 1 2e 2 + 1 

Equation ( 1596b represents the transmission of the field out of the first cavity, transmissions and 
reflections within the host medium, followed by a transmission into the second cavity [type (ii) 
of Fig.[39tb)l; each of the transmissions at the cavity surfaces gives rise to one of the two factors. 
Again, reflections at the outside of either cavity have been neglected [type (iii)]. 

When describing magnetic interactions, the curl of the DGF typically arises. Local-field 
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corrections in such situations cannot simply be obtained by taking the curl of Eqs. ( 15951 ) and 
( 15961 ) above, since the implementation of the local-field correction via the real-cavity model does 
not commute with this vector operation. Instead, the local-field correction has to applied after 
taking the curl. One finds the single-point DGF [58] 



VxG< s )( ri ,r llW )xV'| k 



A*i-1 



3 /i?(5ei-l)-3/xi-l 



2ttc 2 1 2/xi + l w 3 i?3 av 5 (2/xi + l) 2 ujR c 



+ i 



(2/xi + l) 5 



2/ii + l 



VxG^fri.r^JxV (597) 



and the two-point DGFs [58] 



V x G(ri,r 2 , to) 
G(r 1 ,r 2 ,Lu) x V' 

V x G(ri,r 2 , lj) x V' 



3e 2 



loc 2//i + 1 2e 2 + 1 
3ei 3 



loc 2ei + 1 2^ 2 + 1 

3 3e 2 



2/xx + 1 2e 2 + 1 



G(ri,r 2 ,w) forri 7^ r 2 , 
G(ri,r 2 ,w) forri ^ r 2 , 

G(ri,r 2 ,w) forri 7^ r 2- 



(598) 
(599) 
(600) 



The behaviour of the local-field corrected Green tensors under a duality transformation (cf. 
Sec. IA.2l i can be easily derived by combining Eqs. J595t -( l598] l with Eqs. ( 1524l )-( l52"7i i: 



VxG%,r 1(W )x V'| loc 

V x G*(n,r 2 ,w) x V'| loc 
V x G*(r,r',w)L 

v 3 1 / Hoc 



-VxG( s )(n,r llW )x V'l, , 



■7 2 - G ioc( r i> r i> w )- 



(601) 
(602) 



-VxG(n,r 2 ,w) x V'| loc forn ^ r 2 , (603) 



= — 5- Gioc^^ra,^) forri ^ r 2 , 
cr 

= -G(ri,r 2 ,w)x V'| loc forri^r 2 , 
G*(r,r',w) X V'lioc = -VxG(ri,r 2) w)| loc forn^r 2 . 



(604) 

(605) 
(606) 
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